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ABSTRACT 


In  the  continuing  search  for  materials  which  can  withstand  the  grueling 
requirements  of  modem  day  applications,  Functionally  Graded  Materials  (FGMs)  seem  to 
be  a  promising  alternative  to  conventional  materials.  These  nonhomogeneous  materials 
offer  better  interfacial  properties  by  improving  bond  strength  and  reducing  thermal 
mismatch.  Before  putting  these  materials  into  application,  an  important  step  in  the  design 
of  FGMs  is  the  stress  analysis  and  fracture  characterization.  The  fracture  performance  of 
FGM  coatings  on  homogeneous  substrates  is  the  focus  of  this  study. 

In  this  study,  various  internal  and  surface  crack  configurations  in  the  coating  and 
the  substrate  are  subjected  to  mechanical  and  thermal  loads.  The  analysis  is  linear  elastic. 
The  thermo-mechanical  properties  of  the  FGM  coating  are  assumed  to  vary  exponentially 
with  the  spatial  coordinate.  The  equilibrium  equations  are  solved  using  integral 
transforms.  The  resulting  singular  integral  equations  are  solved  using  numerical 
integration.  The  results  of  interest  for  this  mode  I  formulation  are  the  stress  intensity 
factors  and  the  crack  opening  displacements.  The  effects  of  the  nonhomogeneity 
parameter  and  various  dimensionless  length  parameters  are  studied. 

One  of  the  most  important  outcomes  of  this  study  is  the  theoretical  proof  that 
"kink"  in  material  property  at  the  interface  does  not  introduce  any  singularity.  In  the 
numerical  results  it  is  observed  that  generally  the  stress  intensity  factors  tend  to  increase 
with  material  nonhomogeneity.  Also,  it  is  observed  that  the  substrate  thickness  tends  to 
suppress  cracking  in  the  coating.  In  pure  thermal  loading,  the  surface  cracks  may  either  be 
arrested  or  there  might  be  crack  closure.  The  stress  intensity  factors  from  different 
loadings  can  be  added  up  to  obtain  the  resultant  stress  intensity  factor  for  multiple  loading. 

Results  in  this  study  have  wide-ranging  applications.  They  can  be  applied  to 
thermal  barrier  coatings  on  turbine  components,  combustion  chambers,  parts  of  the 
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airframe  for  the  "Space  Plane",  soil  mechanics,  bone  fractures  and  many  more  applications 
where  the  material  is  macroscopically  nonhomogeneous.  Thus  this  study  solves  a  basic 
problem  common  to  a  variety  of  applications  in  diverse  fields. 
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CHAPTER  1 


INTRODUCTION 


1.1  Introduction 

In  modem  engineering  applications,  the  requirements  on  materials  are  becoming 
more  and  more  strict.  The  materials  are  subjected  to  grueling  conditions.  Conventional 
materials  do  not  stand  up  to  the  high  temperatures  or  corrosive  environments.  Hence  there 
is  always  a  need  for  new  materials  or  new  material  systems.  Most  failure  processes  such  as 
corrosion,  wear  and  fatigue  are  surface  related.  Hence  these  failures  can  be  controlled  by 
controlling  the  material  properties  near  and  at  the  surface.  Mere  surface  treatments  such  as 
surface  hardening  or  introduction  of  residual  compressive  stresses  to  prevent  fatigue  crack 
initiation  may  not  be  sufficient  in  many  applications.  In  such  applications  it  is  necessary 
either  to  use  new  materials  to  withstand  the  conditions  or  to  provide  the  bulk  material 
with  a  protective  coating  of  a  more  resistant  material.  Examples  of  such  applications 
include  the  blades  in  gas  turbine  engines  which  are  subjected  to  high  stresses  and  highly 
corrosive  environment,  machine  tools  undergoing  excessive  wear  and  combustion 
chambers  which  have  to  sustain  very  high  temperatures  upto  3000F.  In  all  these 
applications  it  is  necessary  to  develop  new  systems  where  conventional  materials  have 
failed. 

Various  attempts  have  been  made  in  developing  new  materials  for  these 
applications!!].  The  most  promising  ones  are  monolithic  structural  ceramics  and  fiber 
reinforced  ceramics  including  carbon-carbon  composites.  The  disadvantage  of  these 
ceramics  is  their  brittle  nature.  This  can  be  improved  by  introduction  of  fibers  or  whiskers 
to  avoid  catastrophic  cracks.  Carbon-carbon  composites  can  sustain  very  high 
temperatures  but  have  very  low  oxidation  resistance  and  extraordinarily  long 
manufacturing  time. 
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All  the  new  materials  mentioned  are  essentially  homogeneous.  In  today's 
demanding  technological  environment,  these  materials  seem  to  be  inadequate[2].  The  main 
requirement  is  for  materials  to  possess  good  thermomechanical  as  well  as  structural 
properties.  In  high  temperature  applications,  the  material  is  required  to  have  high  heat  and 
corrosion  resistance  as  well  as  high  mechanical  toughness  and  heat  conductivity.  Similarly 
in  gears  and  bearings,  a  high  degree  of  wear  resistance  and  high  toughness  is  required  in 
the  same  material.  Homogeneous  materials  seem  to  fail  in  one  regard  or  another.  Hence 
the  idea  of  coating  the  substrate  with  a  heat-resistant  ceramic  layer,  called  the  thermal 
barrier  coating,  appears  very  attractive.  Here  is  where  composites  seem  to  find  widespread 
acceptance  in  industry. 

One  of  the  earliest  applications  of  composites  may  have  been  in  protection  against 
corrosion  by  using  a  variety  of  coatings.  Composite  materials  are  designed  to  utilize  the 
mechanical  properties  of  the  component  materials  so  as  to  optimize  certain  physical  bulk 
properties  of  the  medium.  In  thermal  barrier  coatings,  a  ceramic  coating  is  applied  to  the 
base  substrate[3, 4].  Usually  the  substrate  is  some  Ni  based  superalloy.  The  coating  may 
be  of  Yittria-stabilized  Zirconia  on  a  bond  coat  of  NiCoCrAlY.  Zirconia  makes  the 
component  heat  resistant  while  the  bond  coat  protects  the  substrate  from  oxidation.  The 
mechanical  strength  is  derived  from  the  substrate.  The  disadvantage  of  such  systems  is  the 
high  thermal  and  residual  stresses  and  relatively  poor  bonding  strength. 

A  novel  way  of  reducing  the  residual  stresses  and  improving  the  bond  strength 
between  the  coating  and  substrate  is  to  layer  the  interfacial  zone  going  from  metal  rich  to 
ceramic  rich  compositions.  For  example,  in  joining  tungsten  to  zirconia  by  introducing 
layers  that  contain  80/20, 60/40,  40/60  and  20/80  percent  W/Z1O2,  it  was  shown  that  the 
peak  value  of  the  residual  stress  becomes  approximately  one-sixth  of  that  obtained  from 
direct  W-Zr02bonding[5].  The  next  logical  step  is  the  processing  of  fully  tailored 
materials  and  interfacial  zones  with  predetermined  continuously  varying  volume  fractions. 
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Such  materials  are  referred  to  as  "  Functionally  Graded  Materials"  (FGM).  These  are 
essentially  nano-composites  or  fine  composites  having  continuously  varying  rather  than 
piecewise  constant  volume  fractions.  Apart  from  increasing  the  bond  strength  [6]  and 
reducing  the  thermal  stresses,  grading  also  seems  to  improve  the  fracture  toughness  and 
fatigue  and  corrosion  crack  growth  resistance  parameters  of  thermal  barrier  coatings.  With 
the  advent  of  new  techniques  in  powder  metallurgy  and  in  surface  chemistry  it  is  indeed 
possible  to  manufacture  fully  tailored  functionally  graded  materials  with  varying  volume 
fractions. 

The  philosophy  of  FGMs  is  different  from  any  of  the  conventional  materials[7]. 
FGMs  are  designed  by  an  "Inverse  Design  Procedure"  [8, 9]  in  which  the  material  is 
tailored  to  suit  the  needs  of  any  application.  Various  ingenuous  techniques  are  used  in  the 
synthesis  of  FGMs  [10].  Vapor  deposition  techniques[ll],  centrifugal  casting[12], 
combustion  sintering[13, 14],  plasma  spraying[15]  are  a  few  of  them.  A  good  review  of 
different  manufacturing  techniques  can  be  found  in  [16].  Chemical  and  physical  vapor 
deposition  techniques  can  produce  very  clean  FGMs  but  the  processes  are  very  slow  for 
industrial  applications.  Plasma  spraying  is  by  far  the  most  popular  technique  because  of  the 
relatively  simple  equipment  and  high  deposition  rate.  The  disadvantage  of  plasma  spraying 
is  the  introduction  of  porosity  at  the  substrate/coating  boundary.  Also,  some  oxide  layers 
can  form  at  the  splat  boundaries. 

As  the  manufacturing  of  these  materials  can  lead  to  some  inherent  flaws,  it  is 
important  to  carry  out  the  fatigue  and  fracture  characterization  of  FGMs.  This  requires  the 
solution  of  some  standard  crack  problems.  Most  flaws  are  surface  related  and  under 
working  conditions  these  flaws  may  give  rise  to  surface  cracks  which  can  run  through  the 
component.  Hence  it  is  important  to  study  the  surface  cracks  and  their  growth  under 
various  loading  conditions.  This  is  the  main  problem  investigated  in  this  study. 


5 


Applications  of  this  study  are  wide-ranging.  The  FGM  layer  may  be  the  thermal 
barrier  coating  of  a  combustion  chamber.  It  could  be  the  thermal  coating  for 
microelectronic  chips.  A  recent  application  is  that  of  a  "Space  Plane"  [17, 18].  The  engine 
components  and  parts  of  the  airframe  of  such  a  plane  will  be  subjected  to  very  high 
temperatures  and  FGMs  are  proposed  as  a  possible  material  system.  Nonhomogeneous 
materials  with  varying  material  properties  also  occur  in  various  non-industrial  applications. 
In  Biomechanics,  the  cement  line  in  secondary  bones  can  be  considered  to  be  an  FGM.  In 
soil  mechanics  and  rock  fracture  analysis,  the  results  from  this  study  could  be  readily 
applicable.  Thermomechanical  material  nonhomogeneity  may  also  be  observed  in 
components  having  highly  temperature  dependent  properties  and  subjected  to  steep 
temperature  profiles.  The  diffusion  layer  between  two  dissimilar  materials  is  usually  non- 
homogeneous[19].  Thus  this  study  is  attempting  to  solve  a  basic  problem  common  to  a 
variety  of  diverse  fields. 

1.2  Literature  Survey 

Literature  on  non  homogenous  materials  can  be  classified  into  3  periods.  The  early 
literature  of  50s  and  60s  is  mostly  about  nonhomogenous  infinite  or  semi-infinite  planes. 

In  the  70s,  emphasis  was  on  layered  structures.  Here  individual  layers  were  homogenous. 
Interest  in  nonhomogeneous  layers  seems  to  have  picked  up  in  the  mid-80s  and  is  being 
vigorously  pursued  in  present  day  research.  Most  of  the  research  in  FGMs  is  from  this 
third  period. 

The  early  extension  of  the  theory  of  elasticity  to  non-homogeneous  bodies  can  be 
found  in  [20].  The  early  research  of  non-homogeneity  in  elasticity  and  plasticity  can  be 
found  in  the  articles  in  [21].  Most  of  the  articles  derive  equations  for  non-homogenous 
infinite  or  semi-infinite  planes.  Some  solutions  for  displacements  and  stresses  are  found  for 
specific  forms  of  elastic  modulus. [22] 
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In  the  70s,  the  focus  was  on  layered  structures.  With  the  advent  of  composites,  a 
great  deal  of  research  was  concentrated  on  ceramic  coatings  for  metal  substrates.  These 
so-called  cermets  are  precursors  to  the  present  day  FGMs  and  hence  are  important  for  the 
purpose  of  comparison  of  performance  with  FGMs.  The  fracture  characterization  of 
layered  structures  can  be  found  in  many  publications[for  exg,  23  -  32].  Interface 
flaws[23  —  27],  flaws  perpendicular  to  the  interface[27  —  31]  and  delamination[32] 
problems  were  extensively  studied.  The  interest  in  layered  structures  in  still  continuing  but 
is  gradually  being  replaced  by  graded  layers. 

FGMs  were  looked  into  vigorously  in  the  90s.  Most  of  the  literature  pertains  to  the 
design  and  manufacturing  of  FGMs.  A  good  review  of  different  manufacturing  techniques 
is  given  in  [16].  As  these  materials  are  macroscopically  non-homogeneous,  material 
characterization  becomes  an  important  issue.  Thermomechanical  characterization  has  been 
attempted  with  various  different  approaches.  A  review  of  the  techniques  on  bounds  of 
effective  properties  is  given  by  Torquato[33].  Among  these  techniques  is  the  variational 
approach  developed  largely  by  Hashin  and  Strickman[34].  The  bounds  on  effective 
properties  have  been  improved  by  various  techniques  such  as  simulation  techniques,  cell 
models  and  experimental  techniques.  Reader  is  referred  to  [33]  for  more  references  on 
each  of  these  techniques. 

There  is  still  a  considerable  ambiguity  regarding  variation  of  properties  within  an 
FGM  layer.  The  correlation  between  composition  and  properties  is  being  studied 
extensively.  Experiments  using  say  indentation  techniques  to  determine  hardness  show 
wide  variations  in  results  depending  on  the  point  of  indentation[3].  Hence  for  analysis 
some  suitable  assumption  has  to  be  made  for  the  variation  of  material  properties.  Another 
restriction  on  the  analysis  is  that  the  governing  equations  for  nonhomogeneous  materials 
tend  to  be  partial  differential  equations  with  variable  coefficients.  This  has  been  explained 
in  detail  in  Chapter  2.  There  is  no  known  technique  for  solving  such  equations  in  a  direct 
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way.  The  best  approach  is  given  by  Varley  et  al[35].  This  approach,  though,  is  suitable 
only  for  a  limited  class  of  variations.  From  the  foregoing  discussion,  it  is  obvious  that  the 
problem  is  analytically  tractable  only  for  specific  variation  of  the  material  properties.  A 
linear  variation  of  material  properties  was  considered  in  [36, 37].  Bakirtas[38]  solved  the 
problem  of  a  rigid  punch  on  a  nonhomogeneous  medium  assuming  an  exponential 
variation.Exponential  variation  seems  to  be  a  more  attractive  assumption  as  the  material 
property  never  goes  to  zero  as  in  the  linear  case.  This  assumption  is  used  in  this  study. 

Before  FGMs  can  be  put  to  use  in  practical  applications  it  is  important  to  study  the 
fracture  properties  of  these  materials  to  predict  the  service  life  of  the  components.  Not 
much  analytical  work  has  been  done  in  fracture  mechanics  of  FGMs.  Introduction  of  a 
crack  in  an  infinite  space  was  done  by  Konda  and  Erdogan[39].  The  problem  of  an 
interface  crack  was  solved  by  Delale  and  Erdogan[40].  Chen  and  Erdogan  solved  the 
interface  crack  problem  for  a  coating[41].  Crack  perpendicular  to  interface  was  first 
studied  by  Erdogan  for  Mode  HI  loading[42].  Same  problem  was  solved  by  Erdogan, 

Kaya  and  Joseph  [43]  for  Mode  I  loading.  Erdogan  et  al  [44]  also  solved  the  problem  for  a 
nonhomogeneous  diffusion  layer  between  homogeneous  half  spaces.  Erdogan  and  Ozturk 
[45]  have  studied  the  axisymmetric  interfacial  crack  problem.  The  edge  crack  problem  for 
an  FGM  layer  was  looked  into  by  Wu  and  Erdogan  [46].  It  can  be  seen  that  the  more 
practical  problem  of  a  crack  perpendicular  to  the  interface  in  an  FGM  coating  on  a 
substrate  is  yet  to  be  solved  for  Mode  I  loading.  This  is  the  problem  that  is  solved  in  this 
study. 

Another  important  part  of  this  investigation  is  the  thermal  loading  and  resulting 
thermal  stresses  in  FGMs.  This  aspect  has  been  studied  in  literature  to  some  extent. 
Thermal  stress  analysis  for  layered  cylindrical  shells  was  done  by  Miller  et  al[47]  while  that 
for  heterogeneous  shells  was  done  by  Venkatarmana[48].  Numerical[49]  and 
experimental[50]  evaluation  of  thermal  shock  and  thermal  stresses  in  FGM  plates  have 
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been  studied.  Analytical  techniques  have  mainly  been  limited  to  beam  theory  sort  of 
approaches  [51, 52].  The  elasticity  problem  of  an  FGM  layer  under  thermal  loading  was 
looked  into  by  Wu  and  Erdogan[53].  Though,  in  most  analyses,  the  temperature  field  has 
been  separated  from  the  displacement  field,  a  coupled  thermoelastic  model  has  been 
developed  by  Zhang  et  al[54].  In  this  study,  the  temperature  and  displacement  fields  are 
uncoupled.  Variations  in  both  the  modulus  of  elasticity  and  the  coefficient  of  thermal 
expansion  are  considered  to  be  exponential. 

1.3  Outline  of  the  dissertation 

Chapter  1  of  this  study  deals  with  the  introduction  of  the  problem,  and  the  related 
literature  survey. 

Chapter  2  contains  the  statement  of  the  problem  and  its  formulation.  The 
governing  equations  have  been  derived  and  boundary  conditions  stated. 

In  Chapter  3  the  singular  integral  equations  are  derived.  For  different  cases,  the 
singularities  and  the  stress  intensity  factors  are  derived. 

Chapter  4  contains  the  numerical  procedure  to  solve  the  singular  integral 
equations.  For  various  cases,  the  numerical  procedures  differ.  These  have  been  covered  in 
this  chapter. 

Chapter  5  gives  the  results  for  some  specific  material  properties.  These  results 
have  been  discussed  and  compared  with  previously  published  results. 

In  Chapter  6,  conclusions  drawn  from  this  study  have  been  listed  and  directions  for 
future  research  are  given. 

The  appendices  A-F  present  some  analytical  details  and  material  properties 
(Appendix  E). 
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CHAPTER  2 

FORMULATION  OF  THE  PROBLEM 


2.1  Geometry  of  the  problem 

The  geometry  of  the  plane  elasticity  problem  described  in  the  introduction  is 
shown  in  Fig  2.1. 


The  functionally  graded  coating  lies  to  the  left  of  the  y-axis  while  the 
homogeneous  substrate  lies  to  the  right.  This  figure  shows  the  most  general  configuration. 
There  can  be  various  different  configurations  depending  on  the  values  of  £>1,01,02  and 
.When  i>i  7^  aiand  <22  =  £>2.  we  obtain  an  internal  crack  in  the  nonhomogeneous  coating 
only.  When  b\  =  ajand  02  7^  62  > obtain  an  internal  crack  in  the  homogeneous 
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substrate  only.  When  both  b\  ^  oiand  02  ±  62 » we  obtain  internal  cracks  in  both  the 
coating  as  well  as  in  substrate.  When  a<i  =  62,  61  =  h\  and  0  <  a\  <  h\,  we  obtain  an 
edge  crack  in  the  nonhomogeneous  coating.  When  b\  =  h\  and  a\  =  0,  =  0  and 

0  <  62  <  h%,  we  obtain  an  edge  crack  crossing  the  interface.  Another  case  would  be  when 
0  <  &i  <  hi,  a\  =  0,  a2  =  0  and  0  <  62  <  ^2-  In  this  case,  there  is  an  internal  crack 
crossing  the  interface.  All  these  cases  have  been  considered  in  this  study. 

2.2  Solution  procedure 

The  solution  procedure  involves  finding  the  governing  equations  for  the 
nonhomogeneous  part  as  well  as  those  for  the  homogeneous  part.  These  governing 
equations  are  then  solved  to  obtain  displacements  subject  to  the  boundary  conditions  at 
the  boundaries  and  the  interface. 

The  displacements  for  the  cracked  nonhomogeneous  strip  and  the  cracked 
homogeneous  strip  are  determined  by  a  superimposition  technique.  For  example,  the 
superimposition  technique  for  the  nonhomogeneous  strip  is  depicted  in  Fig  2.2. 


Problem  1 

Figure  2.2  Superposition  for  the  nonhomogeneous  part 


Problem  2 
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The  strip  is  considered  to  be  a  combination  of  a  cracked  infinite  space  and  an  uncracked 
strip  with  appropriate  boundary  conditions  applied  at  the  boundaries.  This  is  not  the 
conventional  superposition  method  as  the  geometries  of  the  two  problems  are  different. 
But  the  logic  here  is  that,  as  long  as  the  governing  equations  are  satisfied  within  a  domain 
and  the  boundary  conditions  are  satisfied  along  the  boundaries,  the  solution  is  valid.  Hence 
the  problem  reduces  to  finding  the  solution  of  two  simpler  problems  and  adding  up  the 
solutions.  Similar  procedure  is  applied  for  the  homogeneous  strip. 


2.3  Governing  equations 

The  governing  equations  for  the  elasticity  problem  are  the  equilibrium  equations. 
For  the  homogeneous  part  the  development  is  straightforward.  For  the  nonhomogeneous 
part,  though,  the  shear  modulus  is  a  function  of  x.  Hence  certain  assumptions  need  to  be 
made  to  make  the  problem  analytically  tractable. 


2.3.1  Navier's  equations  for  nonhomogeneous  elastic  medium 


The  equilibrium  equations  for  a  two-dimensional  elasticity  problem  can  be 
expressed  as 


XX  .  dTx y 

dx  dy 


=  0, 


(2.1) 


dTsy  ,  d<3yy  = 
dx  dy 


(2.2) 


The  relation  between  stresses  and  displacements  is  given  by  Hooke's  law  which  can 
be  stated  as 


&XX 


Ms) 

K  —  1 


du  dv 

{(3-k) -+(!+«)->, 


(2-3) 
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(2.4) 


'yy 


ll(x)  du  dv 

£i<(3_K)_  +  (1+K)_}, 


Txy 


.  .  ,  du  dv 

* x){ai  +  ai} 


(2.5) 


where  fi(x)  is  the  shear  modulus  and  «  is  given  as 


K 


3-4i/  for  plane  strain 


3-t/ 

l+i/ 


for  generalized  plane  stress 


(2.6) 


Substituting  (2.3  —  2.5)  into  (2.1)  and  (2.2)  and  assuming  that  k  is  constant  we  obtain 
Navier's  equations  for  a  two-dimensional  elasticity  problem, 


d2u 


Mi+K)  gj  +  M(*-l )gj  +  l»gxdy  ■  dx 


d2u  _  d2u  du  du 
—  +  ^(1+^ 


du  dv 
+  =  °' 


d2v  d2v  d2u  du,  dv 


d/j  .  „ .  du 

+  —  (k— 1)^-  =  0. 
dx  dy 


(2.7) 


(2.8) 


It  can  be  seen  that  equations  (2.7)  and  (2.8)  are  coupled  partial  differential  equations  of 
the  second  order  with  variable  coefficients  as  is  an  arbitrary  function  of  x.  There  is  no 
direct  analytical  way  to  solve  these  differential  equations.  As  discussed  in  the  introduction, 
attempts  have  been  made  in  the  past  to  solve  such  equations  using  specific  functions  for 
In  this  study  it  has  been  assumed  that  the  shear  modulus  has  an  exponential 
behavior,  namely. 
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Kx)=  we' 


for  —h\  <  x  <  0 


(2.9) 


Cte) 


where  is  the  shear  modulus  of  the  homogeneous  substrate  and  /?  is  the 
nonhomogeneity  parameter  defined  as 


hi  m 

Substituting  (2.9)  into  (2.7)  and  (2.8)  we  obtain, 


d2u  d2u  d2v  du  dv 

(1+Kte  +  +  2g^  +  =  °- 


d2v  d2v  d2u  dv  du 

(K-lfe  +  (1+">  0?  +  2dxdy  +  +  =  °' 


(2.10) 


(2.11) 

(2.12) 


Equations  (2.11)  and  (2.12)  are  partial  differential  equations  with  constant  coefficients 
and  can  be  solved  analytically. 


2.3.2  Navier's  equations  for  homogeneous  elastic  medium 

These  equations  can  be  immediately  derived  from  equations  (2.11)  and  (2.12)  by 
letting  /5  go  to  zero.  The  equations  can  then  be  obtained  as 


d2u  d2u  d2v 

(1+/c)^Zo  +  +  2-r . t.  =  0, 


'dx  2 


dy1  dxdy 


(2.13) 


/  1N32«  ...  ,d2v  d2u 

<K_1)a?  +  <1+KV  +  2ela~y  =  °- 


(2.14) 


2.3.3  Superposition  technique 

As  discussed  earlier  the  problem  is  split  into  two  problems,  one  with  a  cracked 
infinite  plane  and  other  with  uncracked  strip.  The  solutions  for  each  of  these  problems  are 
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obtained  in  this  section.  The  nonhomogeneous  and  homogeneous  parts  are  considered 
separately. 

2. 3. 3.1  Nonhomogeneous  part 

The  superimposition  technique  is  shown  in  Fig  2.2.  For  problem  1  ( cracked 
infinite  plane)  Fourier  transform  with  respect  to  x  is  the  natural  choice.  For  problem  2 
(uncracked  strip),  it  should  be  noted  that  u  is  symmetric  about  the  x-axis  while  v  is 
antisymmetric  about  the  x-axis.  Hence  Cosine  transform  is  chosen  for  u  while  Sine 
transform  is  chosen  for  v.  In  the  solutions,  the  subscript  c  denotes  cracked  infinite  plane 
while  subscript  uc  denotes  uncracked  strip. 

Solution  of  Problem  1 

Taking  Fourier  transform  with  respect  to  x,  we  can  write  the  displacements  as 


1  f°° 

uic  =  —J^  My,*)  t^da,  (2.15) 

i  r°° 

Vic  =  9iiy,*)  e^d a.  (2.16) 

Substituting  into  equation  (2.11)  and  (2.12)  and  taking  Fourier  inversion  we  obtain 

(K-I)f![-(K+I)(a2+Pia)fi  +  ( (3(3-K)-2ia  )g[  =  0,  (2.17) 

(K+l)g'{-(K-l)(oi2+/3ia)g1  +  ( /3(K-l)-2ia  )f[  =  0.  (2.18) 


These  are  coupled  second  order  differential  equations.  Hence  there  should  be  four 
arbitrary  constants.  The  solutions  for  /land  g\  can  then  be  written  as 
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(2.19) 


4 

9i(y,a)  =  ^2D]<a)  enjy, 

3-1 
4 

/i(y,a)  =  yimjDjCe)  enjV.  (2.20) 

j=  i 


Substituting  (2.19)  and  (2.20)  into  (2.17)  we  obtain 
4  f 

T:  {  (k— l)mjDj(a)nj  -  (k+  1  )(a2 +(3ia)mjDj(ci) 

+  ( f3(3-K)-2ia  )Dj(a)rij  j  en-'y  =  0.  (2.21) 


Since  enjys  are  linearly  independent,  equation  (2.21)  should  go  to  zero  for  each  value  of 
j  =  1, 2, 3, 4.  Hence  we  obtain  mj  as 


( 2ia—/3(3—K)  )rij 

fffl  *  ■—-*  ■  .  .  .  .i.  —  I.  ** 

J  (K-l)rvj  -  (k+1  )(a2+/3ia)’ 


(2.22) 


On  the  other  hand  substituting  (2.19)  and  (2.20)  into  (2.18)  gives, 

4  f 

Y2  {  («+l )Dj(a)n2j  -  (n—l)(Q?+flia)Dj(oi) 
j- 1 

+  ( (3(K-l)-2ia  )m3Dj{a)nj  }  en^  =  0.  (2.23) 


Again  using  the  linear  independency  condition,  we  obtain, 


mj  = 


~(K+l)n2j  +  (K-l)(a2+{3ia) 

(  /3(k— l)—2ia  )rij  ’ 


j  =  1,2, 3, 4. 


(2.24) 


Equating  (2.22)  and  (2.24)  we  obtain  the  characteristic  equation  as 
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(2.25) 


(rtj— (a?+/3ia)j  —  62nj  =  0 
where 


This  gives  4  solutions  of  n3  as 
n\  =  -  "  ^  y/S2+  4(a2  +/3ia), 

n<i  =  ~ 1  y/iP+  4(a2+/3ia), 
nz  =  ^4“  V^2+  4(a2+/?ia), 
n4  =  —+^2+  4(a2+/?ia). 


(2.26) 


(2.27) 


(2.28) 


(2.29) 

(2.30) 


As  the  problem  is  symmetric  about  the  2 -axis,  we  consider  only  the  upper  half  plane.  Now 
for  u  and  v  to  be  bounded  as  y  goes  to  00,  n\  and  712 are  the  only  acceptable  solutions. 

The  expressions  for  displacements  and  stresses  can  then  be  given  as 


1  r°°  ,  ^ 

uic  =  — y  YjnjDjia)  en’yeTlQXd.a, 


(2.31) 


V\c 


i  poo  2 

=  2=/ 

J- »  j.i 


(2.32) 
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°xx\ C  =  K_1  2^  J  22^-(K+l)mjiQ+(3-K)njJDj(a)cTljy&~iaxda, 


(2.33) 


line 1  /*°°  f  1 

Vyyic  =  ~— j-  ^  J  J2 { -@-K)mjia+(K+ 1  )nj jDj(a)en^ye  taxda. 


(2.34) 


i  roo  2 

T'xylc  =  ,oe^~J  ^  |-io+mJnJ|l)7(Qf)enjye~'iQ:Eda. 


(2.35) 


Solution  of  Problem  2 

Considering  the  symmetries  of  u  and  u,  we  take  Cosine  transform  of  u  with 
respect  to  y  and  Sine  transform  of  v  with  respect  to  y. 

2  f°° 

Uhlc=  -J  f2(x,a)cos(ay)da,  (2.36) 

2  f°° 

vi uc  =  -  J  92(x,a)  sin(ay)da.  (2.37) 

Substituting  into  (2.11)  and  (2.12)  and  taking  the  respective  inverse  transforms  we  obtain, 

(K+D/2  +  /5(«+l)/2  -  («-l)a2/2  +  2a&+  /3(3—K)a92  =  0,  (2.38) 

(/c-l)p2+  /?(«-l)^2  ~  («+l)a2/2  -  2a/2~  /3(/c-l)a/2  =  0.  (2.39) 

These  are  coupled  second  order  differential  equations.  Hence  there  will  be  4  arbitrary 
constants.  The  solutions  can  be  written  as 
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(2.40) 


9i(x,a)  =  V.4,(q)  ep’x, 

J=1 

4 

/2(x,a)  =  e***. 

>-1 


(2.41) 


Substituting  (2.40)  and  (2.41)  into  (2.38)  we  obtain, 

4  , 

J2  {  (K+1)qjPj  -  («-l)a2?i  +  2 apj 
i=i 

+  (3 ((3-/c)a  +  (K+l)<jjp^  |  Aj(a)  e^x  =  0.  (2-42) 


Since  ep>xs  are  linearly  independent,  equation  (2.42)  should  go  to  zero  for  each  value  of 
j  =  1, 2, 3, 4.  Hence  we  obtain  qj  as 


_  —(3 —k)/3q  —  2  apj 

qj  ( K+ 1  )(pj  +  /3pj)~  (k-  1  )Qf2  ‘ 


(2.43) 


On  the  other  hand  substituting  (2.40)  and  (2.41)  into  (2.39)  gives 

E  {  (*-l )Pj  ~  (K+l)a2  -  2 aqjpj 

+  P(K-l)(pj-  oiq^j  e?ix  =  0. 

Using  the  linear  independency  condition  we  obtain, 

_  (/c-l)(p^  +  ppj ;)-  (/c+l)a2 

^  2apj+(K-l)fia 

Equating  (2.43)  and  (2.45)  we  obtain  the  characteristic  equation  as 


(2.44) 


(2.45) 
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|(p|  +  Sp})  -  a2)  +lra2  -  0 


(2.46) 


where  6  is  defined  in  (2.26).  The  4  solutions  for  pj  are  as  follows 


Pi  =  —  ^  — /32+4q2+4z<5q, 


(2.47) 


Pi 


—  p-  —  ^\//52+4a2— ^ 4i6a, 


2  2 


(2.48) 


pz=  —  ^+-\//?2+4a2+4i6a!, 


(2.49) 


P4 


=  —  ~+~  \/  /52+4a2— 4z(5q. 


2  2 


(2.50) 


The  displacements  and  stresses  can  then  be  expressed  as 


2 

=  -  /  yqjAj(a)  ePjXcos(Qy)da, 

71  Jo 


(2.51) 


2  /°°A 

^luc  =  -  >  ,A;(a)  ep’Isin(ay)dQ, 

^0 


(2.52) 


~^f0  V.  { (K+ 1  )g;P7+(3-^)Q}^7(Q)&RZCOS(Q;y)dQ, 


Vyyiuc  =  ^  J  {(3“«)9jPj+(«+  l)a} Ai(a)ep^cos(ay)dQ!, 


(2.53) 


(2.54) 
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2  4  r  i 

Txyiuc  =  ftoe{f5x)  -J^  {pj~qjajAj(a)eKXsm(ay)da. 


(2.55) 


By  superimposing  these  two  sets  of  solutions,  we  obtain  the  displacements  and  stresses  in 
the  nonhomogeneous  cracked  strip  as  follows: 


-h  E 

*  n  J  — 00  ^1 


m,jDj(a )  enjye~^axda  + 


2  rA 

-j  y^gjAjCa)  ePjZcos(ay)dQ> 


(2.56) 


i  /»oo  2 

vi  =  —  Y^Dj(a)  enjyt~iQxda.  + 
^  J  —00  ^1 


o  /*oo  4 

—J  y^.Aj(a)  e^sinCaj/jda, 


(2.57) 


1  Z*00  v~-\  f  1 

— — —  —  V  ~(K+l)mjia+(3-K)rij  [Dj^e^e^da  + 

K  J-  wj-oo  j=l  ^  ' 

~JQ  ^{(«+1)9jPj+(3-«)a}Aj(a)eRIcos(ay)da, 


(2.58) 


-_1  2^  J  ^{-(3~«)^i*a+(K+l)rij  |£)j(a)en^e_taa:da  + 

^k-1  ~kJq  S  “  Kte+(«+ 1  )a  }  A  J(a)ePjXcos(ay)da , 


(2.59) 
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\  f°°  2  f  N 

Txyl  =  ^e0x)  —  J  y^| -ia+mjnj  |£>j(a)enjye-,ardof  + 

2  f°°  4  r  i 

fj,oe^x)-  /  S']  <j  Pj—qjCx  MJ(a)ep’:rsin(Qy)da.  (2.60) 

*Jo  1  J 


2. 3. 3. 2  Homogeneous  part 

The  superposition  in  the  homogeneous  part  is  same  as  in  the  nonhomogeneous  part 
and  is  shown  in  Figure  2.3.  Again  the  subscripts  c  and  uc  denote  cracked  infinite  plane  and 
uncracked  strip,  respectively. 


Problem  3  Problem  4 

Figure  2.3  Superposition  for  the  homogeneous  part 


Solution  of  Problem  3 

Taking  Fourier  transform  with  respect  to  x 


u*‘  =  h. 

poo 

/  fs(y,a)  t~iQXda, 

(2.61) 

poo 

/  93(y,a)  e'^da. 

*  -oo 

(2.62) 

Substituting  into  equation  (2.13)  and  (2.14)  and  taking  Fourier  inversion  we  obtain 
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(k-1)/s  -  («+l)(«2)/3  -  =  0, 

(«+l)P3  “  («-l)(a2)^3  -  2*0/3  =  0. 


(2.63) 

(2-64) 


These  are  coupled  second  order  differential  equations.  By,  for  example,  eliminating  <73  we 
can  obtain  a  single  4th  order  differential  equation  in  fo.  The  equation  obtained  is 

/f  -  2a 2/£  +  a4/s  =  0.  (2.65) 

If  we  assume 

h  =  (2.66) 

we  find  that  m  has  double  roots  at  ±  a.  Hence  the  solution  of  fs  can  be  written  as 

/s  =  C\e~\a\y  +  C2ye-My  +  C3e +  CiyeH* .  (2.67) 

Considering  the  symmetry  of  the  problem  about  the  x-axis,  we  solve  it  only  for  the  upper 
half  plane.  The  displacements  need  to  be  bounded  as  y  goes  to  00.  This  necessitates  C3 
and  C4  to  be  identically  zero.  Hence  the  solution  can  then  be  written  as 


h  =  +  C2ye~  Wv,  (2.68) 

93  =  -i  (c^-My  +  C^e-l^+y)).  (2.69) 

The  displacements  and  stresses  are  then  given  as 

u2c  =  ^  J°°  (Cre-Wy  +  C2ye~^ e^d a,  (2.70) 
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^2c  =  ~  (Cie-Wy  +  C2e-^y(^+y)y-iQXda,  (2.71) 

&xx2c  =  t^J  {-2iaCi  +  iC2  ((3-/c)J^-2ay)  je-'^e^da,  (2.72) 

Vyy2c  =  ^J  {2 iaCi  +  iC2 ((l+«)J|!+2oy)  (2.73) 

7*00 

rxy2c  =  7^J  {-2|a|Ci-C'2(2|a|y+K-l)}e-^e-iQXda.  (2.74) 


Solution  of  Problem  4 


Considering  the  symmetries  of  u  and  v,  we  take  Cosine  transform  of  u  with 
respect  to  y  and  Sine  transform  of  v  with  respect  to  y  and  obtain, 


2 

^2 xtc  =  -f  /4(x,a)  cos(ay)da,  (2.75) 

Try  0 


2  z*00 

^2w  =  -J  94(?,a)  sin(ay)da.  (2.76) 


Substituting  into  (2.13)  and  (2.14)  and  taking  the  respective  inverse  transforms  we  obtain, 

(*+l)/4  “  («“  1)q2/4  +  2 ag'4  =  0,  (2.77) 

(«-l )fl4  “  («+l)«2/4  -  20/4  =  0.  (2.78) 


These  are  coupled  second  order  differential  equations.  By,  for  example,  eliminating  we 
can  express  this  as  one  fourth  order  differential  equation  in  g 4.  The  equation  thus 
obtained  is 
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g'l"  -  2a2^l  +  a4g4  =  0. 


(2.79) 


If  we  assume 


94  = 


(2.80) 


we  find  that  n  has  double  roots  at  ±  a.  Hence  the  solution  of  <74 and  /4can  be  written  as 


g4  =  +  B2xeax  +  B3e +  B^e^ , 


(2.81) 


f4  =  —Bie0^  —  B2(x--)eQX  +  B3e_QI  +  B4(x+-)e"ax. 


(2.82) 


The  displacements  and  stresses  can  then  be  written  as 


=  -f{- 

*J0  1 


Bxe°*  -  B2(x - )eQX 

a 


+  B3e  0:1  +  B4(x+—)e  QXlcos(ay)da, 
a  ) 


(2.83) 


V2uc  =  -J  jjBie01  +  B^xe01  +  B3e  01  +  B4xe  QX|sin(o;y)da  (2.84) 


aX3i2vc  =  ~^JQ  {laBxe**  +  B2(2ax-K-l)eQX 

+2aB3e~QX  +  S4(2a:c+K+l)e_^  jcos(a:y)da!,  (2.85) 

(Tygiuc  =  ^-JQ  +  B2(2ax-K+3)eQX 

+2aB3e~ax  +  B4(2ax+K— 3)e-Q::r  jcos(a:y)da:,  (2.86) 
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Tzyivc  =  ~JQ  {2a5ieax  +  S2( 2ax-K+l)eQX 

— 2aBse~ax  —  B^{2ax+K—  l)e~^  jsin(o:y)da.  (2.87) 

By  superimposing  these  two  sets  of  solutions,  we  obtain  the  displacements  and  stresses  in 
the  homogeneous  cracked  strip  as 

t*  =  If  (cie-My  +  C2ye~ M*)  e^da  + 


in- 


Bic"5-  Bo(x--)ec 
a 

K 


+  B$e  ^  +  B4(x+—)e  jcos(ay)da, 


(2.88) 


V2 


•  /»oo 

=  ~J  (Cie-Nv  +  C2e-l^(^+y))e^da.  + 


—J  ^Bie01  +  B2xeCiX  +  B$e  01  +  ^xe  QX|sin(ay)da, 


(2.89) 


<7xx2  =  ^ J  {-2iaCj  +  iC2((3-K)^-2ayj'je  Mye-*a*d a 

J  feaB^01  +  B2(2ax-«-l)eQX 

+2aS3e_QX  +  jB4(2ax+«+l)e~^  jcos(ay)dai,  (2.90) 

ayy2  =  ^  J  ^2icxC\  +  iC2  ((l+/c)-^)+2ay)  t~'lox  &a  + 

^ J  {2 aSie“+  S2(2ax-K+3)e£“ 

+2aS3e-Q:;r  +  B4(2ax+K-3)e-^  jcos(ay)da,  (2.91) 
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Txy2  =  gy'0O{-2|a|C1-C72(2|a|y+K-l)  }e"|Q|i/e^da  + 

**/  ^aSie01  +  B2(2ax— K+l)eQX 

-2 aBae-0*-  B4(2ax+K-l)e~QX^sm(ay)da.  (2.92) 


It  should  be  noted  that  in  equations  (2.56  -  2.60)  and  in  equations  (2.88  —  2.92)  there 
are  totally  12  arbitrary  constants.  These  will  be  evaluated  with  help  of  boundary 
conditions. 


2.4  Boundary  conditions 

The  homogeneous  boundary  conditions  are  given  on  the  free  surfaces  and  the 
interface.  On  the  boundaries,  the  stresses  are  zero  while  on  the  interface  stresses  and 
displacements  are  continuous.  These  can  be  expressed  as 


&xxi  ( —  y)  —  0,  (2.93) 

Txyi(-hi,y)  =  0,  (2.94) 

°xxl  (0,  y)  =  a xx 2  (0,  y),  (2.95) 

Trj/l  (0,  y)  =  'Txy2  (0,  y),  (2.96) 

«i(0,y)  =  1*2(0,  y),  (2.97) 

ui(0,y)  =  u2(0,y),  (2.98) 

(h2,y)  =  0,  (2.99) 

Tx&(h2,y)  =  0.  (2.100) 

On  the  plane  of  the  crack  we  have  symmetry  conditions  given  by 

Txyi  (x,  0)  =  0,  -  h\<  x  <  0  (2.101) 
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(2.102) 


Txy2  (x,  0)  =  0,  0  <  x  <  h2 

and  the  mixed  boundary  conditions  given  by 


<?yy\  (X,  0)  =  Pl(x), 

—  bi<  x  <  —  a\ 

(2.103) 

ui(x,0)  =  0, 

—  h\<  x  <  —  b\  ,  —  a\<  x  <  0 

(2.104) 

Vyi/2  (x,0)  =  P2(x), 

<22  <  X  <b2 

(2.105) 

v2(x,0)  =  0, 

0<  x  <a2  ,  &2<  x  <h2. 

(2.106) 

The  p\(x)  and  P2(x)  are  the  negative  of  the  stresses  at  the  crack  plane  under  external 
loading  in  an  uncracked  specimen.  The  mixed  boundary  conditions  will  give  rise  to  the 
singular  integral  equations  which  will  be  derived  in  Chapter  3. 
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CHAPTER  3 


SINGULAR  INTEGRAL  EQUATIONS  AND 
STRESS  INTENSITY  FACTORS 


In  Chapter  2  the  governing  equations  and  the  expressions  for  stresses  and 
displacements  were  derived.  The  boundary  conditions  were  also  stated.  In  this  chapter,  we 
will  use  those  expressions  to  derive  the  singular  integral  equations.  Also,  the  singularities 
and  the  stress  intensity  factors  at  the  crack  tips  will  be  discussed.  In  deriving  these 
expressions,  the  variable  x  is  replaced  by  x\  and  x%  defined  as, 


( xi,  x  <  0 

(  X2,  X  >  0 


(3.1) 


This  is  useful  in  distinguishing  the  non-homogeneous  part  from  the  homogeneous  and 
hence  the  two  singular  integral  equations. 


3.1  Development  of  singular  integral  equations 

3.1.1  Conditions  along  plane  of  symmetry  y  =  0 

Using  equations  (2.101)  and  (2.60),  it  may  be  seen  that 


1  f°°  2  r  i 

Moe(/3ll)—  /  y^<-ia+mjnj>Dj(a)enjye~ZQXl  da  + 

J-* °  >1  y— >0 


2  f00  4  r  'i 

/jLoe(PXl)—  /  ^  j  pj—qjQ.  \  Aj(a)e?jXl  sin(ay)da 

'Kj— OO  ^  ' 

Taking  now  the  Fourier  inverse  transform,  we  have, 


=  0. 


y— >0 


Di(a) 


(mgn 2  -  ia ) 
(mini  —  ia)  2 


(3.2) 


(3.3) 
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Similarly  from  equations  (2.102)  and  2.92),  it  follows  that, 


/£0 

2tT. 


y’OC{-2|a|Ci-C2(2|a|y+K-l)}e-lQ|j/e^QI1da  |  + 

"X  y—>0 


2^o 
7 r 


J 


+  i?2(2ax— «+l)eQ21 


—2aBze  ^  —  I?4(2a:x+ft— l)e  aXl  |sin(a!y)da!  |  =  0. 

y~*0 

By  taking  the  Fourier  inverse  transform,  we  have, 


(3.4) 


_  .  .  1  —  ft 

c,(a)  =w  2(a)' 


(3.5) 


We  now  define  the  following  new  unknown  functions 


0i(*i)  =  |^(*i»0) 

and 

02(^2)  =  f>,0) 


-  b\  <  xi  <  -  ai 


a2  <  x2  <  62. 


(3.6) 


(3.7) 


By  taking  Fourier  inverse  transform  of  (2.57),  (2.89)  and  using  (3.3),  (3.5),  (3.6),  (3.7), 
we  obtain 


£>2  (a) 
and 

C2(a) 


(mini  —  ia)  i 
(mini  —  m2n2)  a 


yi(t)e*°*dt 


2 

ft  + 1 


dt. 


(3.8) 


(3.9) 
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3.1.2  Homogeneous  boundary  conditions 

By  using  the  eight  homogeneous  boundary  conditions  (2.93  —  2.100)  and  the 
expressions  from  (2.56  —  2.60)  and  (2.88  —  2.92),  after  taking  the  appropriate  Sine  or 
Cosine  transforms  we  can  write  a  system  of  eight  algebraic  equations  in  eight  unknowns  . 
The  system  of  equations  we  obtain  is  as  follows: 


tfiAie-?1'11  +  + 


d3A3e^hl  +  64A4e~^hi 

II 

eO 

1 — 1 

(3.10) 

TiAie~Ph'  +  r2  A2e~^  + 

r3A3e_P3/l1  +  T4A4e~pAh 

=  R%, 

(3.11) 

<f>lA-[  +  <f>2A2  +  4>sA3  +  4>4A4  + 

2o;5i  —  (/c  +  1)^2  +  2aB3  +  (k  +  1)^4 

=  Rh 

(3.12) 

riAi  +  r2^2  +  t3a3  +  t4a4  — 

2olB\  +  (k  —  1)^2  +  2olB3  +  (/c  —  l)j?4 

=  R\, 

(3.13) 

91^1  +  92^2  +  93-^3  +  94-A4  + 

b4-b2--b3-b4- 

OL  CL 

=  Rb 

(3.14) 

A\  +  A2  +  A3  +  A4  —  B\  —  B3 

=  R&, 

(3.15) 

2aB\eCthl  +  (2ah2  —  k  —  l)B2eaf 12  + 

2aB3e~0‘hl  +  (2ah2  +  k  +  l)B4e~ahfi 

II 

ba 

(3.16) 
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Rl  (3.17) 


2ctBieah*  +  (2 ah2  -  k  +  1  )B2eah*  - 

2 aB3e~ah*  +  (2 ah2  +  k  -  l)B4e~ah2 

The  expressions  for  the  variables  used  in  these  equations  are  given  in  Appendix  A.  It 
should  be  noted  that  while  evaluating  the  right  hand  sides,  some  straightforward  contour 
integrals  have  been  evaluated.  The  system  of  equations  in  (3.10  -  3.17)  can  now  be 
solved  for  the  unknowns  4Z  and  Bl  (i  =  1, 2, 3, 4)  giving, 


Bi 

=  *oi  +  ei*3  +  e2  Ba, 

(3.18) 

—  *02  +  /i*3  +  h  Ba, 

(3.19) 

b3 

=  *03  +  -^171  +  -4-272  +  -4373  +  4-474, 

(3.20) 

Ba 

=  Rqa  +  A\8\  +  42^2  +  4363  +  A48a, 

(3.21) 

Ai 

=  +  A3lo3  +  A4u>a, 

(3.22) 

Ai 

=  Rj2  +  43P3  +  44P4, 

(3.23) 
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II 

k) 

CO* 

(3.24) 

Aa 

=  *14- 

(3.25) 

Again  the  expressions  for  the  various  variables  that  appear  in  (3.18  -  3.25)  are  given  in 
Appendix  A. 

3.1.3  Mixed  Boundary  conditions 
The  First  Singular  Integral  Equation 

Using  condition  (2.103)  and  expression  in  (2.59)  we  obtain 
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c p(Pxi)  i  r°° 

-  -  —  J  Y\-(3- K,)mjia+(K+ 1  )rij j D/ope^e"*0^1 da  + 


£?;f£{  (3  -  K)qjpj+(K + 1  )a  |  Aj(a)e?jXl  cos(ay)da  | 


piM- 


Substituting  from  (3.3)  and  (3.22-3.25),  from  (3.26)  we  obtain, 


1  r°° 

± 

2tT  J-oo 


D2da  -) —  /  R^\da  =  p\(x\) 


K  —  1 

/ioePXl 


(3.26) 


(3.27) 


where  $1  and  R^1  are  given  in  Appendix  A.  Using  (3.8)  we  can  write, 


1  /s  /*— ai  9  r — 

—  kie^'da  dtx  +  -  /  R2la(a,xi,ti)da 

v  — oo  i/  — bj  ^  Jo  J  — b\ 

2  f°°  fa  K  —  1 

+  -/  R2ib(a,x1,t2)da  g2(t2)dt2  =  pi(xr) — (3.28) 

7T  Jo  /^oePl 

where  the  function  k\  and  R2\a,  R2\b  are  given  in  Appendix  A.  R2\a,  Ri\b  are  bounded. 

k\,  however,  needs  to  be  investigated  further.  Consider  the  first  integral  in  (3.28).  It  can 

be  written  as. 


i  /*~nGi  r°°  ^ 

b  9\{t\)dtx J  k1(a,y)eia(tl~Xi')da. 


(3.29) 


Separating  ki(a,  y)  into  its  real  and  imaginary  parts  and  noting  that 
k\(  —  a,  y)  =  complex  conjugate  of  k\ (a,  y)  , 

(3.29)  may  be  expressed  as 


(3.30) 
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“  J ^  9i(ti)dtiJ^  |i?e(fei(a,y))cos(a(fi -xi)) 


-  Im(ici(a,y))sm(a(ti  -xi))|da  (3.31) 

In  (3.31)  the  dominant  part  of  the  kernel  may  be  separated  by  taking  the  asymptotic 
expansion  of  k\  (a,  y )  as  a  tends  to  infinity.  We  find  that, 


lim  k\(a,y)  =  kloo 

a— *oo 
— 

where  k  is  given  as 


(3.32) 


k  loo 

where 

/N 

k\ooh 


l/^\3/«-l\2  i//5\4(«-  l)(/c2  +  2/c  -  5) 

2  U/  Vk  +  1/  +  4  \a/  (k  +  l)3 


(3.33) 


i  / /?\  5  (k  -  1)(k2  +  6/t  -  13)  J_ 

**"  8  \a)  (k  +  1)3  + 


(3.34) 


Subtracting  the  asymptotic  limit  in  (3.32  -  3.34)  from  real  and  imaginary  parts  of 
k\(a,  y),  (3.29)  can  be  written  as. 
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4  /At  —  1\  I*-*1 

TT\K  +  l)  J-bi 


(t\  -  *l) 


dt\ 


+ 


b  9i(ti)dti  |iJe(fei(a,y))cos(a(ti  -xi)) 

-  ^Jm(fci(a,y))+4^^-^))sin(a:(ti  -®i))}da 


+ 


lj '  'gxMdhf  {(jJe(fc,(  a,y)  -  i?e(fciooA(a,y)))cos(a(t1  -  xi)) 
-  (jm(k1(a,y))  +  4(^pj)  ~ 

/m(fciooA(a,  y)jj  sin(o;(ti  -  x1  ))}<Ja 


+ 


^Jb  9\{h)dt\  ^  {iJe(fcioo^(a,y))c  os(a(fi  -  xi)) 

-  Im  (kioohia,  y))sm(a(t{  -  xi))  j-da. 


(3.35) 


Thus  we  have  separated  the  Cauchy  singularity.  Using  (3.35)  and  putting  all  the  bounded 
parts  together,  (3.28)  may  be  expressed  as 

-  r  7f-—\dt' + -  si  (*i)^u  (^=i. 

^  J —b\  (^1  ^1/  ^  J  —b\ 

+  ^J  92(t2)Ku(xi,t2)dt2  =  pi(xi)^ .  (3.36) 

The  expressions  for  K\\  (xi,  1 1)  and  K-&  (xi,  t2)  are  given  in  Appendix  A.  In  (3.36),  the 
first  integral  on  the  left  has  the  Cauchy  kernel  whereas  the  other  two  integrals  have 
Fredholm  kernels. 


35 


The  Second  Singular  Integral  Equation 

Using  equations  (2.105)  and  (2.91),  we  find, 


AfO 

27T 


J  $y2iotC\  +  iC<2  ^(l+K)-L^+2ayj  je  ^ve  tQ:r2d a  + 

y— >0 


iri  { 2aBi‘ 


e**2  +  B2(2ax2-K+Z)eaXi 


+2aB^e  ^  +  B4(2ax2+K-3)e  012  jcos(ay)da|  —  P2(x2)  (3.37) 

y— >0 


Also  from  (3.5)  and  (3.18-3.25)  we  obtain, 


If00  2  f°° 

2iLnc*ia+«L  ^ ia 


P2(x 2) 
1^0 


(3.38) 


where  and  R^  are  given  in  Appendix  A.  By  substituting  from  (3.9),  (3.38)  can  be 
written  as, 


1  f°°  _  fa  o  r°°  /■— <21 

—  /  k2e~tQX2da  /  g2(t2)e^dt2  +-  R2la{a,x2,ti)da  /  gi{tx)dtx 

*nJ-oo  Jaz  T\  J 0 


+ 


“/  R2ib(oc,X2,t2)da  f  g2(t2)dt2= 
nJ  0  Jos  /io 


(3.39) 


where  k2  is  given  in  Appendix  A.  Consider  the  first  integral  in  (3.39).  It  can  be  written  as, 


j  roo 

— J  gM2)dt, J  h,(a,  yie^'-^dda . 


(3.40) 


After  some  manipulations,  the  dominant  part  of  (3.40)  can  be  written  as 
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(3.41) 


7TVAC  +  1  /  Joq  (t%  —  X2) 

Substituting  (3.41)  back  into  (3.39)  and  rearranging  we  have, 

~f  7T^  ~\dt2  +  ~  [  9l(ti)K2l(x2,ti)dti 
n  Ja 2  \h  -  X2)  7T  J_bl 

K  +  1 

4mo  ' 


1 

+  -  92(h)K22(^2,h)dt2  =  Pifa) 

K  J„n 


(3.42) 


Expressions  for  K21  fa,  ti)  and  K2zfa,  <2)  are  given  in  Appendix  A.  In  (3.42),  the  first 
integral  has  the  Cauchy  kernel  while  the  other  two  integrals  have  Fredholm  kernels. 


3.2  Asymptotic  analysis  of  Fredholm  kernels 

The  Fredholm  kernels  in  (3.36)  and  (3.42)  may  contain  generalized  Cauchy  kernels 
which  are  singular  and  play  an  important  role  in  determination  of  the  singularity  in  certain 
crack  configurations.  It  is  hence  necessary  to  perform  asymptotic  analysis  on  the 
integrands  of  the  Fredholm  kernels  as  a  tends  to  infinity. 


The  expression  for  K\\  (a,  x\,  t{)  is  given  in  Appendix  A.  For  easy  reference  it  has 
been  reproduced  here. 
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Kn(oc,x1,t1)  =  (4 {j?e(fci(a,y))cos(o:(*i -*1)) 

-  (jm(fei(a,y))+4(^|))sin(a(ii  -*i))}da 


+ 


{  (i2e(fci(a,y)  -  i?e(fci00/l(o:,y)))cos(a(i1  -  *1)) 
-  (jm(ft(a,y))+4(^p|)  - 

Im(kicoh(®i y))J sin(a(f  1  -a?i))}da 


+ 


r  {Re^  00  h  (a,y))cos(a(<i  -  *1)) 

-  Im  (kioohia, y))sm(a(t\  -  *i))}da 


roc 

+  2/  ^210(0:^1. 

Jo 

It  should  be  noted  that  the  first  three  integrals  do  not  lead  to  any  singular  term.  The  terms 
in  these  integrals  are  of  order  higher  than  0(1).  Hence  the  term  that  will  give  rise  to  a 
possible  generalized  Cauchy  kernel  is  Ri\a(a,x\,  t\).  After  some  very  tedious  asymptotic 
analysis,  we  find, 

R2\a{ot,x\M)  ~  (cja2  +  C2Q  +  C3  +  ^e_(ll+<1+2/tl)a 

+  (e\a  +  e*2  +  ^  +  ^)e(ll+<l)o  (3.44) 

where  expressions  for  the  variables  appearing  in  (3.44)  are  given  in  Appendix  B.  It  should 
be  noted  that  c*'s,  and  e*'s  are  not  constants  but  functions  of  x\  and  t\.  Again  here  only 
the  terms  involving  0(a2),  0(a)  and  0(1)  give  rise  to  singular  terms.  Separating  these 
singular  terms  from  K\\  we  can  write, 


t\)da 


(3.43) 
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#n  (a,  xi ,  *1)  =  #11  (a,  *1,  *i)  +  #11  (a,  *1*  *i) 
where 

Kh(°,*i,*i)  =  (0^)1  {(ci*2 +<$«+<*)< 


,-(xi+ti+2hi)a 


(3.45) 


+  (e;a  +  c5)c^-Wl>t|  da 


(3.46) 


#il(a»*i»*i)  = 


4(1 -\))  {Jo  {■Re(^(a>^))cos(a(ti  “  *0) 

-  (jm(ki(a,y))+4:(^— j))sm(o;(ti  -zi))}da 


{  (Re(ki(a,y)  -  i^e(fci00A(o:,y))jcos(a(i1  -  *i)) 
-  (/m(fei(a,y))+4^^-|)  - 

Imfiioohia,  y)))sin(a(ti  -  xi  ))}<*« 

,00  ^ 

|#e(fci00/l(a,2/))cos(a(t1  -  *i)) 

-  Im  (£ioo*(a,y))sm(a(ti  -  *i))}da 


+  2  jf  ^i?2ia(a,xi,  ti)  -  (c|a2  +  cga  +  eg)  e 
—  ^e}a  +  e^1+<1^  da  | . 


-(xi+ti+2/ii)a 


(3.47) 


Upon  integration  #fi(a,  x\,  t\)  can  be  simplified  as 


■  (iSi;) 


{ 


2  cl 


+ 


C2 


(^l  + 1\  +  2hi)3  (xi+fi+2hi)2  (x\-\-t\Jr2h\) 


+  (xi  +  ^i)2  +  (*i  +  *i)  y  ^3'48^ 

This  is  the  generalized  Cauchy  kernel  in  the  sense  that,  even  though  bounded  for  x\  ^  0, 

fi  7^  Oandxi  ^  —  ^i,  5^  -  hi,  it  becomes  unbounded  for  xi  — >  0,ti  — ►  Oor 

xi  — *  —  hi,t\  — >  —  hi  simultaneously.  If  we  let  /?  go  to  0  in  c\,  c%,  C3  we  recover  the 
following  familiar  form  of  the  generalized  Cauchy  kernel  for  homogeneous  materials, 


+ 


6(xi  +  h\) 


4(xi  +  hx)5 


(xi+ii+2hi)  (xi+fi  +  2hi)2  (xi+fi  +  2hi)3 


(3.49) 


iCi2(a,xi,t2)  : 

The  expression  for  if  12  (a,  xi,  f2)  is  given  in  Appendix  A.  For  easy  reference  it  has 
been  reproduced  here 


#12(0,0:1)  =  2 jf  R2ib(a,x1,t2)da  .  (3.50) 

After  performing  some  tedious  asymptotic  analysis  on  R2\b(,a,x\,  i2)  we  find  thatfor 


a 


00, 


R2ib(cx,x1,t2)  ~  iQ  +  7*2  +  ^  +  ^)e(xi~h)a,  (3.51) 

where  expressions  for  the  variables  are  given  in  Appendix  B.  Again  here  only  the  terms 
involving  0(a)  and  0(1)  give  rise  to  singularities  in  the  kernel.  Separating  these  singular 
terms  from  K\2  we  can  write 
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(3.52) 


(3.53) 


(3.54) 


(3.55) 


-K~2i(o;,3;2,*i)  : 

The  expression  for  K2\  ( a. ,  x2,t\)  is  given  in  Appendix  A.  For  easy  reference  it  has 
been  reproduced  here 


poo 

K2\(pL,X2,ti)  =  2/  R22a{a,X2,ti)da 

Jo 


(3.56) 


After  performing  some  tedious  asymptotic  analysis  on  Ri2aipL,x2-,  t\)  we  find  that, 


•R22fl(o:>2;2)^l)  ~  +  ^2  +  ^  +  ^2  +  ^)e  ^  (3.57) 

where  expressions  for  the  variables  are  given  in  Appendix  B.  In  this  case,  too,  only  the 
terms  involving  0(a)  and  0(1)  terms  give  rise  to  singularities  in  the  kernel.  Separating 
these  singular  terms  from  K2\  we  can  write, 
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(3.58) 


K<2i(a,X2,ti)  =  K2i(oc,X2,t\)  +  K2i(a,X2,ti) 
where 

ma,X2,h)  =  (^)/°°{  + 


Kh(a,x2,ti)  =  (^R22a(a,X2,ti) 

-  (tla  +  t^je-^-^da 
Upon  integration  if|i  (a,  X2,t{)  can  be  simplified  as 


K2i(a,X2,ti) 


-{*?) 


(3.59) 


(3.60) 


J_s 

l  (*2  - 


+ 


*l)2  (*2  -  U) 


} 


(3.61) 


The  expression  for  if 22  (<*,  ^2^2)  is  given  in  Appendix  A.  For  easy  reference  it  has 
been  reproduced  here 


K22(a,X2,t2)  —  2 R22b(a,X2,t\)da  .  (3.62) 

After  performing  some  tedious  asymptotic  analysis  on  R2%(a,X2,  £2)  we  find  that 
R22b(a,  X2,t2)~  (ofa  +  °2  +  §  +  %  +  e‘(x2+t2)Q 


+  (p|a2  +  e(x2+<2-2/i2)o: 


(3.63) 
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where  expressions  for  the  variables  are  given  in  Appendix  B.  Again  here  only  the  0(a) 
and  0(1)  terms  give  rise  to  singular  terms.  Separating  these  singular  terms  from  K22  we 
can  write, 


K^(a,X2,h)  =  K$2(a,x2,t2)  +  K%2(a,x2,t2) 
where 

K^(a,x2,t2)  =  + 


(3.64) 


+  ( pW  +P2<*  +p£)e(-x*+ti-2h*^da, 


(3.65) 


K^(a,x2,t2)  =  (R2Zb(a,x2,t2)  -  (oia  +  o^ 


,-(x2+h)ct 


+  (pi^2  +P$a+  e(xt+h-2h)c^  da  | 


(3.66) 


Upon  integration  K^(a,  x2,t2)  can  be  simplified  as 


K22(a,x2,t2) 


■m 

i 


°i _ + _ °L 


(t2  +  X2 )2  (^2  +  X2) 


+ 


2  Pi 


+ 


P2 


+ 


Ps 


(2h2  —  t2—  X2)3  (2/i2  —  1 2  —  X2)2  (2h2  —  t2  —  x2) 


(3.67) 
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If  we  now  substitute  expressions  for  p\,  pf,,  p%  we  obtain  the  following  generalized 
Cauchy  kernel  for  a  homogeneous  medium  with  a  stress-free  boundary  at  rr2  =  /*2- 

4(7*2— £2  )2 _ 6(7*2— £2) _ 1 

(2^*2  ~h~  X2)3  (2h2  —  t2  —  a?2)2  (2/12  -  £2  -  ^2) 

The  kernel  in  (3.68)  is  identical  to  that  found  by  Rizk[67]  and  Wu  andErdogan[46]. 

3.3  Evaluation  of  singularities  and  expressions  for  stress  intensity  factors 

With  the  knowledge  of  the  possible  singular  kernels  obtained  in  the  previous 
section,  we  can  now  evaluate  the  singularities  of  the  unknown  functions  g\  and  §2  at  the 
crack  tips  for  different  crack  configurations.  Each  crack  configuration  has  been  considered 
separately  as  the  singularities  may  be  different  in  each  case. 


3.3.1  Internal  crack  in  the  Nonhomogeneous  material  only 

For  this  case,  the  relevant  integral  equation  is  (3.36)  without  the  coupling  term 
and  can  be  written  as 


I  n  gigo 

*J-bl  (*l-*l)  1 


^J-bi 


=  piM 


K  -f-  1 
4/ioe^i  ‘ 


(3.69) 


From  (3.48)  it  can  be  seen  that  K\\  (rri,  t\)  may  have  singular  terms  only  when  x\,  t\ 
either  both  go  to  zero  or  both  go  to  —  h\.  Hence  in  the  case  of  an  internal  crack 
Ku  (*i,ti)  is  a  bounded  term  and  does  not  contribute  to  the  singularity.  The  only 
singular  term  is  the  dominant  term  — - —r .  The  integral  equation  can  now  be  written  as 


irja&L*. 

KJ-b,  (*i-*i) 


^l(xi) 


(3.70) 
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where  ^i(rri)  represents  all  the  bounded  terms.  To  study  the  behavior  of  g\  (fi)  near  the 
endpoints  —  ai  and  —  b\,  let  us  assume  that 


Sliti)  — 


gf(*i) 

(fi  +  i>i)Ql  (  —  ai  —  t\)& 


(3.71) 


where  g\ (ti)  satisfies  the  Holder  condition  on  —  &i  <  ti  <  -  a\.  Also,  the  singularities 
ai  and  /5i should  satisfy  the  condition  0  <  Re{a\,  (3\)  <  1.  Defining  a  sectionally 
holomorphic  function 


i  rai 

m  =-/  ■ 

7T  J-h 


9i(h) 


dt\ 


(*i  ~  z) 

and  substituting  (3.71)  into  (3.72)  we  obtain, 


(3.72) 


F(z\  =  I  P 

nj-bi  (*i  +  h)ai  (t l  +  ai)A(*i  -  «) 

Following  Muskhelishvili[55],  the  asymptotic  behavior  of  (3.73)  near  the  end  points  -  a\ 
and  —  £>i  can  be  expressed  as 


F{z)  = 


9t(  -  bje™' 


^i(-Qi) 


( b\  —  ai)&  (z  +  6i)Q>sin(7ro;i)  (6j  —  ai)Ql  (z  +  ai)#sin(7r/3i) 


+  F0(z) 


(3.74) 


where  Fq(z)  has  singularity  less  than  Now  using  the  Plemelj  formula,  we  obtain, 

o).  (3.75) 

From  (3.74)  and  (3.75)  it  then  follows  that 
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I  f  a>  9i(ti)  dt  _  Pi(  ~  6i)cot(7rai) _ g\(  -  Qi)cot(7r/?i) 

77 J-bi  (M-Xi)  1  (6i  -  oi)A  (Xl  +  61)^1  (b1-a1)*(-a1-x1)& 


+  -Fo(*i). 

Substituting  (3.76)  into  (3.70)  we  obtain, 


(3.76) 


ffi(  -  jljcotQrai) _ ffx(  -  a1)cot(7T/91) 

(h  -  ai)A  (xi  +  &i))Q>  (£>i  -  oi)°i  (  -  ax  -  xx)A 


(3.77) 


Multiplying  the  entire  equation  by  (sx  +  bx)a'-  and  taking  limit  as  zx-4  -  6X,  we  obtain, 


pi(-&i)cot(7rai) 

(fei  -  oi)A 

(3.78) 

COt(7TQ!l)  =  0. 

(3.79) 

Now  multiplying  the  equation  (3.77)  by  ( -  ax  -  *i)A  and  taking  limit  as  zx->  -  ah  we 
obtain, 


gj(  ~  Qi)cot(7r/?i) 

Q>i  -  ai)Ql 
or 

cot(7r/5x)  =  0. 

The  only  acceptable  roots  of  (3.79)  and  (3.81)  are  ax  =  /3X  =  Thus  we  recover  the 
square  root  singularity  at  the  crack  tips. 


(3.80) 

(3.81) 


Stress  intensity  factor: 


The  stress  intensity  factor  at  irregular  points  -  6xand  -  ax  can  be  defined  as, 
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(3.82) 


h(~h)  =Bm  \/2(  -  i»i  -  xi)  <jyy i(*i,  0), 

Xi  — >~Ol 


fci(-oi)=lim  a/2(xi  +  01)  rn^xi,  0) 


(3.83) 


where  x\  lies  outside  the  crack.  It  should  be  noted  here  that  the  left  hand  side  of  the 
singular  integral  equation  in  (3.36)  is  the  expression  for  <jyy\  (xi,  0)  inside  as  well  as 
outside  of  the  crack.  Using  that  expression  we  can  write  (3.82)  as, 

- h)  =  tti  A,  e*‘  ^  - 61  - li){;  h ' 

+  ~Jb  ^  g2{h)Kn{x\,t2)dt<^  (3.84) 

Since  K\\ (xi,  ti)  and  K-& (xi,  t2)  are  bounded,  we  can  rewrite  (3.84)  as 


h(-h)  = 


— T  Mm  e@Xl  \/2(-i>i-xi) 

K  “t*  1  xl  —>”0i 


I—  [  _)_  bounded  terms} 

U  7T  y_bl  (ti-Xi)  J 


(3.85) 


Defining  a  sectionally  holomorphic  function, 


F(.)  .irji&u 

*V-6, 

and  substituting  the  expression  for  g\(t\)  from  (3.71)  and  noting  that  a\  =  /?i  = 
obtain, 


(3.86) 


5.  we 


£i  ~  &i)e » 

(*>i  -  ai)  2  (2  +  61)2 


(i>i  ~  01)5(2  +  01)2 


r  +  Fq{z)- 


(3.87) 


Using  Plemelj  formulas,  we  can  write 


1 P  9l(t,)  m  g?( - feo _ rf(-° i)  rfl 

^ J-bi  (U  —  ®i)  (61  —  ai)5(  —  i>i  —  £1)2  (61  —  01)2  {x\  +  ai)s 


(3.88) 


In  (3.85)  x\  <  -  61  and  £1  >  -  ai.  Substituting  back  into  (3.82),  we  obtain, 
hi  ~  h)  =  —^7  lim  e^1  y/2(  -  61  -  xi) 

Kj  ~T  I  ®1— 0\ 

r  g?(-&i)  t  gj(-oi) 

''(h  ~  0.1  )i(  -  61  -  xi)*  (6!  -  <21)2  (xj  +  oi)i 

+  6oundedferms|  (3.89) 

which  gives  us, 

(-) 

By  similar  analysis  we  obtain, 

*i(~0i)  =  (3-91) 

3.3.2  Interna]  crack  in  the  homogeneous  material  only 

For  this  case,  the  relevant  integral  equation  is  (3.42)  without  the  coupling  term 
and  can  be  written  as 

lL^i)dt2+li92(t2)K22{X2't2)dt2  =K(i2)^7'  (3-92) 


The  analysis  in  this  case  is  similar  to  that  in  Section  3.3.1  and  we  find  singularities  again 
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3.3.3  Internal  cracks  in  homogeneous  and  nonhomogeneous  materials 

In  this  case  we  obtain  two  coupled  singular  integral  equations  as  follows: 

^  J —b\  (^1  xl)  K  J —b\ 

+  tt/  92(t2)Kn(x1,t2)dt2  =  »  (3-95) 

~[  v/2 ~  ^  H — [  9i(t\)K2i(x2,ti)dti 

nJoi  \h~X2)  tt 7_bl 

1  yk  «  +  1 

+  -/  92(h)  K22(x2,t2)dt2  =  P2(x2)— - .  (3.96) 

*  Jot  4/i0 

From  (3.48),  (3.55),  (3.61)  and  (3.67),  we  see  that  K\\,  K\2,  K2\,  K22  have  no  singular 

terms  as  long  as  the  cracks  are  internal.  Hence  the  only  contribution  to  the  singularity 

1  1 

comes  from  the  dominant  terms  — - r  and  — - - .  The  analysis  is  then  the  same  as 

C t\  ~  *1)  <fo  ~  *2) 
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the  previous  cases  considered  in  Sections  3.3.1  and  3.3.2  and  we  recover  the  expected 
square  root  singularities  i.e.  a\  =  =  /%  =  \ . 

Stress  intensity  factor: 

Since  the  coupling  kernels  do  not  have  any  singular  terms  as  long  as  the  cracks  are 
internal,  this  case  is  similar  to  cases  considered  previously.  Hence  the  expressions  for 
stress  intensity  factors  would  be  identical  to  those  found  in  Sections  3.3.1  and  3.3.2. 


3.3.4  Edge  crack  in  the  nonhomogeneous  material 

For  the  case  of  an  edge  crack  in  the  nonhomogeneous  medium,  the  singular 
integral  equation  be  written  as, 

x/kl  ~0kdt'  +  IL,  S^)Kn(xuh)dtl  =Pl(Xl)^±L.  (3.97) 


In  this  case  apart  from  the  dominant  term  — - -  ,  we  have  contributions  to  the  singular 

(fi  -  *i) 

behavior  of  the  solution  from  the  Fredholm  kernel.  From  (3.48),  if  we  carry  out  a  careful 
analysis,  we  find  that  the  singular  part  of  K\\  (a?i,  f  i)  can  be  written  as, 


_ _l _  6(xi  +  hi) 


4(xx  +  hi)2 


(xi+fi+2hi)  (xi+fi+2hi)2  (xi+fi  +  2hi)3 


(3.98) 


Putting  all  the  bounded  parts  to  the  right  hand  side,  we  can  write  (3.97)  as, 


\  f_h  (('ti  Z  Xl)  +Kh(xi’tij)gi(ti)dt1  =®i(*i). 

Again  assume  that, 

9i(h) 


5f(*i) 


(ii-(-/M))0l((-ai)-t,)A 


(3.99) 


(3.100) 
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Defining  a  sectionally  holomorphic  function 


1  f^1 

F (*)  =  -  / 
nJ-h, 


9i(h) 


dt\ 


-hi  (tl-z) 
and  substituting  (3.100)  into  (3.101)  we  obtain 


(3.101) 


>  _  I  r« _ M 

2  IT  J_ht  (ii  —  (  —  ht))*‘  (t,  —  (  —  ai))«(ti  -  2)  ’■ 

Following  Muskhelishvili[55],  the  asymptotic  behavior  of  F(z)  can  be  written  as 


(3.102) 


F(z) 


9\{-hx)e ^ 


9l(  ~  ai) 


(hi  -  a\)&  (z  —  (  -  ^i))Qlsin(7rai)  (hi  -  ai)ai  (z  -  (  —  ai))^sin(7r/?i) 


+  F0(«)  (3.103) 

where  Fq(z)  has  singularity  less  than  ai,  /3\.  Again  following  the  same  procedure  as 
followed  in  Section  3.3.1,  we  can  write, 

1  p1  0i(*i)  dt  =  Pi(~  fei)cot(7rai) _ gi(  -  oi)cot(7r/3i) 

nj-hi  (*i-*i)  1  (/ii  -  ai)A  (xi  + /ii)Qi  (hi  -  oi)ai  (  -  ^  -  xi)A 

+  Fo(xi).  (3.104) 

1  p  gi(*i)  ,,  =  jjC  ~  M  ,  \ 

TtJ-hi  (ti+xi  +2hi)  1  (/»i  -  oi)A  (xi  +  /ii)Q>sin(7ro:i)  0  1  ‘ 

(3.105) 


Differentiating  (3.102)with  respect  to  xi  we  obtain, 
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I  n _ 9i(h) _ d  =  aifij-ht)  _d_ 

J —hi  (ti  +xi  +  2h-[)2  (/ii  —  ai)^  (xi  + /ii)Qi+1sm(7TQi)  dx\  1 


(3.106) 


Differentiating  once  again,  we  obtain 


i  r  dt  =  Qi(Qi+i)gi(-fei)  ,  Tp*f  s 

nj-hi  (ti+x\+2hi)3  1  2(hi  —  ai)A  (a^  +  /ii)ai+2sin(7ra1)  dx\  0  2:1  ’ 

(3.107) 


Using  (3.104  —  3.107)  and  (3.98),  (3.99)  may  be  expressed  as, 


gj(  ~  hi)cot(7rai) _ gj(  -  ai)cot(7r/?i) _ gj(  -  hj) _ 

(h\  —  ai)A  (xi  +  hi)Q >  (Ai  —  ai)°i  (  —  a\  —  x\)&  (h\  —  aj)A  (xi  +  hi)ai  sin(7ra:i) 


+ 


_ 6aiffi(-/ti) _ 

(h\  —  a{)fo  (x\  +  /ii)aisin(7rai) 


2ai(ai  +  l)gi(-  hi) 

(hi  —  ai)#  (xi  +  fti)Qlsin(7rai) 


(3.108) 


Multiplying  (3.105)  by  ( —  ai  —  xi)&  and  letting  x\  go  to  —  oi,  we  find  the 
characteristic  equation, 

cot(7r/?i)  =0,  =>  Pi  =  §.  (3.109) 

On  the  other  hand,  multiplying  (3.105)  by  (xi  +  hi)ai  and  letting  xj  go  to  -  hi  we 
obtain  the  characteristic  equation, 

cot(7rai)  -  2a1(a1  +  1)  +  6<*i  -  1  =0.  (3.110) 

ai  =  1  is  an  unacceptable  solution.  Only  acceptable  solution  is  a\  =  0. 
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Stress  intensity  factor: 


Observing  that  there  is  no  singularity  at  x\  =  —  hi,  the  only  point  of  interest  is 
x\  =  —  a\.  We  can  define  the  stress  intensity  factor  at  xi  =  —  a\  as 

&l(  — ai)  =  lim  y/2(xi  +ai)<7yyi(xi,0).  (3.111) 

X\ — ► — 

Following  a  procedure  similar  to  that  in  Section  3.3.1,  we  obtain, 

fci(  -  01)  =  -  -  ai).  (3.112) 

/C  i  A 

3.3.5  Crack  in  the  nonhomogeneous  material  terminating  at  the  interface 

The  relevant  integral  equation  for  this  case  can  be  written  as, 

In  this  case  if  we  carry  out  a  careful  analysis  of  K\\  (xi,  ti),  we  find  that  K\\  (x\,  t\) 

contributes  no  singular  terms  to  the  integral  equation.  The  only  singular  term  is  the 

dominant  term  — — - — r .  Hence  this  case  becomes  quite  similar  to  that  discussed  in 
(ti  -  a?i) 

Section  3.3.1  and  if  we  follow  the  same  procedure  as  in  Section  3.3.1  we  find, 

=  A  =  j-  (3-114) 


Stress  intensity  factor: 

The  stress  intensity  factors  are  identical  to  Case  I  and  are  given  as, 


fci(  -  £>i) 


K  -|-  1 


(3.115) 
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(3.116) 


*i(0)  = 


4^o 

K  ■+■  1 


3.3.6  Crack  in  the  homogeneous  material  terminating  at  the  interface 

The  relevant  integral  equation  for  this  problem  can  be  written  as, 


-  [  y92^  dt2 


1  fbl 

+  ~  I  92 (^2 )dt2 K22, (x2 ,  t<2 )  =  £2(^2) 
"  Jo 


K  +  1 

4^o 


(3.117) 


In  this  case  again  we  find  that  ^22(^2,  h)  contributes  no  singular  terms  to  the  integral 

equation.  The  only  singular  term  is  the  dominant  term  - — - — - .  Hence  this  case  becomes 

(*2  “  *2) 

exactly  similar  to  that  discussed  in  Section  3.3.2  and  if  we  follow  the  same  procedure  we 
find, 

c*2  =  /%  =  j.  (3.118) 

This  result  was  also  obtained  by  Erdogan  et.al.[42]. 


Stress  intensity  factor: 

The  stress  intensity  factors  are  identical  to  those  of  Section  3.3.2  and  are  given  as, 


*i(0)  = 


4/^q 

K  +  1 


(3.119) 


fci(&2) 


4^o 

K  +  1 


(3.120) 


3.3.7  Crack  passing  through  interface 

For  this  case,  we  obtain  two  singular  integral  equations  as  follows: 
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~  [  7p\ +  ~  /  9l(t\)Kll(xl>ti)dti 

^  J  — b\  1  7^  J  — b\ 


+ 


1  fa  K  +  1 

-JQ  92(t2)Ku(x1,t2)dt2=p1(xl)-^^, 


(3.121) 


~f  •  2  H —  f  gi(ti)K2i(x2,ti)dti 

irj o  (I2-X2)  tt  J _j,1 


+ 


1  H’2  K  +  1 

-  /  92(t2)K22(x2,t2)dt2  =  P2(x2)-. - • 

7T  7o  4^/0 


(3.122) 


As  was  observed  in  Section  3.3.5  and  Section  3.3.6,  Fredholm  kernels  K\\{x\,  t\)  and 
K22 (x2i  h)  do  not  contain  any  singular  terms.  Hence  they  do  not  contribute  to  the 
singularity.  The  coupling  terms,  however,  do  contain  singular  terms  and  after  some 
analysis  the  singular  terms  are  found  as, 


K^(xi,t2) 

=  , .  ,  -  hi  <  xi  <  0, 0  <  t2  <  h2, 

(h  -  *i) 

(3.123) 

K21  (x2,h) 

=  t,  ,  —  h\  <  t\  <  0, 0  <  x2  <  h2. 

Cm  -  *2) 

(3.124) 

It  should  be  noted  here  that  ti,  x\ are  always  negative  while  t2,x2  are  always  positive. 
Hence  (x\,  t2 )  is  singular  only  when  both  and  t2  go  to  zero.  Similarly,  K2l  (x2,  t{) 
is  singular  only  when  both  x2  and  £1  go  to  zero.  Putting  all  the  bounded  parts  to  the  right 
hand  side,  (3.121)  and  (3.122)  may  be  expressed  as, 


nJ-b^t  i-*i)  nJQ 


92(h) 

0 h  -  ^1) 


dt2  =  Vi(*0f 


(3.125) 
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(3.126) 


1  f  if 

71  J-bl(ti-X2) 


=  ®2(*2)- 


where  ^i(xi)  and  V2(x2)  represent  the  bounded  terms.  Let  us  assume  that 


0i(*i)  = 


9 jfa) 

C*1  +&i)^C-ti)A* 


(3.127) 


92(h)  = 


g*Jh) 

(t2)°2  (62  -  *2)A 


(3.128) 


where  ^J(f  1)  is  Holder  continuous  on  —  £>1  <  t\  <  0  and  g^ih)  is  Holder  continuous  on 
0  <  t2  <  b2.  It  should  be  noted  that  at  zero  there  can  be  only  one  irregular  point.  Hence 
A  =  ct2-  Defining  sectionally  holomorphic  functions 


Fi(z) 


1  f°  gigi) 
nj_bl  (tj  -  2) 


df  1 , 


(3.129) 


*!»(*) 


92(h) 
(h  ~  z) 


dt2 


(3.130) 


and  substituting  (3.127)  and  (3.128)  into  (3.129)  and  (3.130)  we  obtain, 


1  Z*0 

*i(*)  =  £/ 

77  J-h 


9i(h)e 


*101 


(fl  -+-  £>l)Q1(*l)^(*l  “  •2) 


df  1 , 


(3.131) 


F2(z)  =Ub 

77  Jo 


{t2)«hh-h)h{h-z) 


dt2. 


(3.132) 


Following  Muskhelishvili[55],  we  can  write, 
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(3.133) 


Fi(z) 

and 

F2(z) 


fj(  -  h)e™ _ fi(0)  ,  r  ,  , 

(bi)fi(z  +  6i)Glsin(7rai)  (bi)Ql  (2)#  sin(7r/?i)  1 

m _ _ +F2o(z) 

(b2)/32(z)Q^sin(7TQ2)  ( b2)Q2(z  -  i>2)&sin(7T/%) 


(3.134) 


Using  Plemelj  formulas  as  before,  we  obtain, 


1  f°  9i(ti)  _  ff|(  -  bi)cot(Trai)  _  ^(0)cot(7r/?i) 

vj-bit  1-X1)  1  (&i)A(xi +  &i)Ql  (&i)ai(-zi)A  10  1 


1  fh  92(h)  ,  ^(0)cot(7rQ-2)  _  g^(62)cot(7T/%) 

nJo  (<2-«2)  2  (b2)fo(x2)a2  (62)Q2(62  -  x2)^  20  2 


(3.135) 


(3.136) 


Also, 

I  /'°  gl(«l)  ..  _  _  g*(Q)  ,  ™ ,  > 

i~*2)  1  (5i)Ql(x2)Asin(7r/5i)  10  2 

1  92(h)  .. _ gl(0) _  ,  ,  x 

nJo  ($2-*i)  2  (52)^(  —  xi )Q2 sin (7tq!2)  20  1 


(3.137) 


(3.138) 


Substituting  back  into  (3.125)  and  (3.126)  we  obtain, 


gi(  ~  &i)cot(7rai)  _  g|(0)cot(7r/3i)  _ fl|(0) _  *r  . 

(bi)fr(xi  +  bi)a'  (bi)ai(  -  xi)&  (52)&(-x i)°®sin(7rQ:2)  ^  1 

(3.139) 


g^(0)cot(7rQ2)  _  ff2(52)cot(7T/%) _ gj(0) 

(b2)^(x2)°: 2  (62)°:2(62  —  x2)&  (6i)«i(x2)Asin(7r/?i) 
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^2(x2). 

(3.140) 


Multiplying  (3.139)  by  (x2  +  6j)Q;  and  letting  x\  go  to  —  b\  and  multiplying  (3.140)  by 
(62  -  £2)^  and  letting  x2  go  to  62,  we  obtain  the  characteristic  equations, 


cot(7rai)  =  0 

and 

^  ai  —  2 

(3.141) 

cot(7 r/%)  =  0 

=►  /%  =  §• 

(3.142) 

Now  multiplying  (3.140)  by  (  —  x\)&  and  letting  x\  go  to  zero  and  multiplying  (3.141)  by 
X2  and  letting  x2  go  to  zero,  keeping  in  mind  that  (3\  =  02, we  obtain, 


gj(0)cot(7rft)  ^(0) 

(i>i)Ql  (62)^sin(7r/51) 


(3.143) 


ffl(O)  ^(0)cOt(7T/?i) 

(&i)Qlsin(7r/?i)  (62)^ 


(3.144) 


From  equation  (3.143),  we  obtain, 

02(0)  _  Pi(°)c°s(7r/g1) 

(&2)/%sin(7r^i)  (6!)°!  sin(7r/5i)  * 

Substituting  (3.145)  into  (3.144)  we  obtain  the  following  equation 

_  (0)  .  ^(0)c°s2(7r^i)  _ 

(f>i)aisin(7r/5i)  (&i)aisin(7r/?i) 

Solving  (3.146)  we  obtain, 

0i(O)sin(7r/?i) 

0 b i)Ql 

For  non-trivial  solution,  we  obtain, 

A  =  (0,  ±1,  ±2 .... ). 


(3.145) 


(3.146) 


(3.147) 


(3.148) 


The  only  acceptable  solution  is 


We  can  now  write  the  unknown  functions  as 


9l(h) 


am 

(ti  +  h)i 


(3.150) 


92(h)  = 


92(h) 
(i>2  -  h)* 


(3.151) 


From  (3.143)  and  (3.144)  it  is  clear  that  gl(0)  and  g%(0)  are  related  and  the  relation  can 
be  expressed  as 


^(0)cos(7t/3i)  |  flg(0)  _0 

(6i)5sin(7r/3i)  (b2)2sin(7r/3i) 


(3.152) 


_  ffi(0)cos(7T/31)  &(0)  =  Q 

sin(7r/?i)  sin(7r/5i) 

As  /?i  =  0,  (3.153)  will  be  satisfied  if 
92(0)  =  gi  (0). 


(3.153) 


(3.154) 


The  result  (3\  —  0  implies  that  there  is  no  power  singularity  at  the  interface.  There  is  a 
possibility  of  a  weaker  logarithmic  singularity.  We  investigate  it  as  follows.  Defining 
sectionally  holomorphic  functions 


(3.155) 


iV)  =lf 

ft  Jo 


92(h) 

(h-z){ 


(3.156) 
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and  substituting  (3.150)  and  (3.151)  into  (3.155)  and  (3.156)  we  obtain. 


*(*) 


l  f  _ 


+  &l)5(^i  ~  z) 


dt\ , 


1  f 62 
^  Jo 


g*2(h) 


■ dt 2  . 


(62  -  <2)2^2  -  *) 

Investigating  in  the  neighborhood  of  zero,  we  obtain, 


Fi(z)  ~^^log(z) 


(3.157) 


(3.158) 


(3.159) 


F2(z)  ~-^^p-log(z) 


H 


Using  Plemelj  formulas,  in  the  neighborhood  of  zero,  we  can  write 


if 

^  ''—bi 


V 

nj 0 


(<1  -  *1) 

h  92(h) 
(to.  -  XI ) 


dt2 


1  t 

irJo 


92  (t2) 
(to  ~  Xo) 


dt% 


(3.160) 


^  lQg|*i| 
b\  * 

(3.161) 

92(0)1 

1  log  *1 

1.2  7T 
°2 

(3.162) 

b\  * 

(3.163) 

h ?  K 

(3.164) 

It  should  be  noted  that  in  (3.159  —  3.164),  there  are  bounded  terms  that  have  not  been 
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mentioned.  Hence  the  ~  sign  has  been  used.  Substituting  (3.161  —  3.164)  into  (3.125) 
and  (3.126),  we  obtain, 


^  -  log  In  I  -  ^  -  log  Is,  I  =»;(*!),  (3.165) 

b\  w  6|  * 


&21  -  log  M  -  ^  -  log  |*2|  =  1-5(12)  (3.166) 

L2  7T  j,2  7T 

°1  °2 

where  W|(xi),  ^(a^)  are  bounded  functions.  Using  (3.150),  (3.151),  (3.154)  it  can  now 
be  clearly  seen  that  the  coefficient  of  the  logarithmic  term  cancels  out.  Hence  there  is  no 
singularity  at  the  interface.  This  is  an  important  conclusion.  This  should  be  compared  with 
the  case  of  bonded  layers  of  dissimilar  materials.  It  was  shown  by  Ritz[67]  and  Lu, 
Erdogan[26, 27]  that  when  a  crack  crosses  the  interface  between  two  dissimilar  materials 
there  is  an  irregular  point  at  the  interface  apart  from  the  two  irregular  points  at  the  crack 
tips.  The  strength  of  the  singularity  at  the  irregular  point  at  the  interface  has  to  be 
determined  from  the  asymptotic  analysis  at  the  interface.  In  our  study  of  FGMs,  it  was 
found  that  there  is  no  irregular  point  at  the  interface  or,  in  other  words,  there  is  no 
singularity  at  the  interface. 

Since  there  is  no  singularity  at  the  interface,  we  can  pose  (3.121)  and  (3.122)  as 
one  singular  integral  equation  as 


1  /*■  g(t) 

7T y_fol  (t  -x) 


1  rh 

dt  4 —  /  g(t)dtK(x,t) 

*  J-h 


=  p(x) 


K  +  1 

4/20 


where 


t  <  0 
t>  0 


(3.167) 


(3.168) 
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1 


(x,t)  <  0 
x  <  0,f  >  0 


K(x ,  t) 


r#n (x,t) 
#12  OM) 

#21  OM) 
,  #22 (x,t) 


(t-x) 

_1 _ 

(t-x) 


x  >  0,t  <  0 

( x,t )  >  0 


(3.169) 


Pi  (?) 


x  <  0 


p(x)  =  ^  ePx\ 

P2^)  x  >  0 


It  should  be  noted  that  the  only  irregular  points  are  x  — 
write, 

9(t)  = 


(3.170) 


b\  and  x  =  b2.  Hence,  we  can 


9*(t) 


(t  +  61)5(62  ~t)t 
where  g*{t)  is  a  bounded  function. 


(3.171) 


Stress  intensity  factors: 

The  analysis  is  similar  to  previous  cases.  Keeping  in  mind  that  irregular  point 
x  =  —  b\  lies  in  the  nonhomogeneous  material  while  x  =  b2  lies  in  homogeneous  we  can 
write  the  stress  intensity  factors  as, 


ki(-h) 


K  +  1 


b\  +  b2 


9*(  ~  61), 


(3.172) 


h(b2)  = 


4po 
K  +  l 


bl  +  6; 


-9*(b2). 


(3.173) 


3.3.8  Edge  crack  passing  through  interface 

This  case  is  covered  by  some  of  the  previous  cases.  The  analysis  for  the  singularity 
at  x  =  —  h\ is  identical  to  that  for  an  edge  crack.  The  singularity  at  the  interface  can  be 
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evaluated  in  a  manner  identical  to  that  in  Section  3.3.7.  If  we  carry  out  the  analysis,  we 
find  that, 

=  A  =  £*2  =  0,/%  =  (3.174) 

Here  again  we  can  represent  the  problem  by  a  single  singular  integral  equation  as  follows: 


-  P  7pr-,dt  +  i  P  g(t)dtK(x,t) 

J  —hi  \y  % )  IT  J  _ h\ 


=  p(x) 


K  +  1 

4^0 


where 


t  <  0 
t  >0 


f  0 M) 


K(x,t)  =  { 


#12  OM) 


#21  OM) 


(  Kzz(x,t) 


1 

(t-x) 

1 

(t-x) 


x,t  <  0 
x  <0,t  >  0 

x  >  0,t  <  0 
x,t  >  0 


{pi(x)  x  <  0 
P2(x)  x  >  0. 


g(t)  can  be  represented  as 


(3.175) 


(3.176) 


(3.177) 


(3.178) 


9(t) 


g*(*) 

0>2  ~  t)'* 


Stress  intensity  factor: 


The  stress  intensity  factor  at  x  =  62  can  now  be  written  as 


(3.179) 


h(h)  =  - 


(3.180) 
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CHAPTER  4 


NUMERICAL  PROCEDURE 

The  singular  integral  equations  derived  in  Chapter  3  can  not  be  solved  in  closed 
form.  Hence  there  is  a  need  to  use  some  kind  of  numerical  procedure  to  solve  these 
equations.  In  the  literature  various  approaches  have  been  studied.  The  two  most  common 
approaches  are  the  quadrature  method  and  the  series  expansion  method.  In  this  study,  the 
expansion  technique  described  in  [57,58]  is  used.  In  this  chapter  the  method  will  be 
explained  in  some  detail.  For  different  crack  configurations,  the  numerical  procedure  has 
to  be  modified.  These  modifications  will  also  be  explained  in  this  chapter. 

4.1  Numerical  procedure  for  a  general  case 

For  the  most  general  crack  configuration  the  singular  integral  equations  are  given 
by  (3.36)  and  (3.42)  and  may  be  expressed  as 

l  fh  K  +  1 

+  -J  Kn(xi,t2)g2(t2)dt2  =  Pi{x{)---~^ ,  -6i<*i<-oi,  (4.1) 

1  />_ai 

-  K2i(x2,ti)gi(ti)dt1 
nj-h 

02  <  x2  <  b2.  (4.2) 
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In  equation  (4.1),  —  £>i  <  x\  <  —  a\  and  in  equation  (4.2),  <  62.  Now 

normalizing  the  integral  equations  with  respect  to  the  crack  lengths  as  follows, 


h  = 


x-i  = 


h  = 


x2  = 


we  obtain 


fei  —  b\  +  g-i 

2  ri  2  • 

(4.3) 

£>1  -  dl  +  0-1 

2  31  2  ’ 

(4.4) 

62  —  02  ,62  +  02 

2  r2+  2  • 

(4.5) 

^2  ""  a2  .  &2  +  G2 

2  s2+  2  ■ 

(4.6) 

1  yi  r  i  ^  -l  i  /*i  ^ 

-  /  7 - r  +-^n(si,ri)  g1(rl)drl  +-  Kv(si, r2)?2(r2)^2 

ttJ-i  L(^i  -  si)  J  7TJ-1 


=  ^1(Sl)l^7,lSl  <  2I 


“  J ^2i(s2, ri)g1(r1)drl  +  \  j ^ y~Z +  ^22 0*2, ^2)  ?2(r2>^2 


=  P2(S2)^^,|32  <  !|  (4.8) 


where 

Kij(si,rj)  =  (— —  ■)  Kij(x{,tj),  (4.9) 

?j(rj)  =9j(tj),  (4.10) 
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~  /  N  _  PlM 
pi(Sl)  —  ePx 1  ’ 

P2O2)  =P2M- 


(4.11) 

(4.12) 


From  Chapter  3  it  can  be  seen  that  the  density  function  gj(rj)  can  generally  be  represented 
as  a  product  of  a  bounded  function  and  a  weight  function.  Thus  we  can  represent  gArj) 


?i(n)  = 


am 


(4.13) 


~  >, _ jMn) 

511  i;  (l  +  ri)ai(l-ri)A 

where  Fi(ri)  is  a  bounded  function  and 


(4.14) 


fefa)  = 


gate) 

(t2  -  - 1?)0* 


(4.15) 


(4.16) 


Ur2)  =  (4'16) 
where  F2(r2)  is  a  bounded  function.  The  values  of  the  exponents  ai,  (3\,  a2,  /%  depend  on 


the  crack  configuration  as  shown  in  Chapter  3.  As  Fi(ri)  and  F2(r2)  are  bounded 
functions  we  can  expand  them  in  terms  of  orthogonal  polynomials  as 


-Fi(ri)  =  £AnOn(ri)and 

n=0 


(4.17) 


1*2  (Vo )  =  ^T,BnOn(r<i). 

n= 0 


(4.18) 


The  choice  of  the  orthogonal  polynomial  On(r )  which  is  generally  a  special  case  of  the 
Jacobi  polynomials  Pna'®  (r)  depends  on  the  weight  function.  It  can  be  seen  that  on 
substituting  (4.17)  and  (4.18)  into  (4.11)  and  (4.12)  we  have  ( N  +  1)  unknowns 
Aq...An  and  (M  + 1)  unknowns  Bq..Bm •  Hence  we  need  to  have  (N  +  M  +  2) 
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equations.  This  is  usually  done  by  choosing  collocation  points  where  the  equations  are 
satisfied  exactly.  Though  there  is  no  restriction  on  the  choice  of  the  collocation  points, 
usually  the  roots  of  Chebyshev  polynomials  are  chosen  as  collocation  points.  These  seem 
to  give  better  convergence  in  numerical  results.  The  number  of  collocation  points  depends 
on  the  number  of  unknowns.  Evaluating  the  equations  at  the  collocation  points  we  obtain 
a  system  of  algebraic  equations  which  can  be  solved  by  matrix  inversion.  It  should  be 
noted  here  that  particular  crack  configurations  may  dictate  some  physical  conditions  for  a 
unique  solution.  These  conditions  known  as  single  valuedness  conditions  need  to  be 
incorporated  into  the  system  of  linear  equations.  We  will  now  analyze  the  different  crack 
configurations  separately. 

4.2  Numerical  procedure  for  specific  crack  configurations 
4.2.1  Internal  crack  in  the  nonhomogeneous  material  only 

From  (3.69)  we  can  see  that  the  singular  integral  equation  for  this  case  is 

lLt  +IL,  ffi(<.)jru(*i.«i)<tti  =  pi(*i)^r. 

-  b\  <  x\  <  -  a\.  (4.19) 


Normalizing  using  (4.3)  and  (4.4)  the  singular  integral  equation  can  be  written  as 

nf_i  [(rv^aij  +^n(si’ri)]^i(ri)dri  lSl<1l  (4-20) 


where 

£>  /_  _  \  b\-a\  b\+ai  b\-a\  bi+ai 

Au(si,ri)  =( — - — )Kii( — - — si - - — , — - — ri - - — ),  (4.21) 
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—  /  \  fh  -  ai  h+ai 

9l(r\ )  =  9i( — j n - 2 — 


(4.22) 


?l(«i)  = 


PiM 

ef3x\ 


(4.23) 


Observing  from  equations  (3.79)  and  (3.81)  that  oki  =  A  =  § » we  can  write  (4.14)  as 


?i(n)  = 


Fi(r  i) 


(4.24) 


where 


Fi(ri) 


rfcl  —  a\  b\  +  ai 


)fe) 


(4.25) 


As  Pi  (ti)  is  a  bounded  function,  Fi(ri)  is  a  bounded  function.  Hence  it  can  be  expanded 
in  series  of  polynomials.  The  orthogonal  polynomials  corresponding  to  the  weight  function 

1  m  .  .  4  »  1  (.  .t  /■  .  ,  •  1  t.  T  1  /  \  _  _  ,1 _ 


are  Chebyshev  polynomials  of  the  first  kind.  We  now  expand  Fi(rj)  as  follows , 


*i(n)  =  Y^KTni r,). 


(4.26) 


Equation  (4.20)  can  now  be  written  as 

1  pi  r  1  i  AnTn ^  k  +  1 

dn =?iW  w 

Since  the  crack  is  closed  the  single  valuedness  condition  dictates  that 


(4.27) 


/  9i(ti)dti  =  0 

J-h 

which  can  also  be  written  as 


(4.28) 


,4  AnTn(r \) 


■dri  =  0. 


(4.29) 


Chebyshev  polynomials  satisfy  the  following  orthogonality  relations: 
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(4.30) 


-^n^O^roQ*  l) 


dry  = 


0 

7T 

7T 

2 


n  ^  m 
n  =  m  =  0 

n  —  m  >  0 


Using  these  orthogonality  relations  we  clearly  observe  from  (4.29)  that  Aq  =  0.  The 

singular  integral  equation  in  (4.27)  can  thus  be  written  as 

N 


\S-  i[(dbo+*,,(Sl'n! 


,  "^2,  AnTn(r \)  ,  - 


(4.31) 


By  choosing  collocation  points  we  can  convert  this  singular  integral  equation  into  a  system 
of  linear  algebraic  equations  as 


I/"  f _ 1 _ 

7T  J-\  L  (r i  -8k) 


+  ^n(s*.n) 


,  ZAnTn(ri) 

n=l 


dr  i 


=  Pi(s*) 


AC  — 1 

4^0 


(fc  =  l,..,iV) 
(4.32) 


The  collocation  points  can  be  arbitrarily  chosen  but  good  convergences  are  obtained  if  the 
points  are  symmetric  about  the  center  of  the  crack  with  dense  distribution  near  the  ends  of 
the  crack.  The  roots  of  Chebyshev  polynomials  are  used  as  the  collocation  points.  Thus, 

Sl  =  (2fc2JV~1)7r  C*  =  l.-.JV).  (4.33) 

The  Cauchy  part  of  this  equation  is  evaluated  in  closed  form  while  the  Fredholm  part  is 
evaluated  numerically.  In  the  numerical  evaluation  of  the  Fredholm  kernel  Gauss 
Quadrature  method  is  used.  It  should  be  noted  that  the  Fredholm  kernel  contains  a 
logarithmic  singularity  of  the  form 
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i  m  .  J2-A-nTn(ri) 

-/_ i [logh 

Considering  the  numerical  difficulties  involved  in  this  integral  (eg.  Joseph  [56] ),  it  is 
evaluated  in  closed  form.  The  closed  form  expressions  for  some  of  these  integrals  are 
given  in  Appendix  D.  Integrals  of  the  type 

pi  l) 

-  J  (ri  -  Si)m  [login  -  *i|  r2 . dri 

where  m  is  a  positive  integer  emerge  from  the  Fredholm  kernel  too.  But  these  can  be 
easily  evaluated  numerically  by  splitting  the  integral  from  —  1  to  si  and  si  to  1  [41].  Thus 
by  evaluating  the  integrals,  from  (4.31)  we  obtain  an  (N  x  N)  system  of  linear  algebraic 
equations  which  can  be  solved  by  matrix  inversion.  The  right  hand  sides  depend  on  the 
specific  loading  conditions.  The  expressions  for  Pi(ai)  for  different  loading  conditions  are 
given  in  Appendix  C.  Thus  the  solution  to  the  unknown  coefficients  is  obtained.  Some 
useful  expressions  used  in  evaluating  these  integrals  are  given  in  Appendix  D. 


and 


M- 

01) 

4/ixoe  ^ 
AC  +  1 

(4.37) 

where 

N 

N 

*i(- 

1) 

=  ZAnTn(- 

1)  =  EA,(-1)* 

(4.38) 

n=  1 

and 

N 

JV 

fi(i) 

= 

£a,t„(i)  = 

(4.39) 

n=  1  n= 1 


4.2.2  Internal  crack  in  the  homogeneous  material  only 

The  relevant  singular  integral  equation  can  be  written  as 

“  /  7T~~  \ dt2  +  “  /  92(t2)Kv(x2,t2)dh  = 

nJat  (*2  ~  x2j  ^Jai  4/i0 

02  <  X2  <  62.  (4.40) 


Normalizing  (4.40)  using  (4.5)  and  (4.6),  equation  (4.40)  can  be  written  as 


32  <  1|,  (4.41) 

where 


if  \—L 

kJ-\  L(r2  - 


—  +i^22(s2,r2)J?2(r2)^2  =P2(s2)-^“; 


£>  (bi-a2  -  a2  62+02  62-02  62+02^  fA  ir^ 

K 22(32,  r2)  =( — )A22( — o — s2 - o — > — o — r2 - o — ).  (4-42) 


^  /  \  /62  -  02  62+02. 

P2 v  2)  =#>( - ~ - r2 - - - ). 


(4.43) 


From  Chapter  3  we  find  that  a2  =  /%  =  5.  Hence  (4.16)  may  be  expressed  as 
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(4.44) 


?2(r2)  = 


F2(r2) 


where  F2(r2)  is  a  bounded  function.  The  numerical  procedure  in  this  case  is  now  identical 
to  that  in  Section  4.2.1.  By  expanding  F2(r2)  in  terms  of  Chebyshev  polynomials  we  can 
now  write  (4.41)  as 


-t\— 

7^-1  L(r2  - 


M 

EFnTn(r2) 


S2) 


+  -^22  (S2,r2) 


1 


dr-2  =  P!(S2)"~-. l»2  <  1|.  (4.45) 


”*  4W 

Here  we  have  used  the  single  valuedness  condition  to  find  that  Bo  =  0.  Choosing  the 
collocation  points  as  in  (4.33),  we  obtain  the  following  system  of  algebraic  equations  : 


i/‘r  1 

nj-i  l(t*2  -  sk) 


+ 


Kv(sk,r2)] 


M 

YjBnTn{r2) 

n= 1 

v/T^ 


dr2 


=  P2(sfc) 


K  +  1 

4/io 


(fc  =  1..M) 


(4.46) 

After  evaluating  the  integrals  as  in  Section  4.2.1,  this  (M  x  M)  system  can  now  be  solved 
by  matrix  inversion. 


Stress  intensity  factors: 

Using  expressions  for  stress  intensity  factors  (3.93),  (3.94)  from  Chapter  3  and  the 
normalization  in  (4.5),  (4.6)  we  obtain  the  expressions  for  stress  intensity  factors  as 


klM-K  +  l^J  2  ^(-1) 

(4.47) 

and 

“W-.  +  i/T'W). 

(4.48) 

where 

J=2(  - 1)  =  EB„r„(  - 1)  =  £s„(  - 1)" 

n=l  n=l 

(4.49) 
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and 


*2(1)  = 


N 


N 


£>„r„(i)  =  £>„ 

n= 1  n= 1 


(4.50) 


4.2.3  Internal  cracks  in  homogeneous  and  nonhomogeneous  materials 

Using  (4.3  —  4.6),  singular  integral  equations  for  this  case  can  be  written  in 
normalized  form  as 


*/J(^o+i?u(si’ri) 


1  ft  ^ 

di(Ti)drl+-  Ki2(shr2)g2(r2)dr2 
71 J- 1 


j  /*l  ^  |  H  ^ 

“ +  - J  j^— — —  +  AWsj.wj  gjr^dn 


(4.51) 


where 


=  P2WK4+1,N<i|. 

(4.52) 

Kij(su  Tj)  —  (  2 

‘ )  tj)  (i,  j  =  1,2), 

(4.53) 

c? 

II 

(7  =  1,2), 

(4.54) 

-  /  \  pi(*i) 

—  e fa  ’ 

(4.55) 

P2  (s2 )  =J>2(«2)- 

(4.56) 
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Since  the  singularities  in  this  case  are  identical  to  those  in  Section  4.2.1  and  Section  4.2.2, 
the  numerical  procedure  is  identical.  We  find  that  we  obtain  a  system  of  (N  x  M) 
algebraic  equations  as  follows , 


N 


1  rl  *  l) 

+  *u(s*'r,)]  "TTTf  dri 


+ 


-['ft  win,  r2) 
*J- 1 


M 

£BnT„( r2) 

n=l  j..  «  ,«  +  l 


(fc  =  l,..,iV) 


(4.57) 


^  7-1 


E^nTnCrx) 

n=l 


dr! 


(r2  ~  SJk) 


+  K<&(s}.,r<2) 


M 

^SnTn(r2) 

Ti==1  -  =P2(sk)—, -  (fc  =  l,..,M). 


vi3 


4^o 


(4.58) 


The  collocation  points  are  again  chosen  as  roots  of  Chebyshev  polynomials  as  in  (4.33). 
The  resulting  system  of  equations  can  be  solved  by  matrix  inversion. 


Stress  intensity  factors: 

The  expressions  for  the  stress  intensity  factors  at  the  4  irregular  points  can  be 
derived  in  an  identical  manner  to  Section  4.2.1  and  Section  4.2.2  and  are  found  to  be 


,  {  t  v*  4//oe  ^  [b[  —  ai  f 

W  -  >i>  =  -  D". 

n=l 


(4.59) 


h(  -  oi) 


4^oe 

K  +  1 


(4.60) 
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(4.62) 


kith)  = 


4/jq 

K  -\-  1 


4.2.4  Edge  crack  in  the  nonhomogeneous  material 

From  (3.97)  the  singular  integral  equation  for  this  case  is  as  follows: 


1 

nJ-h 


01  (*i) 

(*i  ~  *i) 


dt\  + 


1  />_ai 

nj-h  j 

«  + 1 


=  Pi(*i) 


4/ioe^11  ’ 

Normalizing  the  interval  ( —  h\,  —  oi)  by  using, 


h\  <  x\  <  —  ai, 


(4.63) 


hi  —  a\  h\  +  a\ 
h  =  — o — n - ^ — 


(4.64) 


h\  —  d\  h\  d\ 

x'  =  ^-8'--r- 

the  integral  equation  can  be  written  as 


i  /’[ — 

L(n 


+  -K'n(si,n)l?i(r1)dr1  =  Pi(si)^r^'  lsi  <  !|» 


(4.65) 


(4.66) 


si)  '  -j  —  —  —  4/x0 

From  (3.109)  and  (3.110)  we  see  that,  ai  =  0,/?i  =  The  density  function  can  then  be 
written  as 

Fi(n) 


9i(ri)  = 


V1  -  n 

where  F\(ri)  is  a  bounded  function.  The  corresponding  orthogonal  polynomial  for  the 


(4.67) 


weight  function 


y/l-ri 


(0  -1) 

is  the  Jacobi  polynomial  Pn  ’  2  (r  i ) .  The  use  of  Jacobi 


polynomials  is  quite  cumbersome,  though.  Even  the  closed  form  solutions  involving  Jacobi 
polynomials  contain  serieses  that  have  to  be  summed  up  accurately.  To  avoid  this 
inconvenience,  the  orthogonal  polynomial  used  here  is  again  the  Chebyshev  polynomial  of 
the  first  kind.  Thus, 


?iOr)  = 


J2AnTn(r  i) 

n=0 


V1  ri 

Substituting  into  (4.66)  we  obtain, 


(4.68) 


-C\— 

nj- iL(n  - 


Si) 


+  Kn  (si,n) 


^2AnTn(r-\) 

71=0  drt  =  |si  <  1|, 


V1  -  ri 


(4.69) 


In  case  of  the  edge  crack  there  is  no  single  valuedness  condition.  Hence  there  are  N  +  1 
unknowns  Aq..A^.  We  need  to  choose  N  + 1  collocation  points  as  follows: 


Sk  = 


(2k  —  1)7T 


(k  =  l,..,N  +  l). 


2  (N  + 1) 

We  then  obtain  an  (N  +  1)  x  (N  +  1)  system  of  algebraic  equations  as  follows: 


(4.70) 


N 


1  />i  r  I  ,  ^2AnTn(r\ ) 

(•*.  n)]  ■  *•,  =  : K—  (t  - 1.  « + !)• 


r  i 


4/io 


(4.71) 


It  should  be  noted  here  that  because  of  choosing  Tn(ri)  as  the  orthogonal  polynomial  we 
encounter  integrals  of  the  type, 


I  f1  i  r.(n) 
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We  do  not  have  closed  form  expressions  for  these  integrals.  In  this  case  we  use  some 
algebraic  manipulations  to  get  around  this  problem  as  follows, 


1  r 1 _ 1^ 

*J-i  (n  - 


T„(ri) 


sk)  y/l  ~  H 


dri 


Tn(ri)  -  Tn(sk ) 


sk )  y/l-n 


dr  i 


+ 


Tn(Sk)  r1 _ 1_ 

n  7-i  (n  - 


=  dri 


(4.72) 


3k)  y/l-  n 

The  second  integral  on  the  right  hand  side  has  a  closed  form  solution  and  is  given  in 
Appendix  D.  The  first  integral  on  the  right  hand  side  is  a  smooth  integral  and  can  be 
evaluated  numerically.  The  logarithmic  singularity  coming  out  of  the  Fredholm  kernel  can 
also  be  similarly  handled. 


Stress  intensity  factor: 

The  only  irregular  point  is  x  =  —  a\ .  The  expression  for  the  stress  intensity  factor 
is  found  in  (3.112)  as 


Ki  i"  -L 

which  can  be  written  as 


fci(-ai) 


4/ioe  /3ai 
k  + 1 


_  N 

y/h 1  -  a^An. 

nr=0 


(4.73) 


(4.74) 


4.2.5  Crack  in  the  nonhomogeneous  material  terminating  at  interface 

The  singular  integral  equation  for  this  case  is  as  follows: 


1  f°  9l(ti) 
nj_bl(tl-xl)  1 


=  piM 


K  +  1 

Afioefa'-  ’ 


-  bi  <  x\  <  0 


(4.75) 
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Normalizing  the  interval  (  —  b\,  0)  by  using 


*i  =  — (ri-1),  *i  =  y(«i-l) 

we  obtain, 


(4.76) 


i  A- 1 

nj- iL(n  - 


«i) 


+  -Kn(si,n)  ?i(ri)^rt 


(4.77) 


From  (3.114)  we  find  that  c*i  =  /?i  =  |.  This  case  is  identical  to  that  treated  in  Section 

4.2.1.  The  same  numerical  analysis  is  followed  where  a\  is  now  replaced  by  0.  The 

resulting  system  of  algebraic  equations  is  as  follows 

N 


1  /*[_!_ 

nj- 1  L(r i  -  Sfc) 


+  -^ii(sfc*n) 


^>lnTn(ri) 

7T=1 

v/l-r? 


dri 


K  -I-  1 


(4.78) 


Stress  intensity  factors: 

By  an  analysis  similar  to  that  in  Section  4.2.1  we  find  the  stress  intensity  factors  as 


(4.79) 


fci(0)  = 


4/jp 

K  +  1 


(4.80) 


4.2.6  Crack  in  the  homogeneous  material  terminating  at  interface 

The  singular  integral  equation  for  this  case  is  as  follows: 
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“  [  7T~~{dt^  +  “  /  9‘l{t‘l)K22(x2,t2)dt2  =  P2 (^2 ) ~7~~—  > 

7T70  (52  -  ^2)  TTJo  4/io 


Normalizing  the  interval  (0,62)  by  using 


a2  <  x2  <  b2. 


(4.81) 


t2  =  y(l  +  r2),  ^2  =  y(l  +  s2) 

we  obtain 


(4.82) 


iA— L 

7T  y  _1  L  (7*2  - 


—  +i^22(s2,r2)j?2(r2)^2  =^2(s2)-^-, 


S2  <  1  . 


(4.83) 


From  (3.115),  we  obtain  a2  =  /%  =  Hence  this  case  is  the  same  as  that  treated  in 
Section  4.2.2  where  a2  is  replaced  by  0.  The  system  of  algebraic  equations  we  obtain  can 
be  expressed  as 


-f'\— 

nj-\t(r2  - 


Sk) 


+  -K22(s/C,n) 


M 

Y;BnTn(r2) 

n=l  7 

71  ^ r<1 


V1~4 


P2(sfc)-r—  (k  =  1  v,M). 

4/io 


(4.84) 


Stress  intensity  factors: 

By  an  analysis  similar  to  the  analysis  in  Section  4.2.2  we  find  the  stress  intensity 
factors  as 


*i(0)  = 


4/jq 

K  +  1 


TV 


El 


(4.85) 


h(h)  = 


4/ip  /T 

K  +  1  V  &2 


TV 

£* 


n=l 


(4.86) 
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4.2.7  Crack  passing  through  interface 

As  discussed  in  Chapter  3  there  is  no  singularity  at  the  interface.  Hence  we  can 
write  one  singular  integral  equation  for  the  entire  crack  as, 

1  a(t)  i  fa  k  +  1 

-  /  — — -zdt  +-  g(t)K(x,t)dt  =  p(x)— - ,  -  6i  <  x  <  62.  (4.87) 

nj-b^t-x)  xj-b!  4^o 

Normalizing  the  interval  ( —  &i,  62)  hy  defining, 


62  H-  61  ,62-61 

t  =  — ^ — r+  — - — . 


62  +  61  62  —  61 

I  =  — S  +  ~2- 

we  can  write  (4.87)  as 


[<rbj  +  *(s’r)]?(r)dr = |s  <  ^ 

where 

K(s,r)  =  (*£±*1  )K(x,t), 

9(r )  =  g(t). 


[  K») 

p(s)  =  <  e@x 
\p(x) 


x  <  0 
x  >  0 


From  Chapter  3  it  can  be  seen  that  g(r)  can  be  represented  as, 


?(r)  = 


(t  +  fci)5  (62  —  f)5 


or 


(4.88) 

(4.89) 

(4.90) 

(4.91) 

(4.92) 

(4.93) 

(4.94) 
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(4.95) 


d(r) 


F(r) 
y/l  —  r 


where  F(r)  is  a  bounded  function  which  can  be  expanded  in  terms  of  Chebyshev 
polynomials  of  the  first  kind  as  follows 

F(r)  = 

n=  1 


(4.96) 


Using  the  single  valuedness  condition  and  the  orthogonality  properties  of  Chebyshev 
polynomials  we  find  that  Aq  =  0.  We  can  now  write  (4.90)  as, 


1  ri  r  1  -i  YjAnTn{?) 

_  I  r _ _  ,  t>  /  „\i  n=i 

7T  L  (r  — 


s) 


+  £(«,r) 


dr  =  p(s)~~~ ,  |s  <  1|. 
4^o 


(4.97) 


This  equation  is  identical  in  form  to  (4.27).  It  is  solved  in  an  identical  manner.  Care  has  to 

be  taken  while  evaluating  the  Fredholm  kernel.  The  Fredholm  kernel  is  defined  by  (4.45) 

and  (3.167).  Hence  the  integral  has  to  be  split  into  two  integrals  with  intervals 
b\  -  bz  b\  -  62  . 

(  U  I  L  )  (  L.  ,  L  * 


^1+^2 


'bi  +bo 


Stress  intensity  factors: 

The  stress  intensity  factors  in  this  case  are  expressed  as 


fci(-bi) 


4/uoe 


-A 


K  +  1  V  “h  &2 


9*(~b  1), 


(4.98) 


h(b2) 


4/^o 


-g*(b2). 


K  +  1  V  &1  “f"  ^2  * 

Using  the  normalization  in  (4.76  —  4.77)  we  can  write  these  expressions  as 


(4.99) 


K  +  1 


(4.100) 


n=l 
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(4.101) 


h(b2)  = 


4^o  h  +  b_2 

K  +  lV  2 


4.2.8  Edge  crack  passing  through  interface 

As  in  Section  4.2.7,  we  see  from  Chapter  3  that  we  can  write  a  single  integral 
equation  as, 


+  —  f  g(t)K(x,t)dt  =  —  hi  <  x  <  b%.  (4.102) 

4^0 


i  r  +  ;  /  g(t)K(x,t)dt  = 

ft  J  —hi  \P  ft  J -hi 

Normalizing  the  interval  ( —  h\,  b2)  by  defining, 


,  _  ^2  +  hi  h  -  hi 
2  r+  2  ’ 


(4.103) 


&2  +  h\  b2  —  hi 
we  can  write  (4.50)  as 


(4.104) 


*/J(db)  +^(s>7')]?(o^=p(s) 


« + 1 


,  8  <  1 


where 


?0)  =  $(*). 


(4.105) 


(4.106) 

(4.107) 


f  £(£)  o 

P(s)  =  < 

I  p(x)  x  >  0 


From  Chapter  3  we  find  that,  we  can  write  #(r)  as, 


(4.108) 
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(4.109) 


ff(r)  = 


sr$) 

(h-t)* 


or 


g(r) 


F(r) 


(4.110) 


y/l  ~  r 

where  F(r)  is  a  bounded  function  which  can  be  expanded  in  terms  of  an  orthogonal 
polynomial.  As  in  Section  4.2.4,  we  choose  Chebyshev  polynomial  of  the  first  kind  as  the 
orthogonal  polynomial.  Thus, 


?0) 


E  AnTn{r) 

n= 0 


v/r 


Equation  (4.105)  can  now  be  written  as 


(4.111) 


-  f  f— 

itj- ,l(r- 


*) 


+  K(a,r) 


n= 0  ,  .  AC  +  1  .  | 

dr  =p(s)— - ,|s<l| 


(4.112) 


y/l  —  r  ^  '  4^zo 

Equation  (4.112)  is  identical  in  form  to  (4.69)  and  is  solved  in  an  identical  manner.  It 
should  be  noted  here  that  there  is  no  single  valuedness  condition  in  this  case.  Hence  we 
have  N  +  1  unknowns.  Thus  we  need  to  choose  N  +  1  collocation  points.  The  points  are 
chosen  according  to  (4.70).  We  then  obtain  an  ( N  +  1)  x  (N  +  1)  system  of  linear 
algebraic  equations  which  can  be  solved  for  An  7s  by  matrix  inversion. 


Stress  intensity  factor: 

The  stress  intensity  factor  at  x  =  62  can  be  written  as 
*1(62)  =  -  ■  (4.113) 

Ki  “1  JL 

By  using  the  normalization  in  (4.103  —  4.104)  we  can  write  the  stress  intensity  factor  as 
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+  11? 


(4.114) 


CHAPTER  5 


RESULTS  AND  DISCUSSIONS 

In  Chapter  4,  the  numerical  procedure  for  solving  the  integral  equations  has  been 
explained.  This  chapter  provides  numerical  solutions  for  various  crack  configurations  and 
loading  conditions. 

The  loading  conditions  can  be  divided  into  two  main  classes : 

1 .  Mechanical  loading, 

2.  Thermal  loading. 

In  mechanical  loadings,  three  different  types  of  loadings  are  considered.  First  loading  is  a 
fixed  grip  loading.  In  this,  the  strain  =  €o  is  constant  throughout  the  section.  Second 
loading,  is  the  membrane  loading.  It  is  assumed  that  a  tensile  load  is  applied  at  the 
midplane  of  the  specimen.  Third  loading  is  a  4-point  bending  giving  a  uniform  moment  at 
the  midsection  of  the  specimen.  In  thermal  loadings,  there  are  mainly  two  cases.  In  the 
first  case,  the  specimen  is  uniformly  heated  or  cooled.  In  the  second  case,  one  surface  of 
the  specimen  is  kept  at  a  constant  temperature  while  the  other  is  heated  or  cooled.  The 
second  case  corresponds  to  the  condition,  for  example,  on  a  turbine  blade  where  one 
surface  is  exposed  to  high  temperatures  while  the  substrate  is  kept  at  a  lower  temperature 
by  some  cooling  mechanism. 

The  results  of  interest  are  the  stress  intensity  factors  at  the  crack  tips  and  the  crack 
opening  displacements.  The  stress  intensity  factors  show  the  propensity  for  a  crack  to 
grow.  Higher  stress  intensity  factors  show  a  propensity  for  the  crack  to  grow  in  that 
direction  whereas  lower  stress  intensity  factors  show  a  tendency  for  crack  arrest.  The 
crack  opening  displacements  show  the  shape  of  the  crack.  These  provide  information  on 
whether  the  crack  is  closing  or  not. 

In  this  study,  two  real  material  combinations  are  used.  They  are : 
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1.  Partially  Stabilized  Zirconia  (PSZ)  on  Rene-41  (  —  =  1.45,  (3hi  =  0.375) 

An 

2.  Partially  Stabilized  Zirconia  (PSZ)  on  Ti-6A1-4V  (—  =  0.77,  (3 hi  =  -  0.2576). 

Pi 

The  thermomechanical  properties  for  these  materials  are  given  in  Appendix  E.  Apart  from 
these  material  combinations,  some  hypothetical  material  combinations  are  used  for 
comparison  purposes.  These  combinations  are : 

1.  —  =  5.0,  (3hx  =  1.6094, 

Pi 

2.  —  =  0.5,  phx  =  -  0.6931. 

Pi 

(3h\  =  0.0001  represents  the  homogeneous  case  extremely  closely.  It  is  not  possible  to 
put  (3h\  identically  equal  to  zero  in  this  model  due  to  numerical  difficulties. 

The  results  in  this  study  are  for  plane  strain  case  unless  otherwise  mentioned.  In 
most  cases  results  for  generalized  plane  stress  case  differ  from  those  in  plane  strain  case 
only  in  the  second  or  third  significant  digit.  Also,  in  this  analysis  the  Poisson's  ratio  is 
assumed  to  be  constant  for  all  materials  ( u  =  0.3). 

5.1  Mechanical  loading 

5.1.1  Stresses  in  uncracked  specimens 

Before  the  stress  intensity  factors  and  the  crack  opening  displacements  of  cracked 
specimens  are  studied,  it  is  important  to  investigate  the  stresses  in  uncracked  specimens. 
As  mentioned  before,  the  specimens  were  subjected  to  fixed  grip  loading,  membrane 
loading  and  4-point  bending.  The  stress  profiles  for  each  of  these  loadings  are  shown  in 

/j»2 

Figures  5.1  —  5.3  respectively  for  a  thickness  ratio  of  —  =  1.  Equal  and  opposite  stresses 

h  i 

are  applied  on  the  crack  surfaces  when  studying  the  crack  problems.  Figure  5.1  shows  that 
in  the  fixed  grip  case,  the  stress  profile  in  the  FGM  layer  follows  the  stiffness  profile.  In 
the  pure  bending  case,  the  area  under  the  stress  curve  is  equal  to  zero  (Fig.  5.3).  Also,  the 
point  where  the  normal  stress  is  zero  denotes  the  location  of  the  neutral  axis.  It  can  be 
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seen  that  for  a  softer  coating,  the  neutral  axis  lies  in  the  substrate  whereas  for  a  stiffer 
coating,  it  lies  in  the  coating.  Figure  5.4  and  5.5  show  the  normal  stress  in  the  uncracked 
specimen  for  membrane  loading  and  pure  bending  respectively  when  the  substrate  to 
coating  thickness  ratio  is  5. 

5.1.2  Internal  crack  in  homogeneous  substrate  only 

This  is  the  case  where  an  existing  flaw  in  the  substrate  is  put  under  various  loading 
conditions.  We  are  interested  in  finding  out  the  direction  in  which  the  crack  will  propagate 
if  at  all  it  propagates.  The  problem  of  a  crack  in  a  homogeneous  half  plane  bonded  to  an 
FGM  half  plane  was  solved  by  Erdogan  et.  al.  [43].  The  results  obtained  in  [43]  and  those 
obtained  in  this  work  for  the  special  case  of  h\  =  h%  =  oo  are  tabulated  in  Table  5.1.  It 
can  be  seen  that  there  is  an  excellent  match  in  the  results  of  these  two  independent  works. 

Figures  5.6  —  5.8  show  the  stress  intensity  factors  for  internal  cracks  in  the 
substrate.  The  normalized  stress  intensity  factors  ( see  Appendix  C  for  details  on  the 
normalization  process )  are  plotted  against  the  normalized  crack  lengths.  The  crack  is 
allowed  to  grow  symmetrically  in  the  substrate.  It  should  be  noted  here  that  as  c  — >  0,  the 
stress  intensity  factor  goes  to  some  definite  limit.  In  fact,  the  value  of  the  normalized  stress 
intensity  factor  approaches  the  normalized  stress  in  the  uncracked  specimen  at  that  point 
as  c  — >  0.  In  other  words,  k\  (c  — ►  0)  =  aw(  —  ,0 )y/c.  These  limits  determined 

Z 

analytically  for  fixed  grip  loading  for  (3hi  =  1.6094,  0.375,  0.0001,  —  0.2576  and 
—  0.6931  are  identically  equal  to  1.0.  In  Figure  5.6,  we  can  see  that  the  numerically 
calculated  values  approach  these  analytically  determined  limit.  It  is  also  seen  from  Fig.  5.6 
that  as  the  crack  tip  &2  approaches  the  free  surface  of  the  homogeneous  strip,  the  stress 
intensity  factors  increase  very  fast.  This  indicates  that  the  crack  will  grow  towards  the  free 
surface.  In  fact,  when  62  -*  ^2,  the  stress  intensity  factor  fci(i>2)  —*■  00.  The  other  crack 
tip  <22  growing  towards  the  interface  also  shows  an  increasing  stress  intensity  factor.  But  it 
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is  much  lower  than  that  at  b% .  In  fact,  as  62  — * ►  ^2  and  02  — >  0,  k\  (02)  would  approach  a 
finite  limit  corresponding  to  the  edge  crack  in  homogeneous  material.  That  problem  has 
not  been  solved  in  this  work  as  usually  it  is  the  coating  which  undergoes  surface  cracking. 
It  can  also  be  seen  that  stress  intensity  factors  are  higher  when  the  coating  is  weaker  and 
vice  versa.  This  is  because  the  overall  stiffness  of  the  material  combination  is  lower.  Also, 
when  the  crack  is  very  small  the  effect  of  the  coating  is  negligible  and  the  stress  intensity 
factors  for  all  different  material  combinations  is  the  same  as  that  in  case  of  homogeneous 
materials.  The  increase  in  the  stress  intensity  factor  for  the  end  a<2  can  again  be  attributed 
to  overall  weakening  of  the  specimen. 

In  Fig.  5.7  the  stress  intensity  factors  for  membrane  loading  are  plotted. 
Analytically,  we  find  that  the  limits  for  c  — >  0  for  (3h\  —  1.6094,  0.375,  0.0001,  —  0.2576 
and  —  0.6931  are  0.5041,  0.5001,  0.5,  0.49998,  0.49966  respectively.  In  Figure  5.7,  we 
can  see  all  the  curves  approaching  0.5  as  c  — *  0.  It  is  peculiar  that  for  end  62,  weaker  the 
coating,  lower  is  the  stress  intensity  factor  whereas  for  end  02,  the  trend  is  exactly 
opposite.  But  it  should  be  expected  from  the  stress  distribution  for  an  uncracked  specimen 
under  membrane  loading  shown  in  Fig.  5.2.  The  stress  distribution  is  such  that  for  end 
02 the  stress  is  higher  than  end  62-  Also,  weaker  the  coating,  more  pronounced  is  the 
difference.  Hence  for  (3h\  =  1.6094,  the  stress  felt  by  end  62 is  much  lower  than  the  stress 
felt  by  62  when  /3h\  =  —  0.6931.  Similar  argument  holds  for  the  end  02.  Under  bending 
loading  (Fig.  5.8),  the  limits  as  c  — ►  0  for  f3h\  —  1.6094,  0.375,  0.0001,  —  0.2576  and 
—  0.6931  are  0.7778,  0.7607,  0.75,  0.7415  and  0.7253 respectively.  The  end  62 
experiences  a  much  higher  load  than  the  end  02.  As  a  result  the  stress  intensity  factors  for 
62  are  much  higher.  Another  factor  is  that  the  end  62  is  getting  closer  to  the  edge.  That 
pushes  the  stress  intensity  factors  even  higher.  End  a2experiences  lower  stress  (Fig.  5.3). 
Hence  the  lower  stress  intensity  factors  at  02-  The  slight  increase  at  the  end  is  the  effect  of 
overall  weakening  of  the  material.  It  should  be  noted  that  the  results  are  shown  only  up  to 
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c 

—  =  0.45  because  the  accuracy  of  the  results  was  highly  diminished  as  the  crack 

«i 


approached  the  boundary.  The  numerical  values  of  the  stress  intensity  factors  have  been 
tabulated  for  phx  =  0.375  (PSZ-Rene41)  and  f3hl  =  -  0.2576(PSZ-Ti-6A1-4V)  in 
Tables  5.2  and  5.3  respectively. 

Figures  5.9  —  5.11  show  the  crack  opening  displacements  for  the  three  different 
types  of  mechanical  loadings.  The  crack  runs  from  0.25  h\to  0.75  h\.  In  the  fixed  grip 
case  (Fig.  5.9),  there  is  hardly  any  change  in  the  crack  opening  displacement  though  it  can 
be  seen  that  weaker  the  coating,  higher  is  the  crack  opening  displacement.  In  the 
membrane  loading  (Fig.  5.10)  the  difference  is  more  pronounced.  It  can  also  be  seen  that 
the  crack  surface  is  pulled  towards  the  loading  axis  which  is  in  this  case  at  the  interface  . 
The  loading  intensity  is  higher  on  the  left  part  of  the  crack  than  the  right.  This  effect  can 
be  seen  again  in  the  bending  case.  Figure  5.11  shows  the  crack  surface  displacement  is 
now  higher  away  from  the  interface.  This  is  because  a  bending  moment  causes  higher 
stresses  farther  from  the  interface  than  closer  to  the  interface. 

Figures  5.12  —  5.14  show  the  effect  of  thickness  of  the  substrate  on  the  stress 

ho  ho 

intensity  factors.  The  thickness  ratios  are  varied  from  —  =  1.0  to  —  =  10.0.  The  case 

h\  h\ 


where 


h>2 

hi 


10.0  roughly  matches  the  semi-infinite  substrate  case.  Also,  higher 


thicknesses  for  substrates  is  closer  to  practical  situations.  From  the  plots  it  can  be  seen 
that  the  thickness  of  substrate  reduces  the  stress  intensity  factors.  In  fact,  for  —  =  10.0, 

hi 

there  is  hardly  any  change  in  the  stress  intensity  factors  as  the  crack  grows.  This  result  is 
important  as  it  provides  us  with  an  effective  technique  to  suppress  crack  growth.  If  the 
substrate  is  sufficiently  large,  the  inherent  flaws  in  it  will  be  suppressed.  The  same  effect 
can  be  seen  in  membrane  loading  and  in  pure  bending.  Figures  5.15  and  5.16  show  the 
effect  of  thickness  on  crack  opening  displacement  under  fixed  grip  and  membrane  loading. 
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In  each  case,  the  crack  opening  displacement  is  lower  when  the  substrate  thickness  is 
higher. 

5.1.3  Internal  crack  in  nonhomogeneous  coating  only 

Now  we  assume  that  there  is  a  crack  in  the  coating.  Figures  5.17  and  5.19  show 
the  effect  of  the  material  nonhomogeneity  on  the  stress  intensity  factors  of  such  a  crack.  In 
Fig.  5.17,  the  specimen  is  put  under  fixed  grip  loading.  It  should  be  noted  here  that  for  a 
positive  nonhomogeneity  value,  the  material  stiffens  towards  the  interface  whereas  for  a 
negative  nonhomogeneity  value,  the  material  softens  towards  the  interface.  Again  it  is  seen 
that  as  c  — >  0,  the  stress  intensity  factor  goes  to  some  definite  limit.  In  fact,  the  value  of 
the  normalized  stress  intensity  factor  approaches  the  normalized  stress  in  the  uncracked 
specimen  at  that  point  as  c  — >  0.  In  other  words,  k\(c  — »  0)  =  ayy(  — — ,  0 )\/c.  For 

fixed  grip  loading,  the  limits  for  /3h\  =  1.6094,  0.375,  0.0001,  -  0.2576  and  -  0.6931 
are  0.4472,  0.8290,  1.0,  1.1375  and  1.4142  respectively.  The  numerical  results  match 
closely  with  these  analytical  limits.  In  general,  it  is  seen  that  a  crack  tip  moving  into  a 
stiffening  region  has  an  increasing  stress  intensity  factor.  In  Fig.  5.17,  /3h\  —  1.6094  is 
represented  by  the  lowermost  curves.  The  stress  intensity  factor  at  tip  gi  denoted  by  the 
solid  line  increases  steadily  as  the  crack  grows.  On  the  other  hand,  for  /3hi  =  —  0.6931, 
represented  by  the  topmost  curves,  the  stress  intensity  factor  at  tip  ^decreases  as  the 
crack  grows.  As  the  crack  grows  very  large,  there  is  a  slight  increase  which  is  attributed  to 
the  overall  weakening  of  the  specimen.  These  results  are  consistent  with  previously 
published  results  [43]. 

For  the  crack  tip  b\,  there  are  two  factors  to  be  considered.  One  is  the  material 
stiffness  and  other  is  the  proximity  to  the  edge.  For  /3h\  =  1.6094,  the  stress  intensity 
factor  initially  drops  as  the  crack  tip  is  entering  softer  material.  But  as  the  crack  tip  comes 
closer  to  the  edge,  the  ligament  becomes  very  small  and  the  stress  intensity  factor  starts 
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increasing.  For  j3h\  =  -  0.6931,  both  these  factors  tend  to  increase  the  stress  intensity 

factor  and  as  a  result  the  stress  intensity  factor  rises  very  fast.  This  would  suggest  that 

there  is  a  high  propensity  for  crack  growth  in  that  direction. 
hi 

As  c  — >  — ,  the  internal  crack  approaches  an  edge  crack.  Then  it  should  be 

A 

expected  that  k\(bi)  — >  oo  and  k\ (ai)  should  approach  a  finite  value  obtained  in  the  edge 

crack  problem.  The  results  for  edge  crack  are  given  in  Section  5.1.6.  The  limiting  values 

obtained  for  fci(ai)  for  /3ht  =  1.6094,  0.375,  0.0001,  -  0.2576and  -  0.6931  are 

1.89337,  3.29216,  3.9943,  4.5893  and  5.8545  respectively.  It  is  seen  that  these  values 

are  approached  asymptotically  and  the  slope  of  the  curve  is  infinite  as  the  crack 

approaches  the  interface.  This  phenomenon  can  be  explained  from  the  fact  that  during  the 

transition  from  an  internal  crack  to  an  edge  crack  there  is  a  jump  when  the  ligament 

breaks.  This  causes  the  unbounded  slope.  The  limiting  value  is  still  bounded. 

Figure  5.19  shows  the  effect  of  material  stiffness  under  bending.  It  should  be  noted 

here  that  the  bending  is  reversed  as  compared  to  that  shown  in  Fig.  5.3.  Hence  the  plot  in 

Fig.  5.3  should  be  reflected  about  the  x-axis  to  get  the  proper  stress  distribution.  This  has 

been  plotted  in  Fig.5.18.  Again,  as  c  — ►  0,  the  stress  intensity  factor  goes  to  some  definite 

limit.  These  limits  are  obtained  for  /3h\  =  1.6094,  0.375,  0.0001,  -  0.2576  and 

—  0.6931  as  0.8544,  0.8009,  0.75,  0.7062  and  0.6169  respectively.  Numerical  results 

match  with  these  analytical  limits.  The  trends  in  Fig.  5.19  are  in  agreement  with  the 

observations  made  in  the  fixed  grip  case.  Only  difference  is  seen  in  the  plot  for 

/3/ii  =  1.6094.  The  tip  aigrows  into  a  stiffer  region  and  hence  should  have  an  increasing 

stress  intensity  factor.  From  Figure  5.19,  though,  it  can  be  seen  that  for  —  0.45,  the 

h\ 

stress  intensity  factor  actually  decreases.  This  can  be  explained  on  the  basis  of  Fig.  5.18 
where  there  is  a  decrease  in  the  stress  around  that  region.  Hence  the  stress  intensity  factor 
decreases  inspite  of  increasing  stiffness.  The  numerical  values  of  the  stress  intensity 
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factors  have  been  tabulated  for  f3h\  =  0.375  (PSZ-Rene41)  and  (3hi  =  -  0.2576(PSZ- 
Ti-6A1-4V)  in  Tables  5.4  and  5.5,  respectively. 

A  very  practical  situation  is  where  a  homogeneous  coating  is  laid  on  an  FGM 
substrate.  Under  fixed  grip  loading  it  was  observed  that  crack  propagated  beneath  the 
interface.  To  analyze  this  situation,  a  crack  was  introduced  starting  at  the  interface  and 
growing  into  the  FGM  substrate.  If  the  FGM  substrate  is  softer,  it  was  observed  that  the 
stress  intensity  factor  reduces  ( Fig.  5.20 )  whereas  if  the  FGM  substrate  is  stiffer,  the 
stress  intensity  factor  tends  to  increase  ( Fig.  5.21 ). 

The  crack  opening  displacements  for  an  internal  crack  in  the  FGM  coating  show 
some  interesting  features.  In  Fig.  5.22,  the  specimen  is  put  under  fixed  grip  loading.  It  is 
observed  that  the  crack  opening  displacements  are  lower  when  the  nonhomogeneity 
parameter  is  positive  i.e.  for  coatings  that  stiffen  toward  the  interface.  This  is  due  to  the 
fact  that  under  fixed  grip  loading,  the  load  is  distributed  in  such  a  way  that  the  softer  part 
of  the  coating  bears  lesser  load  (Fig.  5.1).  If,  on  the  other  hand,  the  nonhomogeneity 
parameter  is  negative,  the  coating  softens  towards  the  interface  and  the  load  under  fixed 
grip  loading  is  higher  for  the  coating.  As  a  result,  the  crack  opening  displacements  are 
higher.  For  all  different  nonhomogeneity  parameters,  the  crack  opening  displacements  are 
skewed  towards  or  away  from  the  interface  depending  on  the  sign  of  (3h\.  Figures  5.23 
and  5.24  show  the  crack  opening  displacements  for  membrane  loading  and  pure  bending 
respectively.  Here  again  the  curves  are  skewed  depending  on  the  sign  of  /3h\.  Also,  it  is 
observed  that  the  crack  opening  displacements  are  higher  for  a  softer  coating  than  for  a 
stiffer  coating. 

The  effect  of  thickness  ratio  on  the  stress  intensity  factors  and  the  crack  opening 
displacements  is  shown  in  Fig.  5.25  and  Fig.  5.26  respectively.  It  is  seen  that  the  stress 
intensity  factors  as  well  as  the  crack  opening  displacements  are  lower  with  increasing 
thickness  of  the  substrate.  The  difference  is  even  more  pronounced  for  the  pure  bending 
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case  ( Figures  5.27, 5.28).  Note  that  the  stress  intensity  factors  and  the  crack  opening 
h<2 

displacements  for  —  =  10.0  seem  to  be  almost  negligible.  This  is  due  to  the  fact  that  the 

hi 

stress  intensity  factors  are  normalized  with  respect  to  the  square  of  the  thickness. 

5.1.4  Internal  cracks  in  nonhomogeneous  coating  and  in  homogeneous  substrate 

In  this  section  the  interaction  effects  of  a  crack  present  in  the  homogeneous 
substrate  on  the  stress  intensity  factors  and  crack  opening  displacements  for  crack  in  the 
FGM  coating  are  studied.  First  a  crack  ranging  from  OAhi  to  0.6h\  is  placed  in  the 
homogeneous  substrate.  Then  a  longer  crack  ranging  from  0.2h\  to  0.8h\  is  placed.  The 
effect  of  these  cracks  on  a  crack  located  symmetrically  in  the  FGM  coating  is  studied  by 
comparing  the  stress  intensity  factors  in  these  two  cases  with  those  obtained  from  Section 
5.1.3  where  there  was  no  crack  in  the  substrate.  The  results  are  tabulated  in  Table  5.6. 

It  can  be  seen  from  Table  5.6  that  the  presence  of  a  crack  in  the  homogeneous 
substrate  enhances  the  stress  intensity  factors  at  tip  a\.  As  the  ligament  between  the  crack 
shortens  in  length,  the  effect  is  more  pronounced.  The  stress  intensity  factor  at  tip  feion  the 
other  hand  is  relatively  unchanged.  A  complete  table  of  stress  intensity  factors  for  the 
crack  in  the  FGM  coating  in  the  presence  of  a  crack  in  the  homogeneous  substrate  for 
different  mechanical  loadings  in  given  in  Table  5.7.  It  should  be  noted  that  the  bending  is 
reversed  so  that  the  crack  in  the  FGM  coating  is  open. 

Fig.  5.29  shows  the  effect  of  the  substrate  thickness  on  the  stress  intensity  factors. 
The  stress  intensity  factors  are  lower  for  higher  substrate  thickness.  In  Fig.  5.30,  the  stress 
intensity  factors  are  plotted  when  both  the  cracks  are  growing  symmetrically.  It  is  seen 
that  the  stress  intensity  factors  at  ends  a\  and  a 2  are  rapidly  increasing.  This  is  because  of 
the  close  interaction  of  the  two  tips.  Tips  bi  and  62  show  high  stress  intensities  when  they 
are  very  close  to  the  edges.  The  crack  opening  displacements  for  the  pure  bending  case  are 
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shown  in  Fig.  5.31.  At  higher  thicknesses  of  the  substrate  the  crack  opening  displacements 
are  negligible. 

5.1.5  Internal  crack  crossing  the  interface 

In  Section  5.1.4,  it  was  observed  that  the  stress  intensity  factors  at  tip  a\  and  ao 
tend  to  increase  as  the  ligament  size  in  between  them  decreases.  This  leads  to  the 
possibility  of  the  cracks  fusing  together  at  the  interface  forming  one  single  closed  crack 
crossing  the  interface.  The  stress  intensity  factors  for  such  an  internal  crack  are  discussed 
in  this  section. 

Figure  5.32  shows  the  stress  intensity  factors  for  a  crack  running  from  —  b\  to 
i>2 subjected  to  fixed  grip  loading.  First,  it  should  be  noted  that  as  c  — >  0,  ki(b\)  and  (62) 

both  tend  identically  (Jq^Jc  for  all  different  material  combinations.  The  crack  is  then 
allowed  to  grow  symmetrically.  When  the  coating  has  lower  stiffness  than  the  substrate 
(eg.  j3h\  =  1.6094 ),  end  b\  grows  into  a  softer  region.  In  this  case,  it  is  observed  that  the 
stress  intensity  factors  reduce.  On  the  other  hand,  when  the  coating  is  stiffer  than  the 
substrate  ( eg.  (3h\  —  -  0.6931 ),  end  61  grows  into  a  stiffer  region  and  the  stress 
intensity  factor  increases.  The  numerical  values  for  the  stress  intensity  factors  for  a 
practical  material  ( PSZ-Rene  41,  fih\  =  0.375  )  are  tabulated  in  Table  5.8. 

It  should  be  noted  here  that  the  negative  stress  intensity  factors  under  bending 
indicate  that  the  crack  tip  is  closing  under  the  loading.  The  surfaces  will  be 
interpenetrating.  A  correct  analysis  of  the  problem  should  be  performed  by  considering  a 
crack-contact  problem.  But  the  negative  stress  intensity  factors  are  still  useful  in  the  case 
of  multiple  loading,  where  under  other  applied  loadings  in  conjunction  with  the  bending, 
the  crack  is  completely  open.  The  negative  stress  intensity  factors  are  valid  only  in  such  a 
situation  where  there  is  no  closure  of  the  crack. 
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Figure  5.33  shows  the  effect  of  thickness  ratio  on  the  stress  intensity  factors  under 
reversed  bending  so  that  the  crack  is  completely  open.  The  results  are  in  accordance  with 
those  obtained  in  previous  cases.  The  crack  opening  displacements  in  case  of  membrane 
loading  are  shown  in  Figure  5.34.  These  curves  are  skewed  depending  upon  the  value  of 
fih!. 


5.1.6  Edge  crack  in  FGM  coating  and  passing  through  the  interface 

An  edge  crack  in  the  nonhomogeneous  coating  is  the  most  important  practical 
crack  configuration.  The  stress  intensity  factors  and  the  crack  opening  displacements  for 
this  crack  configuration  will  be  discussed  in  detail  in  this  section.  Also,  the  crack  will  be 
allowed  to  pass  through  the  interface  and  enter  the  homogeneous  substrate. 

The  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  strip  have  been 
tabulated  by  Joseph  [63] .  By  assuming  a  very  small  value  for  the  nonhomogeneity 
parameter  ( /3h\  =  0.0001  ),  those  results  have  been  reproduced  in  this  work.  The 
comparison  of  results  is  given  in  Table  5.9  and  5.10  for  fixed  grip  loading  and  pure 
bending  respectively.  It  can  be  seen  from  the  comparison  that  the  results  in  this  work 
match  very  well  with  published  results.  In  fact,  in  this  work  the  convergence  is  much 
better  for  very  small  ligament  lengths  than  that  in  previously  published  results. 

Figure  5.35  shows  the  stress  intensity  factors  for  an  edge  crack  in  the 
nonhomogeneous  coating  under  fixed  grip  loading.  It  should  be  noted  that  as  c  — ►  0,  the 
stress  intensity  factors  approach  a  definite  limit  which  is  analytically  found  to  be  1.121522 
times  cFyy{  —  h\,Qi)yfc  obtained  from  Koiter's  integral.  The  numerical  limits  approach 
extremely  closely  to  this  analytical  limit  as  seen  in  Figure  5.35.  The  crack  is  then  allowed 
to  reach  up  to  the  interface.  The  effect  of  the  nonhomogeneity  parameter  is  studied.  It  can 
be  seen  that  an  edge  crack  growing  into  a  stiffer  region  (  j3h\  >  0 ),  shows  increasing 
stress  intensity  factors.  For  (3h\  =  1.6094,  the  increase  is  very  rapid.  On  the  other  hand, 
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an  edge  crack  growing  into  a  softer  region  {(3h\  <  0),  shows  a  slight  decrease  before 

increasing  very  slowly.  For  a  homogeneous  medium  (  (3h\  =  0.0001),  the  slope  of  the 

c 

curve  for  - - ►  0,  is  zero.  For  a  nonhomogeneous  material,  whether  it  is  stiffening  or 

hi 

softening,  the  same  should  hold  true.  As  the  crack  is  allowed  to  grow  into  the 
homogeneous  substrate,  the  trends  in  the  stress  intensity  factors  continue.  These  are 
shown  in  Figure  5.36.  The  rapid  increase  in  the  stress  intensity  for  f3h\  =  1.6094  may 
lead  to  catastrophic  failure.  It  should  be  noted  that  the  stress  intensity  factor  curves  are 
continuous  at  the  interface.  This  is  because  there  is  no  singularity  at  the  interface  as  shown 
in  Chapter  3.  A  table  of  stress  intensity  factors  for  different  loadings  for  practical  materials 
is  given  in  Table  5.11  —  5.14. 

Figure  5.37  —  5.40  show  the  variation  of  stress  intensity  factors  for  edge  cracks 
under  membrane  loading  and  pure  bending.  Again  it  is  seen  analytically  that  as  c  — ►  0, 
ki(ai)  -»  0.3037,  0.5042,  0.5608,  0.5942  and  0.6344 for  phx  =  1.6094,  0.375, 

0.0001,  —  0.2576  and  —  0.6931  respectively  for  membrane  loading.  Also, 
fci(ai)  — >  0.7301,  1.41,  1.6823,  1.8893  and  2.2809  respectively  for  pure  bending.  The 
numerical  results  match  with  these  analytical  limits  very  well.  It  can  be  seen  from  Fig.  5.38 
and  Fig.  5.40  that  once  the  crack  enters  the  homogeneous  substrate,  the  variation  in  the 
stress  intensity  factors  with  the  nonhomogeneity  parameter  is  almost  negligible.  This 
suggests  that  under  these  loadings  the  influence  of  the  nonhomogeneity  of  the  coating  is 
minimal  especially  if  the  crack  grows  into  the  substrate. 

Figure  5.41  —  5.44  show  the  effect  of  thickness  ratio  on  the  stress  intensity  factor. 
In  Fig.5.41,  it  can  be  seen  that  the  stress  intensity  factors  are  much  lower  for  higher 
substrate  thicknesses.  From  Fig.  5.42,  it  can  be  seen  that  as  the  crack  tip  reaches  the  edge, 
the  stress  intensity  factors  rise  very  rapidly.  The  numerical  values  for  the  stress  intensity 
factors  for  PSZ-Rene41  for  ~  =  10.0  are  tabulated  in  Table  5.15  —  5.16. 
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Fig.  5.43  and  Fig  5.44  show  the  effect  of  thickness  for  bending  loading.  It  can  be 
seen  that  for  smaller  thicknesses  (  7^  =  0. 1 ),  the  stress  intensity  factors  are  very  high 

hi 

but  for  higher  thicknesses  (  7^  =  10.0 ),  the  stress  intensity  factors  are  very  low.  Again 

hi 

this  is  due  to  the  normalization  used  for  bending  loading.  As  the  crack  passes  through  the 
interface,  the  stress  intensity  factors  are  still  very  low  for  higher  thicknesses. 

The  crack  opening  displacements  for  varying  nonhomogeneity  parameters  are 
plotted  in  Fig.  5.45  -  5.47.  In  Figure  5.45,  an  edge  crack  in  the  nonhomogeneous  coating 
is  considered.  The  crack  opening  displacement  for  softer  coating  ( eg.  (3h\  =  1.6094 )  is 
lowest  because  under  fixed  grip  loading  the  stresses  in  uncracked  specimen  are  lower  for 
softer  coatings  (Fig.  5.1).  The  same  trend  is  seen  for  a  crack  crossing  the  interface  ( Fig. 
5.46).  In  membrane  loading  (Fig.  5.47),  however,  the  stresses  are  higher  for  softer 
coatings  ( Fig.  5.2)  and  we  see  the  expected  result  that  the  crack  opening  displacements 
are  higher  for  softer  coatings.  The  effect  of  stiffness  is  studied  in  Fig.  5.48  —  5.49.  As  the 
thickness  of  the  substrate  increases,  the  crack  opening  displacements  are  reduced. 

5.2  Thermal  loading 

5.2.1  Stresses  in  uncracked  specimens 

The  uncracked  thermal  stress  problem  needs  to  be  solved  for  the  specimen 
geometry.  Because  this  is  a  one-dimensional  heat  conduction  problem  the  thermal  analysis 
is  very  straightforward  and  the  temperature  distribution  and  resulting  thermal  stresses  can 
be  obtained  directly.  If  on  the  other  hand,  if  the  problem  was  two-dimensional,  the  crack 
will  introduce  partial  insulation.  In  that  case  the  diffusion  problem  as  well  as  the  stress 
problem  will  be  analytically  intractable.  The  problem  of  insulation  due  to  the  crack  has 
been  solved  by  Finite  Element  Analysis  [69]. 
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In  the  current  problem,  we  assume  that  the  thermal  and  mechanical  problems  are 
uncoupled.  Hence  it  is  only  necessary  to  solve  the  diffusion  problem  and  the  thermal  stress 
problem  for  an  uncracked  specimen.  The  thermal  stresses  thus  obtained  are  then  applied  to 
the  crack  geometry  to  solve  the  elasticity  problem.  It  should  be  noted  that  the  thermo¬ 
mechanical  properties  of  the  FGM  coating  are  functions  of  the  space  coordinate  in  the 
thickness  direction. 

Temperature  distribution 

The  diffusion  equation  for  the  geometry  shown  in  Figure  5.50  is 

9  rut  \9Ts  ,  v  ,  ,dT 

ei{Hx)fc}  =  p(x  )c{x)m  (5-1) 

where  T(x,t )  is  the  temperature,  k(x)  is  the  coefficient  of  heat  conduction,  p(x) is  the 
density  and  c(x)  is  the  specific  heat.  For  a  steady  state  heat  conduction  problem  (5.1)  can 
be  written  as 


,,  sdT 

«(x)—  =  constant. 


(5.2) 


The  form  of  k(x)  can  be  arbitrary.  But  the  elasticity  problem  can  be  solved  analytically 
only  for  some  specific  forms  of  variations  of  the  elasticity  modulus.  In  this  work,  the 
exponential  variation  of  the  elasticity  modulus  is  considered.  Hence  to  be  consistent,  the 
coefficient  of  heat  conduction  is  assumed  to  be  of  the  form 

k{x)  =  (5.3) 

where  ko  is  the  coefficient  of  heat  conduction  of  the  homogeneous  substrate  and  7  is  the 
nonhomogeneity  parameter  for  heat  conduction  defined  as 


(5.4) 
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where  kc  is  the  coefficient  of  heat  conduction  of  the  coating  and  fcois  the  coefficient  of 
heat  conduction  of  the  substrate.  If  the  temperatures  are  specified  as  T\  and  at  the 
boundaries  as  shown  in  Figure  5.50,  the  temperature  distributions  can  be  expressed  as 

T*h  =  ~  +Ti’  <5-5) 


Tk  =  (T2-Ti)f  +Ti 
ri2 

where  Ti  is  the  temperature  at  the  interface  given  as 


(5.6) 


Ti  = 


Try 

e'lh  -  i 


7 

efl*»  - 1 


(5.7) 


It  should  be  noted  that  for  uniform  temperature  change  from  some  reference  temperature 
To  to  T,  Tnh  =  Th  =  Ti  =  T.  It  is  more  convenient  to  deal  with  temperature  change 
rather  than  absolute  temperatures.  In  this  case  AT  =  T  -  To  =  constant  throughout  the 
specimen.  If,  on  the  other  hand,  we  choose  to  hold  the  temperature  T%  =  To  and  vary  T\ , 
the  driving  force  is  the  potential  ATi  =  Ti  -  To.  We  can  normalize  the  temperature 
distribution  with  respect  to  this  potential  ATi  as 


AT„a  — 


(1  -  A Tt-)e^a  -  1  +  A 
e^i  -  1 


(5.8) 


ATk  =  A7j  (-\— )•  (5-9) 

where 

^=(^)/(^+y.  (5.10) 

The  temperature  distributions  thus  obtained  can  now  be  used  to  find  the  thermal  stresses 
in  the  uncracked  specimen. 
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Thermal  stress  distribution 


For  an  unconstrained  beam,  the  strain  can  be  written  as  [68] 

Cyy  —  —  a(x)AT(x)  +  Ax  +  B  (5-11) 

where  a(x )  is  the  coefficient  of  thermal  expansion  of  the  FGM  coating  assumed  to  be  of 
the  form  of 

a(x)  --  aoe0*  (5.12) 

where  »o  is  the  coefficient  of  linear  expansion  of  the  homogeneous  substrate  and  a  is  the 
nonhomogeneity  parameter  defined  as 

a  =  rl09&  (5-13) 

h\  ao 

where  afis  the  coefficient  of  linear  expansion  of  the  ceramic.  The  stress  can  then  be 
written  as 

ayy  =  y^(-&(x)AT(x)  +  Ax  +  B)  (5-14) 


where  k  =  (3  —  u)/(l  +  v)  for  generalized  plane  stress.  For  plane  strain  k  =  3  —  and 
a  should  be  replaced  by  2(1  +  u).  A  and  B  in  (5.14)  are  evaluated  from  equilibrium 
conditions 


c 

H 

II 

O 

(5.15) 

th. 

/  XGyy{x)dX  =  0. 

J-h , 

(5.16) 

A  and  B  are  obtained  as 

A  R1C4  —  i?2c2 

^  ) 
clc4  —  C3C4 

(5.17) 

g  R%c\  -  R1C3 

C1C4  —  C3C4 

(5.18) 

where  c*  (i  =  1..4)  and  R\,  R%  are  given  in  Appendix  C.  The  thermal  stress  in  uncracked 
specimens  under  uniform  temperature  change  and  under  steady  state  heat  conduction  is 
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shown  in  Figures  5.51  to  5.56.  It  should  be  noted  that  the  plots  are  for  positive  AT.  If 
AT  is  negative,  the  plots  should  be  just  reversed  about  the  zero  stress  axis  (i.e.  stresses 
are  exactly  equal  and  opposite).  This  is  important  when  determining  if  a  surface  crack  will 
open  or  close  under  a  given  loading. 

It  can  be  seen  from  Figures  5.51-5.56  that  the  stress  profiles  are  such  that  the  area 
under  the  curve  is  zero.  This  is  indeed  expected  as  pure  thermal  loading  is  self- 
equilibriating  since  there  is  no  external  load  acting  on  the  specimen.  The  consequence  is 
that  the  loading  changes  from  compression  to  tension  and  vice-versa.  This  represents  a 
peculiar  situation.  The  loading  maybe  compressive  on  part  of  the  crack  while  it  could  be 
tensile  on  the  rest  of  the  crack.  This  could  lead  to  smooth  crack  closing.  For  example,  for 
an  edge  crack,  if  the  loading  is  tensile  on  the  surface,  the  crack  tends  to  grow  until  it 
reaches  the  compressive  zone  where  it  gets  arrested.  On  the  other  hand,  if  the  loading  is 
compressive  on  the  surface,  it  closes  the  crack  at  its  mouth.  A  long  crack,  however,  might 
open  up  in  the  interior  region  under  tensile  loading.  This  leads  to  the  crack-contact 
problem  which  needs  to  be  looked  into  carefully. 

5.2.2  Crack-Contact  problem 

A  possible  crack-contact  problem  is  shown  in  Figure  5.56.  We  can  see  here  that 
the  crack  is  closed  from  x  =  —  h\ to x  =  —  a.  It  is  open,  however  from  to  a:  =  —  a  to 
x  =  —  ai.  Atx  =  —  a,  the  crack  will  form  a  cusp  i.e.  there  is  smooth  closure  of  the 
crack.  The  stress  intensity  factor  atx  =  -  a  will  be  zero.  At  x  =  -  ,  though,  the 

crack  is  open  and  the  stress  intensity  factor  is  positive. 

This  problem  is  again  solved  by  a  superposition  technique.  The  thermal  problem  is 
first  solved  as  described  in  Section  5.2.1.  Equal  and  opposite  stresses  are  applied  on  the 
crack  face.  If  the  stresses  are  compressive  along  the  entire  length  of  the  crack,  the  stress 
intensity  factor  at  the  crack  tip  will  be  negative.  This  indicates  a  complete  closing  of  the 
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crack.  If  the  crack  is  deep  enough  so  that  the  crack  tip  is  in  the  tensile  zone,  the  stress 
intensity  factor  at  the  crack  tip  will  be  positive  but  there  will  be  interpenetrating  crack 
opening  displacements  indicating  that  the  crack  is  closing  at  some  point.  This  is  the  case 
shown  in  Figure  5.57.  The  point  where  the  crack  closes  (x  =  —  a)  is  unknown  but  can  be 
found  out  by  satisfying  the  condition 

^(-^=0,  (5.19) 

where  k\  is  the  Mode  I  stress  intensity  factor.  Equation  (5.17)  is  highly  nonlinear  in  a.  We 
employ  a  simple  iterative  procedure  to  determine  a.  We  first  introduce  an  internal  crack 
instead  of  an  edge  crack.  Reader  is  referred  to  Chapter  4  for  details  of  the  internal  crack 
problem.  The  crack  tip  a  is  varied  to  determine  the  exact  location  where  (5.17)  is  satisfied. 
Usually  a  good  four  significant  digit  accuracy  is  obtained  within  five  iterations. 

5.2.3  Surface  cracks  and  cracks  passing  through  the  interface 

The  important  results  in  this  fracture  characterization  are  the  stress  intensity 
factors  under  given  loadings.  Two  different  kinds  of  loadings  are  considered.  In  one 
loading,  the  temperature  of  the  entire  specimen  is  raised  ( or  lowered )  uniformly.  In  this 
case,  there  is  no  conduction.  In  another  type  of  loading,  the  temperature  of  the  substrate  is 
maintained  at  a  fixed  temperature,  whereas  the  surface  of  the  coating  is  raised  ( or 
lowered)  to  a  higher  temperature.  The  latter  case  in  which  the  surface  of  the  coating  is 
heated  while  the  substrate  is  maintained  at  a  lower  temperature  is  observed  in  thermal 
barrier  coatings.  The  stress  intensity  factors  are  evaluated  under  each  of  these  loadings  for 
two  different  practical  material  combinations : 

I.  Partially  Stabilized  Zirconia  (PSZ)  on  Rene-41  (  —  =  1.45,  (3h\  =  0.375) 

Fl 

II. Partially  Stabilized  Zirconia  (PSZ)  on  Ti-6A1-4V  (—  =  0.77, 6h\  =  -  0.2576). 

Pi 
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First  we  consider  Material  combination  I.  Figure  5.58  shows  the  stress  intensity 

factors  under  uniform  temperature  rise.  It  can  be  seen  from  Figure  5.51  that  the  thermal 

stresses  are  tensile  near  the  surface.  At  x  =  —  0.4073,  the  stresses  become  compressive. 

The  crack  initially  grows  under  the  tensile  loading  which  can  be  seen  by  the  positive  stress 

intensity  factors.  As  the  loading  becomes  compressive,  the  stress  intensity  factors  drop  and 

almost  go  to  zero  as  the  crack  approaches  the  interface.  It  can  clearly  be  seen  that  the 

crack  will  soon  be  arrested.  In  fact,  we  can  see  in  Figure  5.59,  that  the  stresses  do  in  fact 

go  to  zero  when  the  crack  passes  through  the  interface.  The  numerical  values  of  these 

stress  intensity  factors  are  tabulated  in  Table  5.17.  Next,  the  same  edge  crack  has  been 

subjected  to  heat  conduction  ( Figure  5.52).  Here  the  stresses  are  compressive  at  the 

surface  and  tensile  within.  Hence  we  have  a  crack-contact  problem.  It  can  be  seen  from 

Figure  5.52  that  the  stresses  change  from  compressive  to  tensile  at  x  =  —  0.6873.  Hence 

it  is  expected  that  cracks  within  the  compressive  zone  will  tend  to  remain  closed.  This  is 

seen  in  Figure  5.60  which  shows  the  point  of  closure  of  the  crack  as  the  crack  tip 

x  =  —  a\  grows  towards  the  interface.  If  the  normalized  crack  tip  distance  from  the 
h\  —  a\ . 

surface  ( — - )  is  smaller  than  0.3127,  both  the  crack  tips  coincide.  In  other  words,  the 

h\ 

crack  is  completely  closed.  If,  on  the  other  hand,  it  is  larger  than  0.3127,  the  tip  where  the 
crack  closes  (crack  tip  b\  in  this  case)  moves  towards  the  surface.  From  the  plot  it  can  be 
seen  that  crack  tip  bi  moves  closer  and  closer  to  the  surface  and  is  expected  to  reach  the 
surface  when  the  open  crack  tip  cuts  completely  through  the  entire  thickness  of  the 
specimen.  This  can  be  seen  in  Figure  5.61.  As  the  open  crack  tip  reaches  the  surface 
x  =  h.2,  b i  reaches  the  surface  x  =  —  hi.  Figures  5.62  and  5.63  show  the  stress  intensity 
factors  for  the  conduction  case  for  a  surface  crack  in  the  coating  and  crack  crossing  the 
interface  respectively.  When  the  normalized  crack  tip  distance  is  less  than  0.3127,  the 
stress  intensity  factors  are  zero  as  the  crack  is  completely  closed.  For  a  longer  crack,  the 

stress  intensity  factors  increase  initially  but  as  the  loading  drops,  the  stress  intensity  factors 
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drop.  In  fact,  the  stress  intensity  factors  eventually  go  to  zero  as  seen  in  Figure  5.63.  As 
the  loading  is  self-equilibrating,  this  is  an  expected  result.  It  should  be  noted  here  that  as 
the  open  crack  tip  reaches  x  =  /12,  the  numerical  accuracy  drops  rapidly  due  to 
convergence  problems.  Best  values  were  obtained  when  the  open  crack  tip  was  at  least  0.1 
hi  distance  away  from  the  free  surface. 

When  the  loadings  considered  in  the  two  previous  cases  are  reversed,  we  see  that 
the  compressive  loadings  become  tensile  and  vice  versa.  Consider  the  case  of  uniform 
cooling  ( Figures  5.64  -  5.65).  Because  the  loading  is  a  mirror  image  of  that  in  Figure 
5.51,  the  stresses  are  compressive  near  the  surface  and  tensile  within.  Hence  there  will  be 
crack-contact.  The  critical  point  where  the  stresses  become  tensile  is  x  —  —  0.4073.  In 
Figure  5.64,  we  can  see  that  the  two  crack  tips  coincide  before  this  point.  After  that  point, 
the  crack  tips  separate,  the  crack  tip  b\  slowly  moving  toward  the  surface.  Figure  5.65 
shows  the  stress  intensity  factors  for  this  case.  The  open  crack  tip  in  this  case  is  only 
allowed  to  reach  the  interface  in  this  case  but  it  can  be  easily  let  to  cross  the  interface 
following  the  same  procedure  as  in  the  conduction  case  discussed  before.  If  we  now  hold 
T2  =  To,  and  let  T\  decrease,  the  stress  intensity  factors  are  tensile  near  the  surface.  The 
crack  grows  initially  but  as  it  grows  into  the  compressive  zone,  it  gets  arrested  (Figure 
5.66). 

The  same  kind  of  analysis  is  done  for  Material  combination  H  Figures  5.53  —  5.54 
show  the  uncracked  stresses.  Consider  the  uniform  heating  case.  The  stresses  from  Figure 
5.53  show  that  this  indeed  will  be  a  crack-contact  problem.  Figure  5.67  shows  the  location 
of  the  closed  crack  tip.  The  critical  point  is  x  =  —  0.532.  The  stress  intensity  factors  in 
Figure  5.68  show  correspondingly  that  stress  intensity  factors  are  zero  when  normalized 
crack  tip  distance  is  less  than  0.468.  In  the  case  of  uniform  cooling,  the  crack  is  open  at 
the  surface.  It  will  grow  under  the  tensile  loading  but  will  be  arrested  in  the  compressive 
zone  (Figure  5.69).  Figures  5.70  —  5.71,  show  the  crack-closure  lengths  and  stress 
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intensity  factors  respectively  for  surface  heating.  Figure  5.72  shows  the  stress  intensity 
factor  for  surface  cooling  case.  The  results  in  Figure  5.72  can  now  be  compared  with  those 
obtained  in  Figure  5.66  as  the  loading  conditions  are  identical.  That  comparison  is  showed 
in  Figure  5.73. 

In  all  the  cases  considered  above,  the  thickness  ratio  was  7^  =  1.0.  In  practical 

hi 

situations,  the  substrate  is  usually  very  large  as  compared  to  the  coating.  Hence  a 
h<i 

thickness  ratio  of  —  =  10.0  is  considered.  Figures  5.55  —  5.56  show  the  thermal  stresses 
h\ 

under  uniform  temperature  change  and  under  conduction  respectively.  Under  uniform 
temperature  rise,  the  coating  is  under  very  high  tensile  stress  whereas  the  thick  substrate  is 
under  somewhat  uniform  stress.  The  stress  intensity  factors  under  this  tensile  stress  are 
positive  (Figure  5.74)  but  as  the  crack  approaches  the  interface,  the  loading  changes  sign 
and  the  stress  intensity  factor  drops.  This  crack  is  arrested  by  the  compressive  zone 
(Figure  5.75).  For  uniform  cooling  the  problem  becomes  one  of  crack-contact  type.  The 
crack-tip  locations  and  the  stress  intensity  factors  are  shown  in  Figures  5.76  and  5.77 
respectively.  For  surface  heating  we  observe  the  arrest  of  an  open  crack  as  it  enters 
compressive  zone  (Figure  5.78). 

Some  general  comments  can  be  made  from  these  results.  It  is  observed  that 
thermal  loading  either  tends  to  arrest  the  crack  growth  or  create  a  crack-contact.  If  the 
crack  lies  completely  in  the  compressive  zone,  the  crack  will  remain  closed.  This  means 
that  every  particular  loading  is  associated  with  a  critical  length.  If  the  surface  flaw  is 
smaller  than  this  critical  length  the  flaw  will  not  propagate.  If  the  surface  flaw  does 
propagate,  it  will  be  arrested  by  the  compressive  zone  under  pure  thermal  loading.  If  other 
kinds  of  loadings  are  superposed  along  with  the  thermal  loading,  the  stress  intensity 
factors  of  the  different  loadings  can  be  added  up  to  obtain  the  resultant  stress  intensity 
factor.  Care  should  however  be  taken  to  determine  if  the  crack  is  open  or  closed  under  the 
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combination  of  loadings.  If  the  resultant  stress  intensity  factor  is  positive,  the  crack  is  open 
and  the  analysis  is  straightforward.  If,  on  the  other  hand,  the  loadings  tend  to  close  the 
crack,  the  procedure  described  in  Section  5.2.2  should  be  followed  to  find  the  point  of 
crack  closure  and  the  correct  stress  intensity  factors. 
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Present  work 


1.0172 


1.0820 


1.1546 


Ref  [43] 


1.0172 


1.0820 


1.1546 


CTQy/c 


Present  work  Ref  [43] 


1.0087 


1.0377 


1.0654 


1.1042 


1.1669 


1.0087 


1.0377 


1.0653 


2.3283 


3.248 


1.8785 


2.3105 


3.103 


1.2641 


1.3521 


Table  5.1  Comparison  of  normalized  stress  intensity  factors  for  internal  crack  in 
homogeneous  half  plane  bonded  to  nonhomogeneous  half  plane  subjected 


to  fixed  grip  loading.  Published  results  are  from  Erdogan  et.  al.[43]. 

ai  =  0.0,  c  =  (bl~Ql\  h1  =  h2  =  1000.0,  cto  = 

4  (1  +  K) 
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Fixed  grip 

Membrane  loading 

Pure  bending 

c 

£1(02) 

h(h) 

^1(02) 

h  (62) 

fcl(02) 

h  O2) 

hi 

^0  y/c 

cro^/c 

°t\Tc 

*ty/c 

°b\Tc 

vb\Tc 

0.0001 

1.0000 

1.0000 

0.5001 

0.5000 

0.7606 

0.7608 

0.05 

1.0043 

1.0047 

0.5048 

0.4998 

0.7208 

0.8074 

0.1 

1.0167 

1.0199 

0.5136 

0.5049 

0.6871 

0.8623 

0.15 

1.0370 

1.0483 

0.5263 

0.5165 

0.6596 

0.9274 

0.2 

1.0656 

1.0938 

0.5431 

0.5366 

0.6387 

1.0065 

0.25 

1.1037 

1.1633 

0.5646 

0.5686 

0.6256 

1.1059 

0.3 

1.1538 

1.2692 

0.5921 

0.6185 

0.6227 

1.2376 

0.35 

1.2213 

1.4365 

0.6282 

0.6985 

0.6349 

1.4262 

0.4 

1.3184 

1.7284 

0.6788 

0.8393 

0.6736 

1.7344 

0.45 

1.483 

2.383 

0.763 

1.157 

0.775 

2.402 

Table  5.2  Normalized  stress  intensity  factors  for  an  internal  crack  located  in  the 

(62  —  02)  (&2  +  02)  _  hi  h,2  n  n 
homogeneous  layer,  c  = - - - , - - - =  "2"’  ~~  i‘u> 

8uo  N  M 

Bh\  =  0.375 (PSZ-Rene41), op  =  — — eo,  at  =  ab  =  To- 

1  +  K  fl  1  ft} 
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Fixed  grip 

Membrane  loading 

Pure  bending 

c 

fcl(02) 

hi  (62) 

hi (02) 

h(h) 

h  (02) 

k,  (h) 

&0y/c 

croy/c 

at\Tc 

<ytyfc 

<?by/c 

crb\/c 

0.0001 

0.5000 

0.7414 

0.7416 

0.05 

1.0042 

1.0046 

0.5006 

0.7106 

0.7790 

0.1 

1.0163 

1.0194 

0.5051 

0.5127 

0.8241 

0.15 

1.0360 

0.8790 

0.2 

1.0636 

1.0915 

0.5258 

0.5518 

0.6531 

0.9470 

0.25 

1.1001 

1.1593 

0.5426 

0.5873 

0.6471 

1.0338 

0.3 

1.1480 

1.2621 

0.5651 

0.6405 

0.6504 

1.1502 

0.5958 

1.3181 

0.4 

1.3036 

1.7052 

0.6403 

0.8671 

0.7081 

1.5933 

0.45 

1.4573 

2.3304 

0.7164 

1.1855 

2.1866 

Table  5.3  Normalized  stress  intensity  factors  for  an  internal  crack  located  in  the 


homogeneous  layer,  c  = 


(62-02)  (&2  +  02) 


h  ^  =  10 

2  ’  hi 


(3h\  =  -  0.2576  (PSZ-Ti-6Al-4V),  a0 


8//0  N  M 
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Fixed  grip 

Membrane  loading 

Pure  bending 

c 

h(h) 

ki(ai) 

h(h) 

fcl(oi) 

h(h) 

fcl(oi) 

hi 

0Oy/c 

<?0\/c 

°t\Tc 

°t\fc 

ob\Tc 

<?b  \Tc 

0.0001 

0.8290 

0.8290 

_ 

0.4997 

0.4997 

0.05 

0.8275 

0.7768 

0.1 

0.8202 

0.8621 

0.15 

0.8302 

0.8893 

0.9570 

0.7369 

Itggl 

0.5712 

0.3 

0.9579 

lltVJfeM 

0.5890 

0.6011 

1.1223 

0.7440 

1.064 

1.090 

0.6399 

1.262 

0.765 

0.4 

1.251 

1.182 

■ 

0.693 

1.49 

0.81 

0.45 

1.66 

1.32 

1.03 

0.78 

1.98 

0.90 

Table  5.4  Normalized  stress  intensity  factors  for  an  internal  crack  located  symmetrically 

•  i  i  i  (6i  —  cb\ )  (b\  +  d\ )  h\  ho 

m  the  nonhomogeneous  layer,  c  =  - — - ,  - — - — -  =  — ,  —  =  1.0, 

2  2  2  h\ 

Phi  =  0.375 (PSZ-Rene41),  a0  =  crb  =  . 

1  +  k  h\  h\ 
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Fixed  grip 

Membrane  loading 

Pure  bending 

fciOi) 

ki(a,i) 

ki(oi) 

hih)  fei(ai) 

&0  y/c 

o'O  y/c 

ot\Tc 

1.1375 

1.1374 

0.5001 

0.5001 

0.7063  0.7061 

1.1546 

1.1334 

0.5059 

0.5000 

0.7562  0.6655 

1.1840 

1.1380 

0.5170 

0.5037 

0.8148  0.6314 

1.2288 

1.1510 

0.5349 

0.5110 

0.8846  0.6038 

1.2944 

1.1726 

0.5618 

0.5220 

0.9693  0.5829 

1.3894 

1.2038 

0.6014 

mmol 

1.0757  0.5698 

1.5295 

1.2473 

0.6607 

1.2165  0.5665 

1.747 

1.309 

0.7532 

0.5858 

1.418  0.578 

2.12 

1.40 

0.914 

0.627 

1.746  0.614 

2.96 

1.57 

1.28 

0.70 

2.46  0.71 

Table  5.5  Normalized  stress  intensity  factors  for  an  internal  crack  located  symmetrically 

.  ^  ,  Oi  -  ai)  Oi+ai)  hi  in 

m  the  nonhomogeneous  layer,  c  = - - - ,  - - - =  — ,  —  =  1.0, 

2  2  2  hi 

8un  N  M 

Phx  =  -  0.2576 (PSZ-T1-6A1-4V),  a0  =  at  =  —,ab  = 

1  +  k  hi  h\ 
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fci(oi) 

y/c 

Case  II 

Caselll 

0.01 

0.8290 

0.8333 

0.05 

0.8420 

0.8448 

0.8721 

0.10 

0.8621 

0.8656 

0.9004 

0.15 

0.8893 

0.8936 

umi 

0.20 

0.9239 

0.9291 

0.25 

0.9672 

0.9735 

1.0393 

0.30 

1.0212 

1.0287 

1.1093 

0.35 

1.090 

1.0991 

1.1985 

0.40 

1.182 

1.194 

1.318 

0.45 

1.32 

1.34 

1.50 

Table  5.6  Normalized  stress  intensity  factors  for  an  internal  crack  located  symmetrically 

(fej  —  CL  i) 

in  the  nonhomogeneous  layer  subjected  to  fixed  grip  loading,  c  = - - - , 

z 

_  t±t  h  -  i.o,  fihi  =  0.375  (PSZ-Rene41),  a0  = 

2  2  h\  1  +  k 

Case  I :  Crack  only  in  nonhomogeneous  layer 

Case  II :  Crack  in  presence  of  another  crack  in  homogeneous  substrate, 


<i2  =  0Ah\,  62  =  O.6/11 


Case  ffl :  Crack  in  presence  of  another  crack  in  homogeneous  substrate, 


02  —  0.2fei,  62  =  0.8/&1 
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Fixed  grip 

Membrane  loading 

Pure  bending 

c 

Hh) 

fcl(ai) 

ki(bi) 

fci(ai) 

fci(ai) 

hi 

<70 

at\Tc 

crty/c 

<?by/c 

ab\Tc 

0.8286 

0.8333 

Itflgl 

0.8218 

0.8448 

0.4973 

0.5067 

0.10 

0.8212 

0.8656 

0.4990 

0.5169 

0.7556 

0.15 

0.8307 

0.5313 

0.7415 

0.20 

0.8530 

0.9291 

0.5221 

0.5502 

0.9571 

0.7331 

0.25 

0.8925 

0.9735 

0.5743 

1.0273 

liWtiEl 

0.30 

0.9571 

1.0287 

0.6049 

1.1229 

0.35 

1.0629 

1.0991 

0.6547 

0.6445 

1.2624 

0.7592 

0.40 

1.250 

1.194 

0.770 

0.698 

0.802 

0.45 

1.66 

1.34 

1.02 

0.78 

0.89 

Table  5.7  The  normalized  stress  intensity  factors  for  the  internal  crack  located 

symmetrically  in  the  nonhomogeneous  material  in  presence  of  a  crack  in  the 

1.0,  CL2  =  0.4/&1, 


t  ,  Oi  -  Oi)  (h  +  Oi)  hi  h2 

homogeneous  matenal.  c  = -  - —  —  — 


2  ’  2 
i>2  =  0.6hi,  /3h\  =  0.375  (PSZ-Rene41),  <7q  = 


8^o 

1  +  K 


2  ’  hi 


€0,  at  = 


N 

i  >  <7& 
hi 


M 

hr 
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Fixed  grip 

Membrane  loading 

Pure  bending 

c 

ki(bi) 

h(h) 

fc,(h) 

ki(b2) 

ki(bi) 

h(h) 

hi 

% 

O 

b 

% 

O 

b 

°t\Tc 

0.0001 

0.9998 

1.0001 

0.5513 

-0.1035 

-0.1018 

0.1 

0.9836 

1.0115 

0.5473 

0.5526 

-0.1843 

-0.0174 

mm 

0.9799 

1.0351 

0.5499 

-0.2615 

0.0672 

mm 

0.9903 

1.0724 

0.5599 

-0.3352 

0.1514 

n 

1.0178 

1.1272 

0.5792 

-0.4074 

0.2366 

mm 

1.0682 

1.2060 

0.6111 

WiM££A 

-0.4816 

0.3258 

0.6 

1.1528 

1.3222 

0.6622 

0.6957 

-0.5647 

0.4253 

1.2958 

1.5034 

0.7465 

-0.6708 

0.5488 

1.5606 

1.8216 

0.9007 

-0.8364 

0.7333 

2.1983 

2.5595 

1.2705 

-1.2073 

1.1269 

3.1251 

3.6099 

1.8077 

1.8686 

-1.7429 

1.6789 

Table  5.8  The  normalized  stress  intensity  factors  for  a  crack  symmetrically  located  with 


.  h<i 

respect  to  the  interface.  — 

hi 


<70 


8^0 


N 


If.  (frl  +  h)  Rh 
1.0,  c  ^  >  ph l 

M 


0.375  (PSZ-Rene41), 


l  +  «€°’  at~  hi,(7b~  h\ 
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h 

croy/b 

b/h 

Present 

Ref  63 

0.001 

1.12152 

1.121531 

0.05 

1.1399 

0.1 

1.1892 

1.1892 

1.2652 

iw 

1.3673 

1.3673 

0.25 

1.4975 

0.3 

1.6599 

0.4 

2.1115 

2.1114 

0.45 

2.4254 

0.5 

2.8247 

2.8246 

0.55 

3.3430 

0.6 

4.0333 

4.0332 

0.65 

4.9846 

0.7 

6.3552 

6.3549 

0.75 

8.4536 

0.8 

11.9559 

11.955 

0.85 

18.6285 

18.628 

0.9 

34.6336 

34.632 

0.95 

99.121 

99.128 

0.98 

394.3 

Table  5.9  Comparison  of  results  for  stress  intensity  factors  for  edge  crack  in 

homogeneous  strip  under  fixed  grip  loading.  In  solution  in  this  work, 

B  h  8uo 

nonhomogeneity  parameter  —  =  0.0001,  &o  =  .  —  .€o. 

2  i  K>) 
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b/h 

h 

ah\Tb 

Present 

Ref  63 

0.001 

1.12161 

1.1202 

0.05 

0.1 

1.0472 

0.15 

1.0433 

1.0553 

FUMIM 

1ECT1 

1.1241 

him 

Ilia 

0.4 

1.2609 

1.2606 

0.45 

0.5 

FIEEW 

1.4972 

0.55 

1.6749 

0.6 

1.9140 

0.65 

Em 

2.7252 

3.4588 

0.8 

4.6772 

0.85 

IdkkVIilM 

0.9 

1 

0.95 

34.307 

34.306 

Table  5.10  Comparison  of  results  for  stress  intensity  factors  for  edge  crack  in 

homogeneous  strip  under  pure  bending.  In  solution  in  this  work,  nonhomogeneity 

(3h  n  _  6M 

parameter  —  =  0.0001,  =  —5-. 

2  hz 
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c 


hi 


0.0001 


Fixed  grip 

Membrane 

ki(ai) 

fci(oi) 

<70  y/c 

otyfc 

1.1215 

0.5042 

1.2142 

1.3091 

1.4356 

0.6214 

1.5974 

0.6855 

1.8016 

2.0591 

2.3868 

1.0010 

3.3871 

1.4027 

Bending 


fci(ai) 


ObJc 


1.9805 


2.2418 


Table  5.11  The  normalized  stress  intensity  factors  for  an  edge  crack  located  in  the 

h>2 

nonhomogeneous  layer,  c  =  (b\  —  ai),  6i  —  hi,  —  =  1.0, 

h\ 


(5  hi  =  0.375  (PSZ-Rene41),  cto 


8 /ii  N  M 

1+/C60’  at~  h,'ah~  h\- 
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Fixed  grip 

Membrane 

Bending 

c 

h(b2) 

ki(b2) 

h(b2) 

hi 

°0y/c 

°t\Tc 

°b\fc 

1.1 

4.0518 

1.6674 

1.2 

4.9205 

2.0146 

1.3 

PgUlM 

1.4 

7.8292 

1.5 

10.4546 

4.227 

1.6 

14.8381 

5.978 

1.7 

23.191 

9.314 

1.8 

43.241 

17.317 

1.9 

124.063 

49.554 

51.444 

Table  5.12  The  normalized  stress  intensity  factors  for  an  edge  crack  crossing  the 

h2 

interface,  c  —  (b\  +  62),  i>i  =  hi,  —  =  1.0, 

h\ 

/3hi  =  0.375  (PSZ-Rene41),  a0  =  -8;--  €0,  ot  =  t~,  &b  =  To- 

1  +  n  h\  h\ 
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Fixed  grip 

Membrane 

Bending 

c 

fci(ai) 

fci(ai) 

ki(ai) 

h 

0.0001 

1.1215 

0.5942 

1.8892 

0.1 

1.1337 

0.6060 

1.7897 

0.2 

1.1744 

1.7346 

0.3 

1.2392 

1.7117 

0.7261 

1.7146 

1.4384 

0.7924 

1.7404 

0.6 

1.5765 

0.8736 

1.7886 

0.7 

1.7456 

1.8600 

0.8 

1.9522 

1.9574 

0.9 

2.2044 

1.2401 

2.0842 

1.0 

2.5082 

1.4168 

2.2413 

Table  5.13  The  normalized  stress  intensity  factors  for  an  edge  crack  located  in  the 

]%2 

nonhomogeneous  layer.  c  =  (&i  —  ai),  61  =  hi,  —  =  1.0, 

.  8/ii  N  M 

/3hi  =  -  0.2576  (PSZ-Ti-6Al-4V),  a0  =  -^-e0,  at  =  —,ab=  -5. 

1  +  k  h\  h\ 
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Fixed  grip 

Membrane 

Bending 

c 

ki(b2) 

ki(b2) 

mm 

°t\fc 

KDI 

2.5166 

1.2 

3.5489 

2.0178 

2.8731 

1.4 

5.5566 

mmmi 

1.5 

7.3730 

4.227 

1.6 

10.4032 

5.977 

BUM 

9.314 

1.8 

30.018 

mm 

Table  5.14  The  normalized  stress  intensity  factors  for  an  edge  crack  crossing  the 

interface,  c  =  (&i  +  b2),  h  =  hh  ^  =  1.0,  f3hx  =  -  0.2576  (PSZ-Ti-6Al-4V), 

hi 

8/xi  N  M 


Fixed  grip 

Membrane 

Bending 

c 

fci(ai) 

fci(ai) 

fci(ai) 

hi 

ob\Tc 

0.0001 

1.1215 

0.07434 

0.04040 

0.1 

1.1344 

0.07516 

0.04042 

0.2 

1.1551 

0.07649 

0.04069 

0.3 

1.1816 

0.07821 

0.04114 

0.4 

1.2132 

0.08026 

0.04175 

0.5 

0.08263 

0.04249 

0.6 

0.08530 

0.04336 

0.7 

1.3365 

0.08829 

0.04436 

0.8 

1.3879 

0.09165 

0.04550 

0.9 

1.4462 

0.09545 

0.04682 

1.0 

1.5170 

0.1001 

0.04850 

Table  5.15  The  normalized  stress  intensity  factors  for  an  edge  crack  located  in  the 

/&2 

nonhomogeneous  layer,  c  =  (&i  —  ai),  b\  =  hi,  —  —  10.0, 

hi 

Phi  =  0.375  (PSZ-Rene41),  a0  =  .^~-  €i,  at  =  <?b  = 

1  +  K  n\ 
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Fixed  grip 

Membrane 

Bending 

c 

ki(b2) 

ki(b2) 

h(h) 

hi 

V0  y/c 

Vty/c 

Vbyfc. 

2.0 

1.8149 

0.1192 

0.05122 

3.0 

2.1694 

0.05420 

4.0 

2.6796 

0.1747 

0.05946 

5.0 

wmmmm 

0.06807 

6.0 

0.08211 

7.0 

6.5909 

0.4265 

8.0 

10.3839 

9.0 

19.5218 

1.2593 

0.2645 

10.0 

56.381 

3.6327 

0.7078 

Table  5.16  The  normalized  stress  intensity  factors  for  an  edge  crack  crossing  the 


interface,  c  =  (£>i  +  62),  &i  =  hi,  —  =  10.0,  /3h\ 

hi 

8/xi  N  M 

vq  1  .  eo,  Vt  ,  >  —  ,0  • 

1+ «  hi  h{ 


=  0.375  (PSZ-Rene41), 
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c 

ki 

hi 

aTy/c 

0.0 

0.0708 

0.1 

0.0673 

0.2 

0.0649 

0.3 

0.0628 

0.4 

0.0607 

0.5 

0.0580 

0.6 

0.0545 

0.7 

0.0497 

0.8 

0.0430 

0.9 

0.0336 

1.0 

0.0203 

1.1 

0.0097 

1.2 

0.0049 

1.3 

0.0026 

1.4 

0.0008 

1.5 

-  0.0003 

1.6 

-  0.0005 

1.7 

-  0.0009 

1.8 

-  0.0001 

Table  5.17  Normalized  stress  intensity  factors  for  an  edge  crack  passing  through 
the  interface  under  uniform  temperature  rise. 


/l2 

c  =  crack  length,  b\  =  h\,—  =  1.0,  a?  = 

h\ 


8/^qq;o(T  —  Tp)(l  +  v) 
(1  +  k) 


Vo  . 
Vi 

Ph 


>Ph  ,  —  =  eahl ,  ^  =  e^\Ti  =  T2  =  T  >  T0, 
ar  ki 

0.375,  a/ii  =  0.513,  7/11  =  2.5  (PSZ-Rene41). 
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X 

hi 


Figure  5.2  Stress  distribution  for  an  uncracked  specimen  under  membrane 

i  h2  N 

loading.  —  =  1.0,  at  =  — . 

fti  ri\ 
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hi 


Figure  5.3  Stress  distribution  for  an  uncracked  specimen  under  pure  bending. 
h2  M 

h  =  10>  04 = W 
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factors  for  an  internal  crack  symmetrically  located  in  the  homogeneous  layer 


subjected  to  fixed  grip  loading. 
(p2  —  Q%)  (&2  +  U2)  ^1 


-  —  —  =  1.0,  CT0  = 


8M0 

(1  +  k) 


e0- 


29 


c 

~h\ 


Figure  5.7  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 

factors  for  an  internal  crack  symmetrically  located  in  the  homogeneous  layer 

subjected  to  membrane  loading. 

(b 2  —  a<i)  (p2  +  02)  hi  h%  N 

C  =  2  ’  2  =  ~2’1h=  1'0,  ^  =  ft? 
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Figure  5.8  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 

factors  for  an  internal  crack  symmetrically  located  in  the  homogeneous  layer 

subjected  to  pure  bending. 

(62-02)  (62+02)  h\  h 2  in 


M 


Figure  5.9  Effect 


displacement  for  an  internal  crac] 
grip  loading,  v*  =  2h1e0,  ^  = 


ameter  on  normalized  crack  opening 


in  the  homogeneous  material  subjected  to  fixed 


Figure  5.10  Effect  of  nonhomogeneity  parameter  on  normalized  crack  opening 

displacement  for  an  internal  crack  in  the  homogeneous  material  subjected  to 
membrane  loading,  v*  —  2 N  ^  =  1.0,  02  =  0.25hi,  i>2  =  0.75hi. 

hi 


Bmf] 


Figure  5.11  Effect  of  nonhomogeneity  parameter  on  normalized  crack  opening 

displacement  for  an  internal  crack  in  the  homogeneous  material  subjected 

M  (\  K/')  ho 

to  pure  bending,  v*  —  2  —  —  =  1.0,  a<i  =  0.25/ii,  &2  =  0.75/&1. 

h\  8/io  ^1 


hi 


Figure  5.13  The  effect  of  thickness  ratio  on  normalized  stress  intensity  factors  for  an 

internal  crack  located  in  the  homogeneous  layer  subjected  to  membrane  loading. 
_  (p2  -  Q2)  (62  +  02)  _  hi  Rh  _  ,  fincu  _ 

c  —  2  )  2  ~  2  ’  ~  <Jt  —  ^  . 
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Figure  5.14  The  effect  of  thickness  ratio  on  normalized  stress  intensity  factors  for  an 


internal  crack  located  in  the  homogeneous  layer  subjected  to  pure  bending. 


-  1-6094,  o»  =  §. 


c  — 


hi 


Figure  5.15  Effect  of  thickness  ratio  on  normalized  crack  opening  displacement  for  an 


internal  crack  in  the  homogeneous  material  subjected  to  fixed  grip  loading. 
=  2hieo,  0hi  =  1.6094,  a2  =  Q.25hh  h  =  0.75hh 


Figure  5.16  Effect  of  thickness  ratio  on  normalized  crack  opening  displacement  for  an 
internal  crack  in  the  homogeneous  material  subjected  to  membrane  loading. 

V *  =  2iV  ^  — ,  0hi  =  1.6094,  02  =  0.25hi,  b2  =  0.75hv 


Figure  5.17  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 
factors  for  an  internal  crack  symmetrically  located  in  the  nonhomogeneous 


layer  subjected  to  fixed  grip  loading,  c  = 
h<i  _  _  8^q 

l.Oj  (7q  .  \  €0* 

hi  (1  +  k) 


(bi  —  ai)  (&i+ai)  hi 
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Figure  5.19  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 
factors  for  an  internal  crack  symmetrically  located  in  the  nonhomogeneous 

i  i  •  t  ,  ( b\  —  ai)  (b\  -f~  oi)  h\ 

layer  subjected  to  pure  bending,  c  =  - — - — - — - — -  =  — , 


////////////////// 


c 

~hi 


Figure  5.21  Normalized  stress  intensity  factors  for  an  internal  crack  located  in  the 


nonhomogeneous  layer  touching  the  interface  and  subjected  to  fixed  grip  loading. 

8^0 


C  =  —  9-ai\  01  =  0.0,  ^  =  0.2,  a0  = 

2  h  i  (1  +  k) 


EQ,(3h\=  —0.2576. 
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Figure  5.22  Effect  of  nonhomogeneity  parameter  on  normalized  crack  opening 

displacement  for  an  internal  crack  in  the  nonhomogeneous  material  subjected  to 

ho 

fixed  grip  loading,  v*  =  2hieo,  —  =  1.0,  b\  =  0.75Ai,  a\  =  0.25hi. 

hi 
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Figure  5.23  Effect  of  nonhomogeneity  parameter  on  normalized  crack  opening 

displacement  for  an  internal  crack  in  the  nonhomogeneous  material  subjected  to 

membrane  loading,  v*  =  2 N  ^  ^  =  1.0,  2>i  =  0.75/z-i ,  a\  =  0.25hi. 

8/io  m 


Figure  5.24  Effect  of  nonhomogeneity  parameter  on  normalized  crack  opening 

displacement  for  an  internal  crack  in  the  nonhomogeneous  material  subjected  to 
pure  bending,  v*  =  2^-  ^  +  — ,  ^  =  1.0,  bi  =  0.75/ii,  ai  =  0.25/ii. 

hi  8uq  hi 
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Figure  5.26  Effect  of  thickness  ratio  on  normalized  crack  opening  displacement  for  an 
internal  crack  in  the  nonhomogeneous  material  subjected  to  fixed  grip  loading. 
v *  =  2frie0,  &h\  =  1.6094,  i>i  =  0.75/*i,  ai  =  0.25/i2. 
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Figure  5.29  The  effect  of  thickness  ratio  on  normalized  stress  intensity  factors  for  an 

internal  crack  symmetrically  located  in  the  nonhomogeneous  layer  subjected  to 

r-  j  •  i  j-  (Pi  —  a\)  (&i+ai)  hi  o*7c 

fixed  grip  loading,  c  = - - - , - - - =  —  ,ph\  =  0.375, 

ao  =  &2  =  OAhi ,  62  =  0.6hj. 


internal  crack  in  the  nonhomogeneous  material  subjected  to  pure  bending. 
M  (1  +  k) 

v*  =  21Tc - '*  =  0.375,  bi  =  0.75/ii,  Gl  =  0.25hh  a2  =  0.4hh 

hi  o//o 


b2  =  0.6/ii. 


Figure  5.33  The  effect  of  thickness  ratio  on  normalized  stress  intensity  factors  for  a  crack 
symmetrically  located  with  respect  to  the  interface  and  subjected  to  pure  bending. 
Phi  =  0.375,  ab=^. 


Figure  5.34  Effect  of  nonhomogeneity  parameter  on  normalized  crack  opening 

displacement  for  a  crack  crossing  the  interface  and  subjected  to  membrane 
loading,  v*  =  2 N ^o+  — ,  ^  =  1.0,  &i  =  0.2hi,  &2  =  0.2/zi. 


////////////////// 


c 

hi 


Figure  5.35  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 

factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer  subjected  to  fixed 

Su  i 

grip  loading,  c  =  (&i  -  oi),  h  =  hh  —  =  1.0,  cr0  =  -  —  -reo- 

n\  (1  +  Kj 
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Figure  5.36  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 


factor  for  an  edge  crack  crossing  the  interface  subjected  to  fixed  grip  loading. 


c  =  crack  length,  £>i  =  hi,  —  =  1.0,  a0  = 


8/zi 


(1  +  k) 


€0- 


hi 


Figure  5.37  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 

factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer  subjected  to 

membrane  loading.  c=  (b\-  ai),  b\  —  hi,—  =  1.0,  at  =  — . 

hi  hi 
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Figure  5.38  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 

factor  for  an  edge  crack  crossing  the  interface  subjected  to  membrane  loading. 

ho  N 

c  —  crack  length,  i>i  =  h\,  —  =  1.0,  at  =  — . 


Figure  5.39  The  effect  of  nonhomogeneity  parameter  on  normalized  stress  intensity 

factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer  subjected  to  pure 

.  ,,  ..  .  _  _  h%  ilf 

bending,  c  =  (&i  -  01),  t>i  =  fti,  —  =  1.0,  ab  = 


IKII 


Figure  5.41  The  effect  of  thickness  ratio  on  normalized  stress  intensity  factor  for  an  edge 
crack  located  in  the  nonhomogeneous  layer  subjected  to  fixed  grip  loading, 
c  =  (pi  —  &i),  bi  =  hi,  (3 h\  —  0.375,  <Jq  =  ~  :  r^o- 
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Figure  5.43  The  effect  of  thickness  ratio  on  normalized  stress  intensity  factor  for  an  edge 

crack  located  in  the  nonhomogeneous  layer  subjected  to  pure  bending. 

M 

c  =  (bi-  01),  £>i  =  -  0.375,  ah  = 


////////////////// 


muWwww 


x 

hi 

Figure  5.45  Effect  of  nonhomogeneity  parameter  on  normalized  crack  opening 

displacement  for  an  edge  crack  in  the  nonhomogeneous  material  subjected  to 
fixed  grip  loading,  v*  =  2/iieo,  7^  =  1.0,  =  hi,  ai  =  0.6hi. 

hi 
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h~i 

Figure  5.47  Effect  of  nonhomogeneity  parameter  on  normalized  crack  opening 

displacement  for  an  edge  crack  in  the  nonhomogeneous  material  subjected  to 

membrane  loading,  v*  =  2  N  ^  =  1.0,  bi  =  hi,  ai  =  0.6hi. 

8/io  hi 
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Figure  5.55  Thermal  stress  distribution  for  an  uncracked  specimen  under  uniform 

temperature  change.  =  8*>°»Ar(l  +  ,)  W  =  ^ 

(1  +  k)  Ml 


—  =  eahl ,  ^  =  e7/il ,  phi  =  0.375,  ahi  =  0.513,  7/11  =  2.5, 


,7*1 
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Figure  5.58  The  effect  of  uniform  temperature  rise  on  normalized  stress  intensity 

factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer. 

_ fh  n  \  h  h  in-r  Shqcxo(T  —  Tq)(1  +  v) 

C  =  (01  -  Oi),  b  1  -  hi,  —  -  1.0,  aT  = - 7— — r - , 

hi  (1  +  k) 

—  =  e^1 ,  —  =  eah\y-  =  e7Al,Ti  =T2=T>T0, 

Ml  <*i  kc 

Phi  =  0.375,  ahi  =  0.513,  7^1  =  2.5  (PSZ-Rene41). 
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Figure  5.61  Crack  closure  length  as  a  function  of  the  location  of  the  open  crack  tip 


for  an  edge  crack  passing  through  the  interface  loaded  thermally. 

T2  =  T0,Ti  >  T0,  d  =  kl~b\  ^  =  1.0,  ^  ^  =  eah\ 

fl\  h\  jJL\  OL\ 

kf) 

=  e7hl ,  f3h\  =  0.375,  ahi  =  0.513,  7/11  =  2.5  (PSZ-Rene41). 


c\  =  hi  —  ai  for  crack  in  coating,  c\  =  h\+  b%  for  through  crack 
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Figure  5.62  The  effect  of  temperature  rise  at  one  boundary  on  normalized  stress 


intensity  factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer. 

(0i  —  Ci)  ho 

c  — - - — b\  determined  from  Fig.  5.60,  —  =  1.0, 

2  h\ 

8(xoao(Ti  -  T0)(l  +  u)  Mo  0h,  a 0  0 

crT  = - .  -  r - ,  T2  =  T0,  Ti  >  T0,  —  =  ep  1 ,  —  =  e 

(1  + «)  Mi  ai 

T2  =  e7/l1 ,  f3h\  =  0.375,  ah\  =  0.513,  7 =  2.5  (PSZ-Rene41). 


a 
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Figure  5.63  The  effect  of  temperature  rise  at  one  boundary  on  normalized  stress 


intensity  factor  for  an  edge  crack  passing  through  the  interface. 

ho 

c  =  half  crack  length,  b\  determined  from  Fig.  5.61,  —  =  1.0, 

hi 

8//QQ0 (T\  -  To)(l  +v)  n  _m_  rr.  ^  m.  ^  Oq  _  h 


ot  — 

^2  _  cTlh 

be 


72  =  To, Ti  >  To,  —  = 
(1  +  k)  /xi 


ai 


Y  =  e7  1 ,  (3h\  =  0.375,  ahi  =  0.513,  7hi  =  2.5  (PSZ-Rene41). 


ci  =  hi  —  ai  for  crack  in  coating,  ci  =  hi  +  62  for  through  crack 
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Figure  5.65  The  effect  of  uniform  temperature  fall  on  normalized  stress  intensity 


factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer. 

(b i  —  Gi )  ho 

c  = - - — -,bi  determined  from  Fig.  5.64,  —  =  1.0, 

2  h\ 

=  8/XQQo(To-r)(l-h^)  fi o  _  nhl  ao  _  h.  ko  _  h 
T  (1  +  k)  ’m  ’ar"6  ’ kc~e 

Ti  =  T2  =  T  <  To,/3hi  =  0.375,  a/ii  =  0.513,  7/11  =  2.5. 
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Figure  5.66  The  effect  of  temperature  fall  at  one  boundary  on  normalized  stress 


intensity  factor  for  an  edge  crack  passing  through  the  interface. 


i  i  l  i.  h  2  in  8^oQo(^b  ~T\)(l  +  u) 

c  =  crack  length,  bi  -  hi,  —  =  1.0,  aT  = - 7— — r - , 

hi  (1  +  k) 


T2  =  T0,  Ti  <  T0,  ^  =  e0hi ,  —  =  eG\  ^  =  e<h, 

Hi  ai  kc 


phi  =  0.375,  cthi  =  0.513,  7/11  =  2.5  (PSZ-Rene41). 


Figure  5.67  Crack  closure  length  as  a  function  of  the  location  of  the  open  crack  tip 


for  an  edge  crack  under  uniform  temperature  rise.  T\  =  T2  =  T  >  Tq, 


d  =  hl7b\  h  =  10,  ^  ^  =  eahl  ^ 


hi 


hi  hi 


Ql 


'  kr 


=  e' 


(3hi  =  -  0.2576,  ahi  =  -  0.0513,  7/11  =  1.2777  (PSZ-Ti-6Al-4V). 


Figure  5.68  The  effect  of  uniform  temperature  rise  on  normalized  stress  intensity 
factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer. 
c  =  —  -  Ql^ ,  i>i  determined  from  Fig.  5.67, 

^  =  1.0,T1  =  T2  =  T>T0, 

hi 

_  8/ioa0(r  -  Tp)(l  +  v)  vo  _  nhl  ao  _  afh  h  _ 

-  /1.\  *  -  C 

(l  +  «)  Ml  OL\  kc 

fihi  =  -  0.2576,  ahi  =  -  0.0513,  7 hi  =  1.2777  (PSZ-Ti-6Al-4V). 


cjtVc  u-u^ 


c 


hi 


Figure  5.69  The  effect  of  uniform  temperature  fall  on  normalized  stress  intensity 

factor  for  an  edge  crack  passing  through  the  interface. 

ii  -■»  t  ,  h%  ~  n  SfiQao (To  —  T)(l  +  u) 

c  =  crack  length,  bx  =  hh  —  =  1.0,  aT  =  —  . - ^ - -, 

(1  +  k) 

l ^2  _  pPh  _  cthi  ko  _ 

CL  1  fcc 


e"'"1  ,  —  =  e“'n ,  e  =  e7/l1 ,  Ti  =  T2  =  T  <  T0, 


/3h!  =  -  0.2576,  afti  =  -  0.0513,  7/^  =  1.2777  (PSZ-Ti-6Al-4V). 


hi  —  Oj\ 

^1 


Figure  5.70  Crack  closure  length  as  a  function  of  the  location  of  the  open  crack  tip 

for  an  edge  crack  loaded  thermally.  T2  =  To,  T\  >  To, 
hi-bi  fc2  m  =  e»  «o  =  <*  to  =  * 

7  5  7  ■L,w)  J  ? 


/Shi  =  -  0.2576,  a/»i  =  -  0.0513,  7/4  =  1.2777  (PSZ-Ti-6Al-4V). 
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Figure  5.71  The  effect  of  temperature  rise  at  one  boundary  on  normalized  stress 

intensity  factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer. 

(oi  —  0,] )  ho 

C  - - — - ,  b\  determined  from  Fig.  5.70,  —  =  1.0 ,T2  =  To,Ti  >  T0 

l  hi 

aT  =  8^qQq(^i  -Tq)  (1  +  t/)  ^  _  phl  £o  _  ah  h  _  ^ 

(1  +  k)  fix  a\  kc 

(3hi  =  -  0.2576,  ahi  =  -  0.0513,  7/4  =  1.2777  (PSZ-Ti-6Al-4V). 


Figure  5.73  Comparative  effect  of  temperature  fall  at  one  boundary  on 


normalized  stress  intensity  factor  for  an  edge  crack  passing  through  the 
interface. 

11  ,  ,  h%  8fj,octo(To  —  Ti)(l  +  v) 

c  =  crack  length,  h  =  h1,  —  =  1.0,  aT  =  - ^7— — ^ - -, 

hi  (1  +  k) 

F'Q  _  31t  go  _  o/l;  fcQ  'y/lx  rp  rp  rT  ^  HT 

—  =  e  >  —  =e  :,  —  =  e7  sr2  =  T0,Ti  <  T0. 


c 

h[ 


Figure  5.74  The  effect  of  uniform  temperature  rise  on  normalized  stress  intensity 


factor  for  an  edge  crack  located  in  the  nonhomogeneous  layer. 


C  =  (bi-  o,),  h  =  hu£  =  10.0,  oT  =  8Wtto(:r _Ml+y) 

h\  (1  +  k) 


^  =  e^1 ,  —  =  eahl ,  ^  =  eyhl ,  T\  =  T2  =  T  >  T0, 
Pi  <*1  fcc 


Phi  =  0.375,  a/ti  =  0.513,  7/11  =  2.5  (PSZ-Rene41). 
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Figure  5.75  The  effect  of  uniform  temperature  rise  on  normalized  stress  intensity 

factor  for  an  edge  crack  passing  through  the  interface. 

c  =  crack  length,  i>i  =  hi,  ^  =  10.0,  aT  =  ^hoQo(T — jo](l  +  v) 

hi  (1  +  k) 

—  =  ephl ,  —  =  eQh\  ^  =  e7/l1 ,  Ti  =  T2  =  T  >  T0, 

Ml  <*1  &C 

Phi  =  0.375,  ahi  =  0.513,  yhi  =  2.5  (PSZ-Rene41). 
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Figure  5.76  Crack  closure  length  as  a  function  of  the  location  of  the  open  crack  tip 

for  an  edge  crack  passing  through  the  interface  under  uniform  temperature 

fall.Ti  =  T2=T  <T0,  d  =  h1-b1,^  =  10.0, 

h\ 

E2.  —  pph\  _  p&hi  h  _  yhi 

Vl  ctl  kc 

/3hi  =  0.375,  ahi  =  0.513,  7/11  =  2.5  (PSZ-Rene41), 


ci  =  hi  —  a\  for  crack  in  coating,  ci  =  hi  +  b2  for  through  crack. 


9* 


h\ 

Figure  5.77  The  effect  of  uniform  temperature  fall  on  normalized  stress  intensity 
factor  for  an  edge  crack  passing  through  the  interface,  c  =  half  crack 
length,  i>i  determined  from  Fig.  5.76,  —  10.0,  Tj  =  T%  =  T  <  To, 

hi 

aT  —  8^Gq0(Tq  —  T)(l  +  v)  Vo  _  phl  ao  _  ki  ko  _  hi 

T  (1  +  k)  >i  ’a!-6  'kc~e  ’ 

P hi  =  0.375,  ahi  =  0.513,  'yhi  =  2.5  (PSZ-Rene41), 

ci  =  hi  —  ai  for  crack  in  coating,  ci  =  hi  +  for  through  crack. 
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Figure  5.78  The  effect  of  temperature  rise  at  one  boundary  on  normalized  stress 


intensity  factor  for  an  edge  passing  through  the  interface. 


c  =  crack  length,  6,  =  fta,  £  =  10.0,  <rr  =  -  To)(l  ±g), 

hi  (l  +  «) 

—  =  e^1  ,  —  =  eahl ,  =  e7\T2  =  T0,Ti  >  T0 

Mi  <*i  ke 


Ph  =  0.375,  a/ii  =  0.513,  7/4  =  2.5  (PSZ-Rene41). 
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CHAPTER  6 


CONCLUSIONS  AND  FUTURE  WORK 


6.1  Conclusions 

In  the  work  presented  in  this  dissertation,  the  basic  mode  I  crack  problems  of 
Functionally  Graded  Material  coatings  were  studied.  The  substrate  was  assumed  to  be 
homogenous  and  of  finite  thickness.  This  allowed  a  parametric  study  to  determine  the 
influence  of  the  substrate  thickness  on  crack  propagation.  The  model  developed  in  this 
work  can  accommodate  internal  as  well  as  edge  cracks.  Also,  the  problem  of  crack 
crossing  the  interface  has  been  investigated. 

From  a  theoretical  point  of  view,  the  most  important  feature  of  this  study  is  the 
determination  of  the  singularity  at  the  interface  when  a  crack  either  touches  the  interface 
or  when  a  crack  passes  through  the  interface.  In  all  previous  works  published  in  literature 
on  this  subject,  it  was  assumed  that  there  is  no  singularity  at  the  interface.  The  assumption 
is  based  on  the  physical  nature  of  FGMs.  Because  the  material  property  remains 
continuous  through  the  interface,  it  was  assumed  that  there  will  not  be  any  singularity  at 
the  interface.  Finite  element  programs  based  on  this  assumption  lead  to  logical  results. 
Hence  the  assumption  has  been  widely  accepted.  But  no  mathematical  proof  of  it  was 
provided.  The  "kink"  in  the  material  property  arising  from  the  discontinuity  in  the  first 
derivative,  though,  always  raises  some  questions  as  to  whether  the  assumption  is  well- 
founded.  In  [43],  it  was  shown  that  if  a  crack  in  the  homogeneous  medium  touches  the 
interface,  there  is  no  singularity  due  to  the  interface.  This  is  only  one  half  of  the  problem, 
though.  The  real  challenge  is  if  the  crack  lies  in  the  nonhomogeneous  medium.  In  this 
study,  an  elaborate  asymptotic  analysis  was  performed  at  the  interface.  After  rigorous 
manipulations  and  after  some  very  meticulous  algebra,  it  was  shown  conclusively  that,  in 
fact,  there  is  no  singularity  at  the  interface.  It  does  not  matter  whether  the  crack  is  in  the 


202 


homogeneous  medium  or  the  nonhomogeneous  medium  or  whether  the  crack  is  crossing 
the  interface. 

From  a  practical  point  of  view,  the  most  important  question  is  the  prediction  of 
the  life  of  the  coating  under  various  loading  conditions.  Hence  a  comprehensive  fracture 
study  of  the  FGM  coating  bonded  to  the  homogeneous  substrate  was  performed.  The 
model  provided  is  very  flexible  for  the  fracture  study.  It  can  accommodate  a  range  of 
crack  configurations.  This  has  facilitated  a  systematic  parametric  fracture  study  for  internal 
cracks,  surface  cracks  and  combination  of  two  collinear  cracks.  Also,  a  range  of  loading 
conditions  were  considered.  The  applied  loads  were  chosen  so  as  to  model  the  loading 
conditions  in  practical  situations.  Both  mechanical  and  thermal  loadings  were  considered. 
The  results  of  interest  are  the  stress  intensity  factors  and  the  crack  opening  displacements. 

It  was  shown  that  the  stress  intensity  factors  tend  to  increase  as  the  crack  enters 
the  stiffer  part  of  the  graded  material.  This  effect  was  previously  observed  in  [42, 43, 46]. 
The  increase  in  the  stress  intensity  factor,  though,  should  not  be  automatically  assumed  to 
lead  to  crack  propagation.  Crack  propagation  depends  not  only  on  the  stress  intensity 
factor  (driving  force)  but  also  on  the  material  toughness  (resisting  force).  In  the  case  of 
metal-ceramic  Functionally  Graded  Materials,  the  material  toughness  is  thought  to 
increase  with  the  metal  content.  An  elaborate  metallurgical  analysis  of  the  material 
toughness  is  required  to  determine  the  exact  relation  between  the  driving  force  and  the 
resisting  force. 

Another  parameter  of  interest  was  the  substrate  thickness.  It  was  observed  that  the 

stress  intensities  were  considerably  lower  for  thicker  substrates.  This  indicates  that  the 

coating  could  be  made  to  last  longer  if  the  substrate  thickness  is  relatively  higher.  In  the 

case  of  edge  cracks,  it  was  observed  that  doubling  the  coating  thickness  reduced  the  stress 

intensity  factor  for  a  crack  reaching  the  interface  to  almost  half.  Hence  it  is  suggested  that 

h<2 

the  ratio  of  substrate  thickness  to  coating  thickness  —  should  be  maintained  at  safe  levels 
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of  4  or  higher.  Fortunately,  this  is  the  case  in  most  practical  situations.  The  ratio  ■—  is 

hi 

greater  than  5  in  most  industrial  applications. 

During  thermal  loadings,  a  peculiar  feature  of  crack  closure  was  observed.  It  is  an 
important  feature  when  considering  pure  thermal  loadings.  It  gives  a  characteristic  length 
of  the  compressive  zone  where  an  existing  flaw  can  remain  dormant.  Surface  cracks  in  this 
zone  will  always  remain  suppressed.  The  length  of  the  compressive  zone  is  determined  by 
a  simple  stress  analysis  of  the  uncracked  specimen.  Details  of  this  analysis  are  provided  in 
Section  5.2.  If  the  characteristic  surface  flaw  size  is  smaller  than  the  length  of  the 
compressive  zone,  an  elaborate  fracture  analysis  may  not  be  even  necessary.  It  should, 
however,  be  noted  that  there  might  be  other  loadings  superposed  over  the  thermal 
loadings  when  considering  practical  situations.  In  that  case,  the  flaw  might  propagate.  The 
resultant  stress  intensity  factor  due  to  all  loadings  should  be  considered.  In  any  case,  the 
compressive  part  of  the  thermal  loading  will  help  reduce  the  stress  intensities  for  surface 
cracks. 

The  analysis  performed  in  this  study  is  linear  elastic.  As  such,  the  stress  intensity 
factors  due  to  various  loadings  can  be  added  to  obtain  the  resultant  stress  intensity  factors. 
In  this  study,  the  stress  intensity  factors  for  individual  loadings  are  provided.  To  apply 
these  results  in  practical  situations  with  multiple  loadings,  one  needs  to  only  add  up  the 
different  stress  intensity  factors.  It  should,  however,  be  noted  that  the  results  in  this  study 
are  given  in  the  normalized  form.  Care  should  be  taken  to  take  these  normalizations  in 
account  before  adding  up  the  stress  intensity  factors. 

In  this  study,  it  is  assumed  that  the  crack  always  propagates  perpendicular  to  the 
nominal  interface.  This  assumption  has  to  be  considered  carefully  in  case  of  FGMs.  For 
homogeneous  materials  under  symmetric  loading,  the  cleavage  plane  would  be  the  plane  of 
symmetry.  This  is  the  plane  where  the  strain  energy  release  rate  G  is  maximum.  As  Gc,  the 

critical  strain  energy  release  rate  is  same  is  all  directions,  the  plane  of  symmetry  would  be 
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the  preferred  plane  of  crack  propagation.  In  FGMs,  though,  the  situation  is  different.  G 
again  is  maximum  along  the  plane  of  symmetry.  Gc  has  to  be  determined  carefully.  This 
problem  can  be  looked  in  two  different  ways.  In  the  first  case,  we  assume  that  the  crack 
growth  is  infinitesimal.  In  this  case,  the  material  being  isotropic,  critical  strain  energy 
release  rate  Gc  is  the  same  in  all  directions.  Hence,  again,  the  crack  propagation  would  be 
along  the  plane  of  symmetry.  The  assumption  of  infinitesimal  crack  growth  is,  however, 
not  very  practical.  At  the  crack  tip,  the  crack  propagation  will  be  in  finite  jumps  which 
usually  are  of  the  dimension  of  the  grain  size.  Though  the  grain  size  of  ceramics  can  be 
very  small,  the  grain  size  of  metals  can  be  large.  If  we  consider  the  average  grain  size  to  be 
about  10  fj, ,  we  can  no  longer  assume  infinitesimal  crack  growth.  This  leads  to  the 
possibility  of  varying  Gc  in  different  directions.  Depending  on  the  material,  Gc  might  be 
lower  in  a  direction  not  perpendicular  to  the  interface.  Crack  propagation  will  take  place 
in  a  direction  where  G  —  Gc  is  maximum  and  it  is  possible  that  this  direction  might  not  be 
the  plane  of  symmetry.  Again,  a  careful  metallurgical  analysis  is  required  to  determine  Gc 
accurately. 

Another  limitation  of  this  analysis  is  that  plasticity  has  not  been  considered.  It  is 
readily  acknowledged  that  plasticity  may  play  a  major  role  in  the  metal-rich  region  of  the 
FGM  and  in  the  substrate.  However,  it  is  assumed  that  the  crack  growth  considered  here 
is  subcritical  where  the  small  scale  yielding  at  the  crack  tip  does  not  affect  the  large  scale 
K -field.  This  analysis  is  useful  in  creep  and  fatigue  studies.  Rapid  crack  propagation 
affected  by  plasticity  will,  however,  be  out  of  consideration  here.  The  reason  plasticity  has 
not  been  included  in  the  analysis  is  because  of  the  complexity  it  introduces.  A  pure 
continuum  analysis  including  plasticity  for  a  problem  of  the  scale  considered  in  this  study 
is  analytically  intractable.  It  is  suggested  that  such  problems  be  solved  using  Finite 
Element  Analysis.  The  results  of  this  study  would  be  good  benchmarks  for  Finite  Element 
community. 


205 


6.2  Future  work 


The  discussions  in  the  previous  section  leads  to  the  future  work  to  be  done  in  this 
field.  The  following  problems  are  some  of  the  problems  to  be  considered. 

L  Introduction  of  plasticity 

As  discussed  before,  plasticity  would  play  a  major  role  in  stress  analysis  and 
fracture  characterization  of  FGMs.  It  is,  however,  very  difficult,  if  not  impossible,  to 
introduce  plasticity  in  a  continuum  model  for  this  problem.  The  purely  linear  elastic 
problem  itself  is  highly  complex.  The  analysis  is  very  involved  and  time-consuming.  It  is 
therefore  suggested  that  Finite  Element  Analysis  be  used  to  address  the  problems 
involving  plasticity.  The  results  in  the  work  presented  in  this  dissertation  will  be  of  use  in 
elastic  FEM  analysis.  For  example,  the  crack  tip  singularities  are  only  square  root 
singularities.  Hence  only  such  elements  need  to  be  used  at  the  crack  tip.  The  bimaterial 
interface  is  only  a  nominal  one.  No  singularity  needs  to  be  embedded  in  the  elements  at  the 
interface. 

2.  DetermMtXon  of  direction  of  propagation  of  crack 

This  problem  has  also  been  discussed  in  Section  6.1.  For  every  crack  length,  it  is 
necessary  to  assume  a  finite  crack  increment.  For  every  crack  increment,  the  difference 
G  -  Gc  can  be  evaluated.  The  maximum  value  will  determine  the  direction  of  the  crack 
increment.  Though  this  problem  can  be  solved  analytically,  it  is  much  easier  to  attempt  it 
by  Finite  Element  Analysis.  The  problem  will  be  a  mixed  mode  problem.  That  will 
introduce  more  complexity  but  will  be  a  better  model  of  practical  situations. 
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APPENDIX  A 


EXPRESSIONS  FOR  VARIOUS  KERNELS 


ff32+  4a2 +  Aio.8 


f 3  x//?2+4a2-  Aia6 

2~  2 


(3  •y//92+  4a2 +  4*0:5 

2  +  2  ’ 


(3  y//?2+  4a2—  4ia<5 

~  2  +  2  ' 


(A2  -  A5) 


(«-l)(P?  +  ftpi)  -  («+l)Q2 
2api  +  (k— l)a/3 

(«-l)(gj  +  /?P2)  ~  («+1)q2 
2ap2  +  («-l)a/5 

(«-l)(pj  +  ffftO  ~  («+l)a2 

2apz  +  (i K-l)aP 


(«- l)(pj  +  /?P4)  ~  (k+1)q2 
2ap4  +  (k— l)a/? 


(A6  -  A9) 


=  (/c+l)gipi  +  (3-«)a, 


=  («+l)92P2  +  (3 -K)a, 
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h 


(K+l)qsP3  +  (3-K)a, 


64  =  (n+l)q4P4  +  (3 -K)a. 


ri  =  Pi-aqi, 


r2  =  P2  -  ocq2, 


p3  =  P3  -  czqz, 


r4  =  P4  -  aq4, 


01 

02 

03 

04 

C* 


_A_ 

(K-iy 


h 

(K-iy 


63 

(«-!)’ 


Oa 


(A10  -  A13) 


(A14  -  A17) 


(A18  -  A21) 
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d*-  c* 

~T~' 


d*+c* 


(A22  -  A25 ) 


f-G-l 

=  /  9\(t\)R\(0L,ti)  dtx  , 

J-b i 


«~1  „,2  -c/?cos(c(ii+fei))+(2c2+2d2-d^)sin(c(ti+/i1))  ^(_(fl+^l)(d_|)) 


cd(c2+d2) 


(A26  -  A27) 


5i(*i)#2(o^i)^i, 


-c(d2+c2+^)cos(c(ti+/ii))+d(c2+d2-^-)sin(c(ti+/ii)) 


cd(c2+d2) 


r~ ai  /*^2 

/  9i(ti)R3a(a,ti)dtl+  g<2{h)Rzb(ate)dh, 
J  —b\  J  Qo 


o(Hh+hi){d-h) 


(A28  -  A29) 


1  2  c/?cos(cti )-  (2c2 +2<j2+d/3)sin(cti) 

k+1  °  cd(c2+d2) 


2oft2  _c 


(A30  -  A32) 


/*“^i  /»b2 

=  /  9l(ti)R4a(a,ti)  dt!  +  92{t2)R4h{Oi,t2)dt2, 

J  — 6i  «/  02 


2  c(<f2+c2+^)cos(cti)+d(c2+d2-^)sin(cti)  (ti(<f+|)} 

«+l  a  cd(c2+d2)  6 
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(A33  -  A35) 


—2(at^  —  l) 
/c+1 


m 


/— ai  pbi 

9i(ti)R5a(a,ti)dtl+  /  g2(to)Rsb(ptfo)dt2, 
-b{  J  02 


Rsa 


g(<l(CM-f)) 


/  {(c3+cd2+3c^)(^f  )+c<52+ca2  >cos(cti) 
\  cd(c2+d2) 


{(-c2d-d3+d2  f3+c2  (3+Zdtj)(^)+d62+da2  }sin(cti) 

cd(c2+d2) 


2at2-l+«  Qt2 
2a(«+l) 


(A36  -  A38) 


■^6  —  f  9l(tl)R6a(aJl)  dt\  +  f  92(t2)R6b(oc,t2)dt2, 

J  —b\  J  02 


i?6o  —  -Rfol  +  Rqo2  +  RfaZ, 


R&al 


RSa 2 

RooZ 

R% 


a(52  [  a2')(4cd+c/3)cos(cii)+(2c2— 2d2— d/?)sin(cti)  e«l((i+|)) 
8cd(c2+d2)2+(c2+d2)/3(c?+|) 

Q/3  c  cos(cti)  d  sin(cfi) 

2  cd(c2+d2)  6  2  ’ 

a  «+5  c/3  cos(cti)-(2c2+2d2+d/3)  sin(cti)  (t](d+£)) 

~  8  «+l  cd(c2+<£2)  6  2  ’ 


R&>i  +  Rm, 
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#661 

e~ati 

2a  ’ 

#662 

— e~atit2 
/c+1  ' 

#7* 

=  /  32  O2)  #7  (0^2)  <#2, 

Jd2 

#7 

2a(t2-h2 )  alU-ho) 

~K+l 

#8 

rh 

=  /  32  O2  )  #8  (Q  >^2  )  ^2 » 

v  02 

#8 

2(l+a(t2-^2))  ,a(U—h>') 

K+l 

ei 

=  ~(2ah2  -  k)  e(~2ah°\ 

e2 

(1+  (2a/i2  +  k)(2ah2—k) ) 
2a 

#01 

=  /  32  (^2  )  #016  (O!  2 )  ^2  > 

J02 

#016 

(#7+#8)+(2Q!^2— K)(#7~ #8) 
4a 

h  =  2ae(~2ab*\ 


h  =  (2a/i2  +  k )  e("2oJl4). 


0139  -  A46) 


(A47  -  A48) 


(. A49  -  A50) 


(A51  -  A52) 


(A53  -  A54) 


(A55  -  A56) 
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#02  =  /  92(h)  R()2b(a,h)dh, 

Joq 


Rmb  =  {Rs2R7)  e~ah*- 


(. A57  -  A58) 


6*  =  k  —  2ae2  —  «e2/i  +  Kei/2  +  /c/2  +  /cei. 


(A59) 


71 


«(l+/2)  -  a(l+5i)c2 


72 


«(l+/2)  “  a(l+g2)e2 


73 


«(l+/2>  -  a(l+qz)e2 


74 


«(l+/2)  “  a(l+94)e2 


(A60  -  A63) 


<7l(£l)#03a(ct/i)  <#1  +  /  92(h) RQ3b(&M)dt2, 

~bi  J  Q/} 


#03 a  = 


—  «i26a(l+/2)  +  CV^Ci^Sa+iJea) 


•#036  = 


—  ^#016(l+/2)  ~  ^#6b(l+/2)  +  ae2(Rbl >+#66)  +  /Ce2#026 


(A64  -  A66) 


Si 


(gi-l)a  +  a(l+9i)e2  -«/i 
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^2 

(q2-l)a  +  oc(l+q2)e2  -«/ i 

8* 

<53 

(qz- l)o:  +  a(l+qz)e2  -k/i 

5* 

^4 

(^4 — 1)0;  +  a(l+q4)e2  — k/i 

(A67  -  A70) 

5* 

■RSi 

= 

r—ai 

/  9i(ti)Ro4a(aM)  dti  +  92(h)  Ro4b(aM)dh, 

J  — 61  «/  Ct2 

#04a 

^(-R5a"~^6a)  — 

-^04  b 

KR&fi-  &(Rsb—R6b)  ~  ctei(R5b+R®j)  +  i2oib(2a+K/i)  — 

«-Ro5»(l+Cl) 

CJ3 

= 

r2^3  -  ^3^2  (-(j%-Pi)M 

r^2  -  r2*i 

(A71  -  A73) 

U>4 

= 

T204  —  ^4^2  (-(p4-pi)/ii) 

r^2  -  r20i 

(A74  -  A75) 

Ru 

= 

p—ai 

/  Pi(ti)-Ru(a.*i)dti, 

J-bi 

Rlla 

= 

R2&2  ~  R1T2  (aAl) 
ri02  -  r2^i 

(A76  -  A77) 

P3 

= 

■^3^1  —  ^1^3  (-(ps-^hi) 

r  i02  -  r20! 

P4 

= 

r4^i  -  Ti04  (— (p4— p2)/ii> 

r^2  -  r2*i 

(A78  -  A79) 
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R-h  —  f  9i(ti)Rn(,aM)  dti, 

J-bi 

R “■  “  r2-M  £<R',1>-  (^80-^81) 

ii  =  2a  +  2aei-  (k+1)/i, 

j*2  =  («+l)  +  2ae2  —  (k+1)/2, 

*1  =  <f>i  +  ji7i  + 

^2  =  4*2  +  Ji72  +  ^2^2-  (-A82  —  .485) 

✓'"N  ^  ^ 

-^3a  =  R^a  j\Ro3a  ,?2-^04c  ^l^lla  ^2-Rl2a» 

-R&  =  -R&  -  Ji-Ro3fc  -  j^Roib  ~  2ai?oit  +  («+l)i?02a.  (.486  -  487) 

l\  =2 a  -  2aei+  (ac-1)/i, 

^2  =  (k-1)  -  2ae2  +  («-l)/2, 

mi  =  +  Zj7i  +  Z^i, 

m2  =  r2  +  Z|72  + 1*282.  (488-491) 
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/N 

R-4a 


Rtfj  l\Ro&a  ^2  -^04c  ^l-Rlla  rffl't2lRvicLi 


R&  =  —  l\R^sb  —  ^2-^04 6  +  2ai2oife  ~  (k— l)i?026-  (^492  —  A93) 


^3  =  03  +  il73  +  J2^3> 

^4  =  04  +  il74  +  j^4. 


m3  =  73  +  ^73  +  ^3, 


m*  =  7*  +  Z*l74  +  Z^4. 


(.494  -  4.97) 


Pi  =  fcfu^  +  k$p3  +  fc|, 


$  =  fcJo>4  +  k$p4  +  fc|, 


q{  =  m\u 3  +  m^/93  +  m^, 


§2  =  mju>4  +  m2/?4  +  m4. 


(498  -  4101) 


/— Gl 

giit^RiSaiocM)  dt{  +  /  g<2{t2)R\3b(ate)dt2, 

•b\  J  Q% 


R]3a 


R3aQ2  R4aP2 
PlQ2  ~  P2<ll 


Risb 


RmZ  —R-4bp2 

Pitf-Pffl 


(. A102  -  4104) 
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R* 


14 


f—a  1  rbi 

I  9l(tl)Rl4a(a,ti)  dti  +  /  gz(t2)R-],4b(a,t2)dt2, 

J  —b\  J  02 


Rl4a  = 


Rub  = 


RaoP j  -  Rsaqj 


P*l<£  -  PM 


R^bPi  -  -R36gl 

pM  -  PM 


$1  ={—('— (3— K)miia+(K+l)ni)eniy  m277,2 — — 

^  V  /  7711^1  7712712 


+  (3-«)m2«o:+(K+l)ra2^en23/  j-e-10 


3/— >0 


z\  =  <(l+/c)a  +  (3-/c)pigi>  e^2^, 


-^2  =  {(l+«)a  +  (3-/c)p292>e(P2:El), 


23  =  <(l+K)o:  +  (3-K)p3Q'3>e(^a:i), 


•24  —  {(1+k)o:  +  (3— «)P494  }  . 


r—a\  f>k>2 

•^21  =  /  5l(il)-R21a(Q!*^l,^l)  ^1  +  /  g2{t2)R21b(.Ct>Xl,t2)dt2, 

J  —bi  J  G2 

•R21a(^^l,^l)  =  2  Re  +  i?i3a(23+  ^3^1+  P3^2)  j, 


(A105  -  ,4107) 


(A108) 


(A109  -  A112) 
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R2ib(a,xit2) 
kl  ={- 
+ 


Kn  (a,xi) 


Kn(a,x  i) 


{ 


=  2  Re  <  i?l3a(^3+  ^3*1+  P3&2) 


(A113  -  >1115) 


-  (— (3— K)miia+(K+l)ni)eniy  m2Tl2 — — 

\  /  m\Ti\  —  m2n2 


(-(3-«)m2ta+(K+l)n2)en23/}| 

y— *0 


(mini  —  ia)  i 

(mini  —  m2n2)  g: 


(A116) 


=  jjf  {i?e(fci(a,y))cos(Q!(ti -xi)) 


-  (/m(fci(a,y))  +4^^^^sin(a(i1  -xi))}da 


+ 


J r  {  (-Re(fci(a,y)  -  i?e(fci00/l(a,y)))cos(a;(ti  -  xi)) 
-  ^7m(fci(a,y))+4(^Y)  - 


Jm(fci00^(o:,  y))^  sin(a(ii  -  xi))}da 


{Reikioohtd, y))cos(a(ii  -  xi)) 

-  Im  (fci007l(a,y))sin(o:(t1  -  Xi))  }<2a 


+  2^  J?2ifl(a,a:i,ii)da  }  (4^”!))  04117) 

=  2 jf  iJ2ii>(a,xi,i2)dQ!  (>1118) 
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(4119) 


2*(^  +ayy-iQXt~W 


\y 


y->0 


k 2 


4 i  /  |a 


(1  +  k) 


-(^+ay)e"l^ 

k)  V  a  / 


(4120) 


y-*> 


n\  =  2a(  2a(x2-h2)+3)  ea(X2~2h)  +  2aeax\ 


n\  =  ( ( 2ah2+K)(  2a(x2-h2)+3)-l )  eQ(*2-2/l2)  +  (2arr2+K-3)eQX2, 


01  =  ^3(7i^i+^in2)  +  P3(72^i+^2^)  +  73^1  +  hn2, 


O 2  =  k>4(7lni+6ln2)  +  P4  (72^1+62^2)  +  74^1  +  64^2*  (-A121  -  ,4124) 


/‘-ci  pb2 

R^22  =  /  9l(tl)R22a(Oi,X2,t1)  dtx  +  g2(h)R22b(ct,X2,t2)dt2, 

J  —6]  J  02 


R21(a,x2,t\)  =  -Riia(7ini+^ira2)  +  -^120(72^1+^2) 


+  i?13a°!  +  i?14a°2  +  -Ro3anl  +  Ro4an2 » 


(-R7+.R8)+(2Q!/l2~^)(-fi7--R8)  -Q(/i2-i2) 
2  6 


+ 


(2qj;2  — K+3X-R8— -R7) 


+  -Ri3jo|  +  i?i4&o|  +  i?03i,7ij  +  R§^bn2- 


poo 

K2\  (oc,x2)  =2J0  R2ia(a,x2,ti)da\ 


k  +  1' 


(A125  -  A127) 
(A128) 
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Kn(atx2)  =2  jf  R2ib(afx2,t2)da  ^ 


k  + 1 
4 


(A129) 
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APPENDIX  B 


ASYMPTOTIC  EXPRESSIONS  FOR  KERNELS 

The  expressions  for  the  various  variables  used  in  Chapter  3  in  the  aysmptotic  expansions 
of  the  integrands  of  the  Fredholm  kernels  are  given  here. 

•^21a(o!»  ^l) 

Rzia(ot,xi,ti)  ~  (cla2  +  +  4  +  ^e“(zi+tl+2/ll)a 


+  (e;o  +  ^  +  f  +  §)e<w+,‘>“ 

(Bl) 

where 

4 

-  16  (3_k>  p  sm(Ci)sm(C2)e(  2  } 

m 

K  —  l  1  / 

4 

=  4  (3_k)(1+«)  (3  V  2£ism<£i)sm(C2)  +  Acos(C2)sm(Ci)(l+«) 

+  3Asin(C2)cos(Ci)(l+«)  +  6sin(C2)sin«i)(l+«)^e(  2  } 

(S3) 

4 

=  2(3-K)(l+K)2  (d«nK<i»n(6) 

+  £i(H-«){  Acos«2)sin(Ci)+  3Asin(C2)cos(Ci)+  6sin(C2)sin(Ci)  } 
+  (1+k)|  -2(k-3)cos(C2)cos(Ci)  +  A(l+K)cos(C2)sin(Ci) 
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+  5A(l+K)sm(C2)cos(Ci)  +  6(l+K)sin(C2)sin(Ci)  +  16sin(C2)sin(Ci)  j 


-  2^3Asin(C2)cos(Ci)  -  2£2  Acos(C2)sin(Ci))  e(W  '2  1>}  (54) 

k—1  f3 
-  (3— k)(1+«:)3  6 

(2^f  sin(C2  )sin(Ci ) 

+  ( 1 + k)  j  Acos(C2  )sin(Ci )+  3Asin(C2)cos(Ci)+6sin(C2)sin(Ci)| 

+  6^1(1+«)|-2(«-3)cos(C2>cos(Ci)  +  A(l+«)cos(C2)sin(0) 

+  5A(l+K)sin(C2)cos(Ci)  +  2(7+4/c)sin(£2)sin(Ci) 

-  2£3Asm(C2)cos(Ci)  -  2^2 Acos(C2 )sin(Ci )  j 

+  6(1+k)  |Acos(<^2)cos(Ci)  +  16A(l+K)sm(£2)cos(Ci) 

+  (fc+41)(l+K)sin(C2)sin(Ci)  +  £3(/c-3)cos(C2)cos(Ci) 

-  ^2 (« — 3 )sin(C2 ) sin(Ci )  +  3^2(k-3)cos(C2)cos(Ci) 

-  3|3(/c-3)sin(C2)sm(Ci)  -  6^2(l+«)Acos(C2>sin(Ci) 

-  6£3(l+«)Asm(C2)cos(Ci)  J  (55) 
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( 2hi+x\+t\')/3, 


( hi+ti)(3, 

( hi+x{)/3.  (B7  -  B9 ) 

l  A/5(  hi+xi), 

\  A/?(  fti+ti), 

|  A/?(  xi), 

jA/3(ti).  (BIO- £13) 

-8<3 (B14) 
-4  (3-^X1+^ 1 *<“ 

(sin(C3)  ({s+2k)  +  Acos(C3>  (l+(t))e(  (£15) 

«—  1  /? 

~  (3— k)(1+k)3  2 


(sm(C3)sin(C4)  {  2£§+8k£5+k3-k2+7/c+25} 

-cos(C3)cos(C4)  {  (1+k)2(k-3)  >  -  sin(C3)cos(C4)A(l+K)3 
-cos(C3)sin(C4)A(l+K)(K2-2K+5) 

-4sin(C3-C4)A(l+K)/?t1)e(=^)  (£16) 

K — 1  01 

(3— k)(1+k)4  12 

(sin(C3)sin(C4)  {2£i+12/c£§+3/3*i(K3-/c2+7K+25) 

+3/3xi(k3— 5k2+15k+37)+3(k4— 2k3— 4/c2+10k— 21)  j 
-cos(Cs)cos(C4)  {3^i(1+k)(k-3)2+3/3x1(1+k)2(k-3) 

+3(k2  — 1)(k— 3)(1+k)  j 

-sin(C3)cos(C4)A  |l2/52x1t1(l+K)+12/?2^(l+K) 

+3/3ti(l  +  K)(K2+10K+l)+3^5(l+K)3+3(l  +  K)2(K2+3)| 
-cos(C3)sin(C4)A  |6/32x^(1 +k)-  6/32£2(1+k)  +  3/5£i(1+k)(k2-2k+5) 
+3^5(1+k)2(5+k)  +  3(k2-1)(k2+6k+13) 

(517) 
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where 


£5  =  (xi+ti)/3.  (B 18) 

Rwbfa  »i.  *2)  ~  (ria  +  r2  +  ^  +  ^l)  e(xi_t2)£i  (-B19) 

where 

r\  =  ^A(-2+k)2  (sin(C3)(^2/5(«-l)2+2(«2-l))-cos(C3)^2^(«2-l))e(^) 

(B20) 

r2  =  (/5^i{  sin(C3)02/5(«-l)2+2(K2-l))  -cos(C3)^2/5(«2-l)} 

+sin(C3)  {  2i2/?(/c+3)(/c-l)+2(«2-l)(l+K)  } 
-cos(C3)A{2t2/?(«2-l)-4(«2-l)(l+«)  })e(r^>  (521) 

r$  =  2A(1+k)4  (^2j;l{sin(C3)(^2A«-l)2+2(«2-l))  -cos(C3)^2A«2-1)} 

+2/3*! {  sin(C3)(t2/?(K2+3)i^i11+(«2-l)(l+«)) 

+cos(Cs)  (-t2/3(K2-l)^+(«2-l)(l+«))  } 
+sin(C3){2i2/?(«2+3K-2)(«-l)+(K2-l)(l+«)2} 

+cos(C3)  A  |  -2  £2/?(/c2  -  l)(/c-  1)+(k2  -  l)(l+«)2 e(_^) 

(522) 
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rt  =  4A(1+k)5  (/5M{sin(C3)(t2/5(K-l)2+2(K2-l))-cos(C3)i2/5(K2-l)} 
+2/32x\ |  sin(C3)(t2/5(«2+3)^+(K2-l)(l+K)) 

+cos(C3)  (-t2/5(«2-l)^+(«2-l)(l+«))  } 
+/?xi{sin(C3)(2t2A2«2+K-5)(/c-l)-(K2-l)(l+K)(3K-13)) 
+cos(C3)A(4f2A«2-l)(«+l)+(«2-l)(l+«)2 )  } 
+sin(C3){2i2/5(«3+K2-llK-3)(/c-l)-(/c2-5)(K2-l)(l+K)  ] 
-cos(C3)a{  2i2^(«2 — 1)(k— l)2 +(k2 — l)2 (1+k) 

(S23) 


•R22a(o:>a?2^l) 

i?22a(a,  *2,  ti)  ~  +  ^2  +  f  +  ^  +  5)  e~(X2~tl)a  (£24) 

where 

*1  =  -^(k+I)2  (”2cos^)Aa;2/5(K+l)+sm(C4)(2x2/3(«;-l)+4(K+l))) 

(£25) 
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*2  =  ~  A(kH~1.)3  (^{  _2cos(^)Ax2/5(«:+l)+sm(C4)(2a:2/5(/c-l)+4(K+l))J 

+  2cos(C4)A(/c+l)2+sm(C4)(2x2/3+(l-K2»  j  e("^ 

(1326) 

<3  =  ~  4A(f+1)4  (/^  -2cos(C4)Ai2/3(K+l)+sm«4)(2i2/8('t-l)+4(K+l))} 

+  4/3ti|2cos(C4)A(K+l)2+sin(C4)(2a:2/?+(l-«2))  j 

—  2cos(C4)A(k+1)3—  2sin(C4)^2^iA2/52(«+l)2 
+  sin(C4)(4x2/3(K-3)  -  2(k+1)(k2-2k+5)) 

-  cos(C4)Ati/3(K+l)(2a:2/5(K-l)+4(K+l))j  e(^ 

(£27) 

t*4  =  ~gA(^+1^  (^i{  -2cos(C4)Ax2/3(«:+l)+sin(C4)(2rE2/3(«;-l)+4(«+l))} 

+  4/52^  j  2cos(C4)A(«+l)2+sin(C4)(2x2/?+(l-/c2))  | 

+  /0ii{-2cos(C4)A(K+l)3+sm(C4)(4x2/3(K-3)-2(K+l)(K2-2«+5))  } 

+2cos(C4)A(k+1)2(2x2/?+(k2-1)) 

+  sin(C4)(4x2/5(K2  -8«- 1)+2(«2  -1)  (k2  -2k-  11) 
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-sin(C4)(2x2ii/?3A2(«:+l)2-8A2ySt1(K+l)3 


— cos(C4 )  |^  i/52  A  (k+ 1 )  (2^2  /3(k — 1)  +  4(k+1)) 

+  4(K+l)t1/5A(l —k2 +2x2/3)  }  )  e(^)  (B28) 

t*5  =  ~  16A^.+1)6  (/^^i{-2cos(C4)Ax2/3(K+l)+sm(C4)(2a:2/?(K-l)+4(K+l))| 

+  4/33i3{  2cos(C4)A(«:+l)2+sm(C4)(2a:2/5+(l-/c2))  | 

+  /?2^{  -  2cos(C4)A(k+1)3+  sin(C4)(42:2/3(K-3)-2(K+l)(/c2-2K+5))} 
+  /?£i  j  2cos(C4)A(K+l)2(2a:2/5+(«2-l)) 

+sin(<^4)(42:2/5(K2— 8/c— 1)+2(/c2— 1)(/c2— 2k— 11)  j 
+  2cos(^4)A  (  k(k+1)3(k— 3)— 2x2(3(k2— 1)(/c+1) 

+  2sin(C4)(4x2i0(5K2-8K-5)+(K+l)(K4-5K3+3K2-7K+16)) 

+  cos(£i) j  X2t2A3/33(K+l)2+4^A/52(l-/c2+2x2/9) 

+i1A/?  («+!)( 4x2/3(«-2)-(k+1)(«-1)2)  } 
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+  sin(Ci){  t2A2/32(/c+l)2(  X2/3(/c-l)-6(/c+l)  ) 


+2t|A2/?(/c+iy*+2t3X2A2/53(K;+l)2  ^  ^ 


•R22&(«>®2>*2) 


(529) 


R22b(^x2,t2)  ~  (oi<2  +  +  ^|)e 


-(rc2+^)^ 


+  ^Pia2  +  P2<2  +  P2^j  e'X2+<2  2/l2)Q 


where 


(B30) 


2X2^/? 

(K+l)2 


(531) 


2*2/3  X2*2/32(«-l) 

(«+l)2  +  (k+1)3 


(532) 


/3  *2/32(l-«)  x2*2/33«(ac-3) 

2(«+l)  +  (k+1)3  +  2(«+l)4 


(533) 


£2(1-k2)  *2/33(l-«)2 

4(«+l)3  2(/c+l)4 


x2/33  X2*2/34(«— 1)(«2— 5k— 2) 

2(k+1)3  +  4(«+l)5 


(534) 


/33k(k-3)  t2/34  (k2 —1)(k2 —4k  - 1)  x2/34(k-1) 

8(k+1)3  4(/c+l)6  4(«+l)4 


X2*2/35(k4-10k3+14k2+10k+1) 

8(k+1)6 


(535) 
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(^2-^2)(^2-^2) 

(«+l) 


Pi 


P2 


2 


(^2—^2)  +  3(^2 — ^2) 

<*+Ij 


P3 


4 

(«+!)* 


(jB36  -  S38) 
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APPENDIX  C 


RIGHT  HAND  SIDES  AND  NORMALIZATIONS 
FOR  DIFFERENT  LOADINGS 


As  can  be  seen  from  equations  (3.28),  (3.36)  in  Chapter  3,  the  right  hand  side  of  the 


singular  integral  equation  is 
RHS  =  pi(x ) 


(g+jO 

Ap(x) 


where  p\(x)  is  the  crack  surface  loading  corresponding  to  different  far-field  loadings.  Due 
to  superposition,  it  turns  out  that  this  p\(x)  is  nothing  but  the  negative  of  the  normal 


stress  produced  at  the  plane  of  the  crack  in  an  uncracked  specimen.  In  this  appendix,  the 
expressions  for  right  hand  sides  and  the  normalization  for  the  stress  intensity  factors  are 


derived  for  fixed  grip  loading,  far-field  tension,  far-field  bending  and  thermal  loading.  All 
strains  and  stresses  are  evaluated  at  y  =  0  plane. 


Figure  C.l  Nonhomogeneous  strip  bonded  to  homogeneous  strip.  Internal  crack  in 
nonhomogeneous  strip. 
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Consider  a  strip  as  shown  in  Fig  ((7.1).  The  strip  is  homogeneous  on  positive  side 
of  x  and  nonhomogeneous  on  the  negative.  Let  it  have  a  crack  extending  from  a  to  b.  An 
edge  crack  can  be  obtained  by  letting  a  go  to  zero.  Let  the  unknown  function  (derivative 
of  the  displacement)  be  g(t)  and  the  corresponding  function  in  normalized  coordinates  be 
g(r).  Let  g(r)  be  defined  as 

g(r)  =  F(r)w(r )  ( C2 ) 

where  w(r)  is  some  weight  function.  F(r)  in  general  will  be  a  summation  of  the  sort 

F(r)  =  Y/AnX„  (C3) 

n 

where  Xn  is  some  orthogonal  polynomial  and  An's  are  the  unknowns. 


Fixed  grip  loading 

In  this  case,  we  apply  a  uniform  strain  eo  at  infinity.  Stress-strain  relation  gives  us, 

_  («+l)  /  _ 

€yy~  8fl(x)[  yy  1 +K* 

For  this  loading, 

£zz  ~~  ^ xx  ~  0  ■ 

Hence,  using  equation  ((74) 

=  ~Pi(x).  ((76) 

Substituting  into  ((71) 

RHS  =  -  2e0  ((77) 

and 


m 

(C5) 
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(08) 


ki(a)  = 


4m(«) 

(k+1) 


<»» 

It  should  be  noted  that  for  an  edge  crack,  the  expression  ^  °  is  repalced  by  \/&. 

Dividing  the  entire  equation  by  2eo,  we  can  now  normalize  the  stress  intensity  factors  with 
respect  to 


8^o  lb  —  a 

ao  =  (k+1)  V  2  f0' 

(CIO) 

Then  we  obtain, 

RHS  =  -1, 

(Oil) 

M<0  =  —  f  ( - 1), 

^0 

(C12) 

hnW  =  ~ 

^0 

(CIS) 

It  should  be  noted  when  a  goes  to  zero,  we  have  an  edge  crack  at  the  left.  For  an  edge 
crack,  the  normalization  is  usually  done  with  respect  to  the  shear  modulus  y,\.  This  yields 
the  benchmark  value  of  1.1215222  for  the  normalized  stress  intensity  factor  in  the  limiting 
case  when  the  crack  length  goes  to  zero.  To  obtain  normalization  with  respect  to 
i,  it  is  necessary  to  multiply  eqn  (C13)  by  a  factor  of  — . 


Far-field  tension  and  bending 

Usually  any  far-field  loading  can  be  resolved  into  a  far-field  tension  (N)  acting 
along  the  midplane  of  the  specimen  and  a  far-field  moment  (M).  If  the  crack  is  sufficiently 
away  from  the  loading,  by  St.  Venant's  principle,  the  stresses  produced  by  the  resultant 
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force  and  moment  will  be  the  same  as  those  produced  by  the  specific  loading  for  all 
practical  purposes.  Under  such  loading,  due  to  compatibility,  we  have, 


6yy  -  Ax  “b  B. 

(C14) 

Equilibrium  equation  states  that 

J  (Jyy{x)dX  =  N 

(C15) 

while  moment  equation  gives, 

/■**,  h-h 

(x  9  )ayy(x)dx  -M. 

J  ~h\  " 

(C16) 

Using  Hooke's  law, 

_  /_x  M*)  £ 

w( )  "  (K+l)  w 

(C17) 

and  equations  {CIA),  ((715)  and  ((716),  some  simple  algebra  yields, 


e-ph i  2e~Phl  2  2e~@hl  h\ 

(3  ~  +hp5_"hp3-  +  3hf)’ 

e~@hl  1  e~Phl  h\ 

h\0>  +  h\p  +2hf}' 


(CIS) 

(C19) 

(C20) 

(021) 

(C22) 

(C23) 
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Solving  (C18)  and  (Cl 9)  we  get, 


(C24) 

(C25) 

(C26) 

(C27) 

(C28) 

(C29) 

(C30) 
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Substituting  ((726)  into  (Cl)  we  have 
RHS  =2  (Ax  +  B). 


(C31) 


In  case  (I)  M  =  0,  dividing  the  entire  equation  by  2  N,  we  can  normalize  the  stress 
intensity  factor  with  respect  to 


8uq  lb  —  a  fv  N  /b  —  a 

ffJ,"(«+l)  V  2  W=/!1V  2 

(C32) 

giving, 

Ax  +  5 

=  Jj  , 

N 

(C33) 

km(a)=  ^F(-l), 

/^o 

(C34) 

M»)  =  - 

H  0 

(C35) 

In  case  (II)  N  =  0,  dividing  the  entire  equation  by  2  M,  we  can  normalize  the  stress 
intensity  factor  with  respect  to 


aM  = 


8^0 

(«+ 1) 


M  lb -a 

hi\l  2 


giving, 


RHS  = 


Ax  +  B 
M 


(C  36) 


(C37) 
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Thermal  Loading 

The  thermal  stress  problem  has  been  discussed  in  detail  in  Section  5.2.  The 
expressions  for  Ci(i  —  1, ..,  4)  are  as  given  in  ((720  —  (723).  For  uniform  temperature 
change  and  for  conduction  case,  the  expressions  for  R\  and  R%are  different.  Hence  the 
two  cases  are  considered  separately. 


Uniform  temperature  change 


R\  and  R%  are  given  as, 


Ri  = 


ftp  AT 
(a  +  P)hi 


(l  -  e"(Q+^Al)  +  a0AT^, 


The  right  hand  side  can  now  be  written  as 


((740) 

(041) 


RHS  = 


2a0A:r(-e“  +  |gfi),i<0 

2a0AT(^),rr>0 


(042) 


Dividing  the  entire  equation  by  2ap  AT  the  stress  intensity  factor  can  now  be  normalized 
with  respect  to 
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C Jj' 


SjUpapA  T 

1  +  K 


(CAS) 


giving, 


RHS  = 


(-«“  +  4g?).*<o 
(£$).*>» 


((744) 


M®)  =  —  F(  -  1), 

Mo 


((745) 


M&)  = 


mW 

MO 


^(1). 


((746) 


Stead#  state  heat  conduction 

R\  and  Ro  are  given  as, 


Ri 


<yp(AT-[  -  AT{)  /I  -  \  /A7^-A7\  wl-e^^N 

/ii(a  +  /5  —  7)  \  e^1  —  1  /  a°  V  e^hl  —  1  1 )  \  h-[(a  +  f3)  ) 


+  a»(  -  ATi)^-  +  ao(ATj)^ 

Z/Zl  All 


((746) 


i?2  = 


Qq(ATi  -  ATj)  /  (Q+^  _  1  -  \ 

(a  +  (3  —  7)(e'>*1  —  1)  \  (a  +  (3  —  7)  / 


+  Qo 


/ATj-ATi 
V  e^i  -  1 


+ 


e-(a+/3)A, 
/*i(a  +  /3) 


1  _  e-(a+0)hi  \ 

/if(a!  +  /5)2  / 


+  <*o(  -  AT*) 


hi 

Sh\ 


+  a0(ATl) 


hi 

2  h\ 
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(C47) 


The  right  hand  side  can  now  be  written  as 


j2coAr1(  + 

RHS  =  <  ) ,  dN  1  J 

|  2ao  ATi  ^  J  i  x  >  0 

(C48) 

Dividing  the  entire  equation  by  2ao  ATi  the  stress  intensity  factor  can 

now  be  normalized 

with  respect  to 

8/ioQ;oATi 
<7T=  l  +  K 

(C49) 

giving, 

((750) 

M«)  =  —  F (  -  1) 

Vo 

(C51) 

M»)=  -  —  ffl) 

Vo 

(052) 
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APPENDIX  D 


SOME  USEFUL  INTEGRATIONS 

During  the  numerical  analysis  performed  for  this  dissertation,  some  integrations 
were  evaluated  in  the  closed  form.  Closed  form  solutions  have  an  advantage  over 
numerical  solutions  in  that  the  time  spent  during  the  numerical  phase  is  considerably  less. 
Also,  the  accuracy  is  much  higher.  In  this  appendix,  the  expressions  for  these  closed  form 
solutions  are  given. 

In  the  analysis  of  internal  cracks  ( Chapter  3 ),  the  following  expressions  were 

used: 


In  the  analysis  of  edge  cracks  ( Chapter  3 ),  the  following  expressions  were  used: 
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dr  =  —  iy/2  +  2y/2 log\l  +  s\  +  2\/l  —  ~slog\B\ , 


s|  <  1 


(C6) 


where 
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APPENDIX  E 


MATERIAL  PROPERTIES 


This  appendix  gives  a  table  of  material  properties  used  in  the  numerical 
calculations  is  this  dissertation.  The  values  are  taken  from  Noda  et.al.  [51] 


k,  W/mK 

E,  MPa 

a  x  106,  l/K 

PSZ 

2.09 

151 

10 

Ti  -  6  Al  -  4  V 

7.5 

116.7 

9.5 

5175304 

11.9 

193 

18.7 

Rene  —  41 

25.5 

219.7 

16.7 

Molybdenum 

14.2 

310 

5.1 

Tungsten 

16.5 

400 

6.0 

Table  E.l  Thermal  and  mechanical  properties  for  some  commonly 
used  materials 
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APPENDIX  F 


INTERNAL  AND  EDGE  CRACK  IN 
NONHOMOGENEOUS  HALF-PLANE 

The  geometry  of  this  problem  is  shown  in  Fig  F.  1.  The  crack  is  located  on  the 
y  =  0  plane  from  x  =  a  to  x  =  b.  If  a  >  0,  we  obtain  an  internal  crack.  If  on  the  other 
hand  a  =  0,  we  obtain  an  edge  crack. 


Figure  F.l  Geometry  of  the  problem 

Formulation  of  the  problem 

The  governing  equations  for  the  nonhomogenous  medium  were  developed  in 
Section  2.3.1.  They  are  as  follows : 
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d^u  d^u  d^v  du  dv 

(1+Kfe  +  (K_1)a^  +  2fad~y  +  Pa+K)ai  +  m~K)ai  =  °’ 


(FI) 


d^v  d^v  d^u 

(k_1)_  +  (1+k)_+2_ 


dv  du 

+  /3(K-1)-+/3(«-1)-=° 


(. F2 ) 


where  /?  is  the  nonhomogeneity  parameter.  We  use  the  superposition  technique  to  find 
expressions  for  stresses  and  displacements  in  a  manner  similar  to  that  discussed  in  Section 


2.3.3.  We  superpose  the  solutions  obtained  from  a  nonhomogenous  plane  with  crack  and 


those  obtained  from  a  nonhomogeneous  half-space  without  a  crack.  The  expressions  for 


stresses  and  displacements  thus  obtained  are  as  follows: 


1  f°°  ^ 

u  =  —  y~'mJDj(a)  e^e-^da  + 

MR  J  —00 

9  r°°  ^ 

-  /  y ^qjAj(a)  ePjXcos(ay)da,  ( F3 ) 

nJo  j-  l 


i  roc  2 

V  =  —  /  ^Dj(a)  e,We“ioxda.  + 

J  _oo 

9  /*°°  _  £ 

-  /  V Ai(a)  ePjIsin(ay)da, 


(F4) 


O rr  — 


rtne(^x)  1  f00  r  'i 

- —  —  /  /  ^  — (M+l)m,io;+(3— «)nj  fjD7(a)eniJ/e-tQ;Eda  + 

K-l  27r J-oo  4-^  t 

^K-l  x/  ^ 1  )gjPj+(3 - «)a | ^4 7(a)eP;j:cos(Q:y)dQ> 


(F5) 
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ayy  ~ 


p{px)  i  r°°  *  r 

^  J  ^  {  -  (3- K)miia+(K+ 1  ^  I  Di(a)en^e^“:da  + 

(#r)  2  Z*00  ^ 


^  J  y;{(3- K)gjPj+(K+ 1  )a | A j(o;)ep,a:cos(ay)dQ;, 


(F6) 


i  r°°  ^  ^  n 

J  y;{-ia+mjnJ}^(a)en^e— da  + 

2  r00  ^  > 

fi oe^—  /  y^^Pj—gjQ  [•Aj(a)eP;3:sin(ay)da.  (F7) 

Expressions  for  p/s,  g/s,  rrij's  and  rij's  are  provided  in  Appendix  A  and  in  Chapter  2.  The 
unknowns  Ai,  A2,  -Diand  F>2are  evaluated  from  the  boundary  conditions  which  can  be 
expressed  as 


<?xx(o,y)  =  0, 

(F8) 

II 

O 

(-F9) 

Txy  (x,  0)  =  0, 

0<  x  <  00, 

(F10) 

&yi /(•c>^)  =  P(*^)» 

a<  x  <b, 

(Fll) 

u(z,0)  =  0, 

0<  x  <a  ,  b<  x  <  00. 

(FI  2) 

Development  of  the  Integral  Equation 

Using  (F7)  and  (F 10),  it  may  be  seen  that, 


1  r°°  ^  (  \ 

fjboe^—  V  -ia+mjnj  F/a^V^da  + 

2*J~°°  J-l  ^  J  y— >0 


2  Z*00  ^  f  N 

poe^x)  -  /  ^  <  Pj-qjCi  }A/a)e**sin(ay)da 

TTJ-oo  j=!  1  J  y^0 


0.  (FI  3) 
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Taking  now  the  Fourier  inverse  transform,  we  have, 


Di(a)  =  - 


(m2n2  —  ia) 


(mini  —  ia) 
We  define  new  unknown  function, 


F>2  (a). 


(F14) 


,  ,  d*u  .  . 

9(x)  =  ^(x’°) 


a  <  x  <  b. 


(F 15) 


Now  taking  Fourier  inverse  transform  of  (F4)  and  using  (F14),  (F15)  we  obtain 


£>2(0;) 


(mini  —  ia)  i 
(mini  —  m2n2)  a 


f 


gfye^dt. 


(F 16) 


By  using  the  two  homogeneous  boundary  conditions  (F 8)  and  (F9)  and  the  expressions 
from  (F5)  and  (FI),  after  taking  the  appropriate  Sine  or  Cosine  transforms  we  can  write  a 
system  of  two  algebraic  equations  in  two  unknowns.  The  system  of  equations  we  obtain  is 
as  follows: 


anAi  +  anA2  =  Rl  (F 17) 

&2lAi  +  <222-42  =  R%  (F 18) 


where 


on  =  (/c+l)giPi+(3— /cK  (FI  9) 

<212  =  (k+1)52P2+(3— K)a,  (F20) 

ail  =  Pi-qia,  (F21) 

G22  =  P2-<72<x  (F22) 


R\=  l  g(t)Ri(a,t)dt, 
Ja 


(F23) 


k— 1  2  (-t(d-i))  (  — c/5cos(ct)+(2c2+2d2— d/3)sin(cf) 

k+1  “  6  2  V  cd(c2+d2) 


(F24) 
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RZ  = 

[  g(t)R.2(a,t)dt , 
Ja 

(F25) 

i?2  = 

2  n  -(-*(<*- f »  ( 

■  — c(d2 +c2 + ^  )cos(ci)+d(c2 +d2 — ^  )sin(ct) 

A 

G:  6  *  1 

AC+l  ' 

cd(c2+<22) 

)' 

(F26) 

C  = 

/  <2*-  c* 

V  4  ’ 

(F27) 

d  = 

/  <2*+  c* 

V  4  ’ 

(F28) 

c*  = 

2  /?2 
cr  + 

A 

(F29) 

d*  = 


(a2  +  ^-)2  +  a262, 
4 


(F30) 


(F31) 


The  pair  of  algebraic  equations  in  (FI 6)  and  (F17)  can  be  solved  using  Cramer's  rule. 


The  solution  thus  obtained  can  be  expressed  as 


Ai  = 


0^2  R\  —  a  12  Fo 
an  022  -  °12a21  ’ 


(F32) 


M  = 


ani?2  —  asiFi 
°11°22  —  012021 


(F33) 


By  using  the  mixed  boundary  condition  in  (Fl  1),  and  using  (FI 4),  (FI 6),  (F31),  (F32) 
we  can  now  write  the  integral  equation  as 
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— J  g(t)dt J  |jRe(/ci)cos(a(t  —  x ))  —  Im(k\)sm(a(t  —  x ))  j da 

j  pb  poo  ^  K  —  1 

+  -/  9(t)dt  k2(a,x,t)da  =  p(x) - g~ 

nja  JO  VOeP 


(F34) 


where 


k\(a,x,t)  =  — 


^(1  +  k)ui  —  (3  —  K)iamij  (m2n2  —  ia ) 


ia(m2n2  —  mini) 


+ 


( (1  +  k)u2  —  (3  —  K)iam2)(mini  —  ia) 
ia(m2n2  —  mini) 


k2  (a,x,t)  = 


((1  +  k)q  +  (3  -  K)qipi)J  (<222 -Ri  -  ai2R2)ePl0 


ail  022  —  012021 


(F35) 


f(l  +  fc)a  +  (3  -  k)?2P2))  (oh-R2  -  a2iRi)e*x 
+  ± - J- - .  (F36) 

011022  —  O12O21 


Extracting  the  singular  part  of  the  kernel  in  a  manner  identical  to  that  in  Section  3.1.3,  we 
obtain, 

<F37> 

where 


K (*,  t)  =  \  (fQ  {-Re(fci)cos(a(t  -  *)) 

—  ( Im(ki )  +  4^ ^  )sin(a(f  —  x))  +2k2(a,x,t)'^daSJ . 

(F38) 
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It  should  be  noted  here  that  K  ( x ,  t )  contains  a  generalized  Cauchy  kernel  as  was  found  in 
Chapter  3.  To  extract  the  generalized  Cauchy  kernel,  we  need  to  perform  asymptotic 
analysis  on  k%  (a,  x,  t)  as  a  — ►  oo.  The  behavior  of  (a,  x,  t)  as  a  — ►  oo  can  be 
expressed  as 

g(a, x, t)  -  (cjo:2  +  c$a  +  4  +  ^) e(x+t)Q •  (^39) 

The  expressions  for  c*'s  are  as  given  in  Appendix  B  where  x\  +  hi  is  now  replaced  by  x 
and  1 1  +  hi  is  replaced  by  t. 


Singularities  at  crack  tips 

Internal  crack 

For  a  >  0,  we  obtain  an  internal  crack.  This  case  is  identical  to  that  discussed  in 
Section  3.3.1.  It  should  be  noted  here  that  the  coordinate  axis  is  shifted  in  this  analysis. 
But  that  does  not  affect  the  evaluation  of  the  singularities  in  any  way.  We  find  that, 

=  |  (F40) 

and  the  unknown  function  can  be  expressed  as 


9(1)  = 


m 


(t  —  a)*(b  —  t)s 
where  g*  (t)  is  a  bounded  function. 


(F41) 


Edge  crack 

For  a  =  0,  we  obtain  an  edge  crack.  This  case  is  similar  to  the  case  discussed  in 
Section  3.3.4.  The  singularities  obtained  are  as  follows: 

Qi  =  0,/3i  =  J.  (F42) 

We  can  now  express  the  unknown  function  as 
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9(t)  = 


m 


(b  -  t)l 

where  again  g*(t)  is  a  bounded  function. 


(F43) 


Numerical  Procedure 

Internal  crack 

We  normalize  the  singular  integral  equation  in  (F36)  using 


b  —  cl  b  +  CL 

t  — - T  H - 

2  2 


(F44) 


b  —  a  b  “h  a 
x=~s+~- 

Thus  we  obtain, 


(F45) 


^J_1[ir-s)+KM\9^dr  -p(s)  4/io 

(F46) 

where 

a  b  +  a  b  —  a  b  +  a\ 

S+  2  ’  2  r+  2  )’ 

(F47) 

?(D=scrr+6r) 

1  , 

(F48) 

rv,  /  x  fb-a  b  +  a\ 

?W  =  P(  2  *  +  2  J 

I- 

(F49) 

Using  (F40),  (F43)  and  {F 44),  we  can  express  the  unknown  function  as  follows: 

~  /  \  F(r) 

9  (r)  =  (™) 
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where  F(r)  is  a  bounded  function.  It  may  be  expressed  as 


F(r)  =  X>„r„(r). 


(F51) 


Substituting  (F49)  and  (F50)  into  (F45)  we  obtain, 


Tn(r)  ,  ~  ,  s  «  +  1 


(F52) 


The  single  valuedness  condition  dictates  that  Ao  =  0.  Converting  equation  (F51)  into  a 
set  of  algebraic  equations  as  explained  in  Chapter  4  we  obtain, 


where  the  collocation  points  are  chosen  as  follows: 


(k  =  1.  JV) 
(F53) 


(2k  -  1)tt 


(fc  =  1.  JV). 


(F54) 


Edge  crack 


Normalizing  the  singular  integral  equation  in  (F36)  using 


;  b  b 

t  =  2r+r 


(F55) 


b  b 

x=r+2 ■ 

we  obtain, 


(F56) 


jr/,[(T^j  +  #(s.r)]sr  M*'  =P  (*) 


v  K  +  1 
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(F57) 


where 


(F58) 

(F59) 

p  W  =  p(j«  +  j)- 

(F60) 

Using  (jF42),  (F55)  and  (F56)  we  can  express  the  unknown  density  function  as  follows: 


g  (r)  =  -^L  (F61) 

where  F(r)  is  a  bounded  function.  It  may  be  expressed  as 

F(r)  =  Y,AnTn{r).  (F62) 

n= 0 

Substituting  (F60)  and  (F61)  into  (F56)  we  obtain, 


Ip-Clth 


(F63) 


We  can  convert  this  integral  equation  into  a  set  of  linear  algebraic  equations  as  discussed 
in  Chapter  4.  Thus  we  obtain 


- ['  [rJ-1  +K(sk,  r)]  dr  =  $  (sk) 

7-iUr-  sfc)  -t  VI  —  r  4a*o 


where  the  collocation  points  are  chosen  as  follows: 


(k  =  1..N ) 
(F64) 


Sk 


(2k  —  1)7T 
2(N  +  1) 


(k  —  1..N  +  1). 


(F65) 
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Numerical  Results 


Table  FA  provides  the  normalized  stress  intensity  factors  for  an  edge  crack  in  a 
nonhomogeneous  half-space  subjected  to  fixed  grip  loading.  The  results  are  evaluated  for 
various  values  of  the  normalized  nonhomogeneity  parameter  /3b  where  b  is  the  crack 
length.  These  results  have  also  been  plotted  in  Fig  FA. 
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h 

ay/b 

1.121523 

-  0.0001 

1.1215 
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qqH 

1.1139 
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1.1293 

0.4 

1.3324 

-0.4 

1.4103 

-0.5 

EEgQBI 

0.7 

-0.7 

0.8 

0.9 

1.2561 

Table  JF.l  Normalized  stress  intensity  factors  for  edge  crack  in 
nonhomogeneous  half  space  subjected  to  fixed  grip  loading. 

8flQ€o 

a  =  - - . 

1  +  K 
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Pb 

Figure  F. 2  Normalized  stress  intensity  factor  for  a  nonhomogeneous  half-plane  as  a 

function  of  normalized  nonhomogeneity  parameter  /3b  under  fixed  grip  loading. 
8/ioeo 

G  1  +  «* 
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